MatLab: сложить графики, один из которых составлен через цикл. цифровой сигнал складываем с аналоговым

@remi · Регистрация: 21.04.2015

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

Здравствуйте! у меня есть задание по тэц - есть сигнал и к нему нужно просто прибавить помеху. Это цифровой сигнал, он идет прямоугольниками

Matlab M

>> % Построение цифрового сигнала
>> % Код из 8 символов
>> K = [0 0 1 1 0 0 1 1];
>> % Уровень 1
>> U1=5;
>> % Уровень 0
>> U2 = -5;
>> % Длительность элементарного символа сигнала
>> dt = 0.1;
>> % Период сигнала
>> T=10*dt;
>> S1 = K*U1;
>> S2 = S1+U2;
>> S = S1+S2;
>> % Временной интервал
>> t = linspace(0,T, 200);
>> for i = 1:200
if t(i)<dt
u(i)=S(1);
end
if t(i)>=dt & t(i)<2*dt
u(i)=S(2);
end
f t(i)>=2*dt & t(i)<3*dt
u(i)=S(3);
end
if t(i)>=3*dt & t(i)<4*dt
u(i)=S(4);
end
if t(i)>=4*dt & t(i)<5*dt
u(i)=S(5);
end
if t(i)>=5*dt & t(i)<6*dt
u(i)=S(6);
end
if t(i)>=6*dt & t(i)<7*dt
u(i)=S(7);
end
if t(i)>=7*dt & t(i)<8*dt
u(i)=S(8);
end
if t(i)>8*dt
u(i)=0;
end
end
>> plot(t*10,u)
>> grid on

вот так мне помогли построить цифровой сигнал

это график помехи

Matlab M

>> syms t1
>> n=14;     'не обращайте на эти строчки внимания, это мой вариант, с помощью этого я задаю промежуток .
>> f1=(15+7*log10(n+1))*10^3; '
>> T1=1/f1;
>> wx = 101585.9080;
>> x1=@(t1)1.5*cos(wx.*t1); 
>> t1=0:0.000001:10*T1; ' вот этот промежуток
>> X1 = x1(t1);
>> plot(t1,X1)

в общем, как я понимаю, та программа, с помощью которой я строю цифровой сигнал, она сразу выстраивает линию из linespace. и я даже начинаю подумывать, что сложить эти графики невозможно. Попытки всунуть в эти строчки построения цифрового сигнала формулу графика помехи ничем хорошим не увенчались - матлаб игнорировал и строил его как обычно. Если вы сможете помочь, можете даже соединить обычный график с этим цифровым сигналом - дальше я разберусь.

@Nick07 · 22.04.2015, 17:25

plot(t*10,u, t1,X1)

Добавлено через 1 минуту

Сообщение от remi

X1 = x1(t1);

А тут МL не ругается?

Добавлено через 5 минут

Сообщение от remi

матлаб игнорировал и строил его как обычно

MATLAB все правильно делает, ошибка у Вас. Посмотрите значения по оси Х, они отличаются на 4 порядка! Растяните ось Х в самом начале и увидите сигнал шума

@remi · 25.04.2015, 17:29 **[ТС]**

нет, он там не ругается
нет, если сделать plot, (ну и выбрать правильный диапазон, как вы сказали

)то эти графики просто нарисуются вместе, а нужно их сложить
и я должна была вообще показать, как именно мне нужно их сложить

ну в общем, все равно спасибо большое!=) я немного попозже догадалась, как сделать)
вот так:

Matlab M

% Построение цифрового сигнала
% Код из 8 символов
 K = [0 0 1 1 0 0 1 1];
% Уровень 1
 U1=5;
 % Уровень 0
 U2 = -5;
 % Длительность элементарного символа сигнала
 dt = 0.1;
 % Период сигнала
 T=10*dt;
 S1 = K*U1;
 S2 = S1+U2;
 S = S1+S2;
 % Временной интервал
 t = linspace(0,T, 200);
for i = 1:200
if t(i)<dt
u(i)=S(1);
end
if t(i)>=dt & t(i)<2*dt
u(i)=S(2);
end
if t(i)>=2*dt && t(i)<3*dt
u(i)=S(3);
end
if t(i)>=3*dt & t(i)<4*dt
u(i)=S(4);
end
if t(i)>=4*dt & t(i)<5*dt
u(i)=S(5);
end
if t(i)>=5*dt & t(i)<6*dt
u(i)=S(6);
end
if t(i)>=6*dt & t(i)<7*dt
u(i)=S(7);
end
if t(i)>=7*dt & t(i)<8*dt
u(i)=S(8);
end
if t(i)>8*dt
u(i)=0;
end
end
 % Зададим параметры для аналогового сигнала и помехи
 n=14;
 % Аналоговый сигнал
 n=14;
 f1=(15+7*log10(n+1))*10^3;
omega1=2*pi*f1;
% Помеха
wx = 101585.9080;
% Временной интервал
t = linspace(0,T, 200);
for i = 1:200
if t(i)<dt
a(i)=S(1)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
if t(i)>=dt & t(i)<2*dt
a(i)=S(2)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
if t(i)>=2*dt & t(i)<3*dt
a(i)=S(3)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i)+ 1.5*cos(wx*i);
end
if t(i)>=3*dt & t(i)<4*dt
a(i)=S(4)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
if t(i)>=4*dt & t(i)<5*dt
a(i)=S(5)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
if t(i)>=5*dt & t(i)<6*dt
a(i)=S(6)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
if t(i)>=6*dt & t(i)<7*dt
a(i)=S(7)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i)+ 1.5*cos(wx*i);
end
if t(i)>=7*dt & t(i)<8*dt
a(i)=S(8)+(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
if t(i)>8*dt
a(i)=(n/9)*cos((1.2*omega1)*i)+(n/7)*cos((1.3*omega1)*i)-(n/4)*cos(omega1*i) + 1.5*cos(wx*i);
end
end
 plot(t*10,u,t*10,a)
 grid on

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	22.04.2015, 17:25	2
	plot(t10,u, t1,X1) Добавлено через 1 минуту* Сообщение от remi X1 = x1(t1); А тут МL не ругается? Добавлено через 5 минут Сообщение от remi матлаб игнорировал и строил его как обычно MATLAB все правильно делает, ошибка у Вас. Посмотрите значения по оси Х, они отличаются на 4 порядка! Растяните ось Х в самом начале и увидите сигнал шума 0

	25.04.2015, 17:29