matlab第5章作业.docx
《matlab第5章作业.docx》由会员分享,可在线阅读,更多相关《matlab第5章作业.docx(10页珍藏版)》请在冰豆网上搜索。
matlab第5章作业
《信号与系统》课程研究性学习手册
姓名吕航(第十一组)
学号09213044
同组成员马琳
刘煜晨
朱志杰
王琛珏
指导教师魏杰
时间2010/11
专题研讨三、信号与系统的复频域分析
研讨题目:
M5-1
题目分析:
(1)分别取RC为0.00250.0250.00001时的值
(2)取RC为0.0025时画波形
(3)用sum求和。
仿真程序:
untiled1.m
RC=0.0025;
RC2=0.025;
RC3=0.00001;
b=[1];
a=[RC1];a2=[RC21];a3=[RC31];
w=0:
3000;
H=freqs(b,a,w);H2=freqs(b,a2,w);H3=freqs(b,a3,w);
subplot(2,2,1);
plot(w,abs(H));
subplot(2,2,2);
plot(w,abs(H2));
subplot(2,2,3);
plot(w,abs(H3));
untiled2.m
RC=0.0025;
b=[1];
a=[RC1];
w=0:
3000;
t=0:
0.00001:
0.1;
H=freqs(b,a,w);
x=cos(200*t)+cos(3000*t);
y=abs(H(200))*cos(200*t+angle(H(200)))
+abs(H(3000))*cos(3000*t+angle(H(3000)));
subplot(2,1,1);
plot(t,x);
subplot(2,1,2);
plot(t,y);
untiled3.m
RC=input('RC=');
t=linspace(0,0.2,1024);
w0=100*pi;
H0=(1/RC)/(1/RC+j*w0);
x=10*abs(sin(w0*t));
y=abs(H0)*10*abs(sin(w0*t));
plot(t,y);
ylabel('y(t)');
xlabel('ʱ¼ä£¨Ã룩');
title(['RC=',num2str(RC)])
untiled4.m
RC=0.01;
t=linspace(0,0.2,1024);
w0=100*pi;
H0=(1/RC)/(1/RC+j*w0);
x1=sum(10*abs(sin(w0*t)));
y1=sum(abs(H0)*10*abs(sin(w0*t)));
subplot(2,1,1);
plot(t,x1);
ylabel('x1(t)');
xlabel('ʱ¼ä£¨Ã룩');
subplot(2,1,2);
plot(t,y1);
ylabel('y1(t)');
xlabel('ʱ¼ä£¨Ã룩');
仿真结果:
图一
图二
图三
(1)
图三
(2)
图四
结果分析:
取不同的RC值可以对不同频率的波进行阻断。
可以从图中很明显的看出影响。
自主学习内容:
陈后金《信号与系统》高等教育出版社编著
研讨题目:
M5-2
题目分析:
模拟滤波器的作用,与理论值比较。
仿真程序:
untiled5.m
t=0:
0.005:
2.5;
x=(t>=0&t<=1)+t.*(t>=0&t<=1).*cos(50*t)+(2-t).*(t>1&t<=2).*cos(50*t);
plot(t,x);
untiled6.m
b=[10000];
a=[1,26.131,341.42,2613,10000];
[H,w]=freqz(b,a);
subplot(2,1,1);
plot(w,abs(H));
xlabel('\omega');
ylabel('|H(j\omega)|');
subplot(2,1,2);
plot(w,angle(H));
xlabel('\omega');
ylabel('\phi(\omega)');title('Magnituderesponse');
untiled7.m
sys=tf([10000],[126.131341.422613.110000]);
x=(t>=0&t<=1)+t.*(t>=0&t<=1).*cos(50*t)+(2-t).*(t>1&t<=2).*cos(50*t);
x1=x.*cos(50*t);
y1=lsim(sys,x,t);
y2=lsim(sys,x1,t);
subplot(2,1,1);
plot(t,y1);
title('y1(t)');
subplot(2,1,2);
plot(t,y2);
title('y2(t)');
仿真结果:
结果分析:
可以从图中看出,原信号和输出信号基本一致,输出信号在幅度上有所改变,在时间上有些延迟。
这说明这2个信号都是由低频分量组成的,它们的傅里叶变换的频率范围都大致位于给出的低通滤波器的通带之内。
自主学习内容:
陈后金《信号与系统》高等教育出版社编著
研讨题目:
M5-4
题目分析:
利用modulate函数调制幅度。
仿真程序:
Untiled8.m
Fm=1;Fc=20;Fs=1000;
N=1000;k=0:
N-1;
t=k/Fs;
x=sin(2.0*pi*Fm*t);
subplot(2,2,1);
plot(t(1:
1000),x(1:
1000));
xlabel('Time(s)');
title('ModulateSignal');
Xf=abs(fft(x,N));
subplot(2,2,2);stem(Xf(1:
1000));
xlabel('Frequency(Hz)');
title('ModulateSignal');
y=modulate(x,Fc,Fs);
subplot(2,2,3);
plot(t(1:
1000),y(1:
1000));
xlabel('Time(s)');
title('ModulateSignal(AM)');
Xf=abs(fft(y,N));
subplot(2,2,4);stem(Xf(1:
1000));
xlabel('Frequency(Hz)');
title('ModulatedSignal(AM)');
仿真结果:
结果分析:
调制后,频域上做相应的搬移。
当由于是利用离散代替连续,因为抽样,而使频谱发生了周期化。
自主学习内容:
陈后金《信号与系统》高等教育出版社