tiaozhi.m

信号与系统时域、频域分析及MATLAB软件的应用及仿真m程序

R=0.005;
t=-1.2:R:1.2;
f=Heaviside(t+1)-Heaviside(t-1);
f1=f.*cos(10*pi*t);%已调信号
subplot(221)
plot(t,f);
xlabel('t')
ylabel('f(t)')
subplot(222);
plot(t,f1)
xlabel('t')
ylabel('f1(t)=f(t)*cos(10*pi*t)');
W1=40;
N=1000;
k=-N:N;
W=k*W1/N;
F=f*exp(-j*t'*W)*R;%求F(jw)
F=real(F);
F1=f1*exp(-j*t'*W)*R;%求F(jw)
F1=real(F1);
subplot(223);
plot(W,F);
xlabel('W')
ylabel('F(jw)');
subplot(224);
plot(W,F1);
xlabel('w');
ylabel('F1(jw)');