Моделирование прохождения сигнала - Matlab - Ответ 7307377

@Олеся7777 · Регистрация: 06.03.2015

Author24 — интернет-сервис помощи студентам

нужно найти ошибку

Matlab M

clc;
clear;
N=60;%Количество отсчетов
R=32;%Скорость передачи в кбит/с
fs=8000;%частота дискретизации
f=[300 660 900 330];%частоты, входящие в сигнал
A=[4 7 10 1];%Амплитуды
t=0:1/fs:N/fs;
K=4095/(A(1)+A(2)+A(3)+A(4));%Нормирование амплитуд
Si=K*(A(1)*cos(2*pi*f(1)*t)+A(2)*cos(2*pi*f(2)*t)+A(3)*cos(2*pi*f(3)*t)+A(4)*cos(2*pi*f(4)*t));
%Начальные установки 
Se(1)=0;
y(1)=0;
b1(1)=1;
b2(1)=1;
b3(1)=1;
b4(1)=1;
b5(1)=1;
b6(1)=1;
td(1)=1;
tr(1)=1;
al(1)=0;
a1(1)=0;
a2(1)=0;
ap(1)=0;
dml(1)=0;
dms(1)=0;
W16=[-1.38 27.44];
F16=[0 7];
W24=[-0.25 1.88 8.56 36.38];
F24=[0 1 2 7];
W32=[-0.75 1.13 2.56 4.00 7.00 12.38 22.19 70.13];
F32=[0 0 0 1 1 1 3 7];
for k=1:N  
d(k)=Si(k)-Se(k);%Расчет разности сигналов
if k>1
   y(k)=al(k)*yu(k-1)+[1-al(k)]*yl(k-1);%Расчет шага квантования
end
%определение квантованного значения исходя из скорости передачи.
Q=log2(abs(d(k)))-y(k);
if R==32
   if q>=3.12
       I(k)=7;
       dq(k)=2^(3.32+y(k))*sign(d(k));
   elseif q>=2.72
       I(k)=6;
       dq(k)=2^(2.91+y(k))*sign(d(k));
   elseif q>=2.34
       I(k)=5;
       dq(k)=2^(2.52+y(k))*sign(d(k));
   elseif q>=1.91
       I(k)=4;
       dq(k)=2^(2.13+y(k))*sign(d(k));
   elseif q>=1.38
       I(k)=3;
       dq(k)=2^(1.66+y(k))*sign(d(k));
   elseif q>=0.62
       I(k)=2;
       dq(k)=2^(1.05+y(k))*sign(d(k));
   elseif q>=-0.98
       I(k)=1;
       dq(k)=2^(0.031+y(k))*sign(d(k));
   elseif q<-0.98
       I(k)=0;
       dq(k)=0;
   end
   yu(k) = (1-2^(-5))*y(k)+2^(-5)*W32(I(k)+1);
   if k>1
       yl(k) = (1-2^(-6))*yl(k-1)+2^(-6)*yu(k);
   else
       yl(k) = 2^(-6)*yu(k);
   end
   dms(k+1) = (1-2^(-5))*dms(k) + 2^(-5)*F32(I(k)+1);
   dml(k+1) = (1-2^(-7))*dml(k) + 2^(-7)*F32(I(k)+1);
end    
if or(abs(dms(k)+dml(k))>=2^(-3)*dml(k),or(y(k)<3,td(k)==1))
   ap(k+1)=(1-2^(-4))*ap(k)+2^(-3);
elseif tr==1
   ap(k+1)=1;
else
   ap(k+1)=(1-2^(-4))*ap(k);
end
if ap(k)>=1
   al(k+1)=1;
else
   al(k+1)=ap(k);
end
%расчет нулей предсказателя
if k>7
   b6(k)=(1-2^(-8))*b6(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 6));
   b5(k)=(1-2^(-8))*b5(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 5));
   b4(k)=(1-2^(-8))*b4(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 4));
   b3(k)=(1-2^(-8))*b3(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 3));
   b2(k)=(1-2^(-8))*b2(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 2));
   b1(k)=(1-2^(-8))*b1(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 1));
elseif k>6
   b6(k)=2^(-7)*sign(dq(k))*sign(dq(k - 6));
   b5(k)=(1-2^(-8))*b5(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 5));
   b4(k)=(1-2^(-8))*b4(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 4));
   b3(k)=(1-2^(-8))*b3(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 3));
   b2(k)=(1-2^(-8))*b2(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 2));
   b1(k)=(1-2^(-8))*b1(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 1));
elseif k>5
   b6(k)=0;
   b5(k)=2^(-7)*sign(dq(k))*sign(dq(k - 5));
   b4(k)=(1-2^(-8))*b4(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 4));
   b3(k)=(1-2^(-8))*b3(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 3));
   b2(k)=(1-2^(-8))*b2(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 2));
   b1(k)=(1-2^(-8))*b1(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 1));
elseif k>4
   b6(k)=0;
   b5(k)=0;
   b4(k)=2^(-7)*sign(dq(k))*sign(dq(k - 4));
   b3(k)=(1-2^(-8))*b3(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 3));
   b2(k)=(1-2^(-8))*b2(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 2));
   b1(k)=(1-2^(-8))*b1(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 1));
elseif k>3
   b6(k)=0;
   b5(k)=0;
   b4(k)=0;
   b3(k)=2^(-7)*sign(dq(k))*sign(dq(k - 3));
   b2(k)=(1-2^(-8))*b2(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 2));
   b1(k)=(1-2^(-8))*b1(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 1));
elseif k>2
   b6(k)=0;
   b5(k)=0;
   b4(k)=0;
   b3(k)=0;
   b2(k)=2^(-7)*sign(dq(k))*sign(dq(k - 2));
   b1(k)=(1-2^(-8))*b1(k-1)+2^(-7)*sign(dq(k))*sign(dq(k - 1));
elseif k>1
   b6(k)=0;
   b5(k)=0;
   b4(k)=0;
   b3(k)=0;
   b2(k)=0;
   b1(k)=2^(-7)*sign(dq(k))*sign(dq(k - 1));    
end    
if k>6
   Sez(k)=b1(k-1)*dq(k-1)+b2(k-1)*dq(k-2)+b3(k-1)*dq(k-3)+b4(k-1)*dq(k-4)+b5(k-1)*dq(k-5)+b6(k-1)*dq(k-6);
elseif k>5
   Sez(k)=b1(k-1)*dq(k-1)+b2(k-1)*dq(k-2)+b3(k-1)*dq(k-3)+b4(k-1)*dq(k-4)+b5(k-1)*dq(k-5);
elseif k>4
   Sez(k)=b1(k-1)*dq(k-1)+b2(k-1)*dq(k-2)+b3(k-1)*dq(k-3)+b4(k-1)*dq(k-4);
elseif k>3
   Sez(k)=b1(k-1)*dq(k-1)+b2(k-1)*dq(k-2)+b3(k-1)*dq(k-3);
elseif k>2
   Sez(k)=b1(k-1)*dq(k-1)+b2(k-1)*dq(k-2);
elseif k>1
   Sez(k)=b1(k-1)*dq(k-1);
elseif k==1
   Sez(k)=0;
end
p(k)=Sez(k)+dq(k);
%Расчет полюсов предсказателя
if k>1
   a1(k)=(1-2^(-8))*a1(k-1)+(3*2^(-8))*sign(p(k))*sign(p(k-1));    
end
if k>1
   if abs(a1(k-1))<=0.5
       f=4*a1(k);
   else
       f=2*sign(a1(k-1));
   end
end
if k>2
   a2(k)=(1-2^(-7))*a2(k-1)+2^(-7)*(sign(p(k))*sign(p(k-2))-f*sign(p(k))*sign(p(k-1)));
elseif k>1  
   a2(k)=(1-2^(-7))*a2(k-1)+2^(-7)*(0-f*sign(p(k))*sign(p(k-1)));
end
if k<2
   Se(k+1)=a1(k)*(Se(k)+dq(k))+Sez(k);
else
   Se(k+1)=a1(k)*(Se(k)+dq(k))+a2(k)*(Se(k-1)+dq(k-1))+Sez(k);
end
if a2(k)<-0.71875
   td(k+1)=1;
else
   td(k+1)=0;
end
if and(a2(k)<-0.71875,abs(dq(k))>24*2^yl(k))
   tr(k+1)=1;
else
   tr(k+1)=0;
end
end
for i=1:N
I(i)=I(i)*sign(d(i));%Учитываем знак изменения
end;
for i=1:N
Sr(i)=Se(i)+dq(i);%Восстановление сигнала
end
subplot(2,2,1);
plot(Si(1:N));
grid on;
title('Входной аналоговый сигнал');
subplot(2,2,3);
plot(Sr(1:N));
grid on;
title('Выходной декодированный сигнал');
subplot(2,2,2);
stem(I(1:N));
grid on;
title('Квантованный разностный сигнал в канале связи');
subplot(2,2,4);
plot(y(1:N));
grid on;
title('Коэффициент адаптации шага');

Вернуться к обсуждению:
Моделирование прохождения сигнала Matlab

Следующий ответ

Другие темы раздела
Matlab "свернуть" все возможные сегменты кода Здравствуйте! Напротив операторов цикла for и while, а так же напротив function слева в редакторе м-функций есть кнопки + или - чтобы свернуть или развернуть кусок кода, который входит в цикл или в функцию. Вопрос - есть какая-то опция, позволяющая развернуть и свернуть все такие участки кода? https://www.cyberforum.ru/ matlab/ thread1389018.html	Matlab Отделение корней с использованием графической оценки Плоскости располагаются вертикально. Во всех плоскостях определить прямоугольную систему координат. Построить графическое отображение по заданным нелинейным уравнениям.
Matlab Расчет свертки. Undefined function 'mtimes' for input arguments of type 'cell' Выдает ошибку не понимаю почему function y = conv3 (u, v) m = length (v); n = length (u); y=zeros(1,m+n-1); for k=1:m for i= 1:n a = k+i-1; https://www.cyberforum.ru/ matlab/ thread1388954.html	Matlab Построение кохлеоиды в полярных координатах(уравнение) https://www.cyberforum.ru/ matlab/ thread1388937.html Построить кохлеоиду, уравнение которой в полярных координатах. r=(asinFi)/Fi, где Fi принадлежитa=1; f = 0:0.1:1000; r=((asin(f))/f); subplot(1,2,1) polar(f,r,'r') При выполнении необходимо создать две координатные плоскости. Плоскости располагаются горизонтально. В первой плоскости определить прямоугольную систему координат, а во второй полярную. Построить заданную линию в двух...
*Matlab Прошу помощи в построении графиков уравнений подобного типа: x+xy/((1+x)(1+y)) - y=0* Доброго всем времени суток! Прошу помощи в построении графиков уравнений подобного типа: x+xy/((1+x)(1+y)) - y=0 Вместо единиц могут быть и другие значения, но мне важна сама форма графика. Я пересмотрела форум, но только еще больше запуталась, к сожалению, я не математик. Буду рада, если кто-нибудь поможет, так как нелинейных уравнений у меня еще много, и мне хотелось бы на примере понять...	Matlab Построить трехмерный объект со сглаженной поверхностью https://www.cyberforum.ru/ matlab/ thread1388603.html Построить фигуру по данным формулам x = cos(u/c) * cos(u/b) * (a + cos(v)) + sin(u/b) * sin(v) * cos(v) y = sin(u/c) * cos(u/b) * (a + cos(v)) + sin(u/b) * sin(v) * cos(v) z = -sin(u/b) * (a + cos(v)) + cos(u/b) * sin(v) * cos(v) где 0 <= u <= 2bcpi и 0 <= v <= 4pi Заранее спасибо!
Matlab Как отключить Warrings found в редакторе m-функций? Объединил несколько функций в одну получилось семь тысяч строк, теперь редактировать в редакторе матлаб не получается, потому что после каждого изменения кода запускается поиск ошибок, и все виснет намертво на несколько секунд. Так невозможно нормально работать. Как с этим бороться подскажите! (копировать в блокнот, там что-то делать потом обратно - не лучший вариант, наверное можно как-то... https://www.cyberforum.ru/ matlab/ thread1388462.html	График поверхности, Matlab Нужно построить поверхность z^2=xy. С этим я разобрался, вроде, получился примерно такой код. = meshgrid(-10:0.1:10); F = z.^2 - x.y; hpatch = patch(isosurface(x,y,z,F,0)); isonormals(x,y,z,F,hpatch); set(hpatch,'FaceColor','g','EdgeColor','none'); camlight left; lighting phong;
Matlab Алгоритм, что использовать? (рекурсия/индексы?) Добрый день! Вопрос следующий. Пусть есть n одинаковых векторов с действительными числами размерностью (1,m). Также существует n-1 выражений, которые нужно посчитать, используя вектора a,b,...,n, перебирая все возможные кобинации компонент и сравнить с нулем. Например, дано два вектора: a=-1:1:1; b=-1:1:1; Выражение, которое необходимо рассчитать в этом случае a+b (в результате всех... https://www.cyberforum.ru/ matlab/ thread1388422.html	Matlab Как правильно пользоваться функцией lqr или Откуда брать матрицы Q,R,N ? Добрый день! Функция lqr в математическом пакете Matlab осуществляет синтез оптимального регулятора для непрерывной системы. Функция lqr возвращает матрицу K- оптимальных коэффициентов усиления. Для исполнения функции lqr ей необходимо передать соответствующие матрицы: A,B,Q,R,N Q,R - это весовые матрицы, которые отвечают за "штрафование" системы за использование энергии управления и... https://www.cyberforum.ru/ matlab/ thread1388136.html
Matlab Интегрирование в Матлаб. Найти площадь под кривой Здравствуйте, ребята! Появилась такая сложная задача, вроде бы не очень сложная, а на самом деле оказалась сложной задачей. Есть два массива, первый Func c размером 1 x 50 000, и другой E с размером 1 x 50000; при plot(E,Func) видите график, точнее пик такой, мне нужна площадь под этой кривой, пробовал разными методами, все возвращают NAN; Что делать? ======== Спасибо.	Matlab Вычислить значения функции на отрезке в заданном числе N равномерно отстоящих друг от друга точек https://www.cyberforum.ru/ matlab/ thread1388043.html Помогите пожалуйста

	06.03.2015, 14:22