信号研究学习5.doc

信号与系统课程研究性学习手册姓名 任叶 学号 09212049 同组成员 曹冰清 09212032 陈金 09212033 指导教师 魏杰 时间 2010年10月21日 专题研讨二、信号与系统的频域分析研讨题目：1 求M5-1电路（1）不同RC值的幅度响应（2）当x（t）=cos（100t）+cos（3000t）时，确定合适的RC值滤去其中的高频信号，并画出图形。（3）求x（t）=10|sin（100t），取不同值，计算并画出x（t）通过该系统的响应y（t）。2 （1）计算信号x（t）=x1（t）+x2（t）cos(50t)的波形。（2）画出模拟M5-2中的H（jw）（3）用lsim函数求出信号x（t）和x（t）cos（50t）通过系统的响应y1和y2。3 调制信号为1Hz，载波信号频率为20Hz，分析比较信号调制前后的频谱。题目分析：1（1）分别取RC为0.0025 0.025 0.00001时的值 （2）取RC为0.0025时画波形 （3）同上，用sum（x，2）求和。 Y（t）=|(H(314)|*10|*(cos(100*pi*t2+angle(H(314)| 仿真程序：1（1）RC=0.0025;RC2=0.025;RC3=0.00001;b=1;a=RC 1;a2=RC2 1;a3=RC3 1;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); (2)RC=0.0025;b=1;a=RC 1;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); (3)RC=1; （在此处改变不同的RC值即可）b=1;a=RC 1;w=0:3000;t=0:0.00001:0.1;H=freqs(b,a,w);x=abs(10*cos(100*pi*t);y=abs(H(314)*10*abs(cos(100*pi*t+angle(H(314);subplot(2,2,1);plot(t,x);subplot(2,2,2);plot(t,y);t1=linspace(0,0.02,100);x1=0.0002*abs(10*cos(100*pi*t1);s=sum(x1,2)/0.02;subplot(2,2,3);stem(1,s);t2=linspace(0,0.02,100);x2=0.0002*abs(H(314)*10*abs(cos(100*pi*t2+angle(H(314);s2=sum(x2,2)/0.02;subplot(2,2,4);stem(1,s2);RC为0.0025时的情况：RC为1时的情况：RC为0.00001时的情况：2t=0:0.005:2.5;x=(t=0&t=0&t1&t=2).*cos(50*t);x1=x.*cos(50*t);subplot(2,2,1);plot(t,x);b=10000;a=1 26.131 341.42 2613.1 10000;w=0:30;H=freqs(b,a,w);subplot(2,2,2);plot(w,abs(H);subplot(2,2,3);plot(w,angle(H);sys=tf(10000,1 26.131 341.42 2613.1 10000);y1=lsim(sys,x1,t); （在此处修改x即可得到不同波形）subplot(2,2,4);plot(t,y1);3Fm=1;Fc=10;Fs=1000;N=1000;k=0:N-1;t=k/Fs;x=sin(2*pi*Fm*t);subplot(2,2,1);plot(t,x);axis(0,1,-1,1);xlabel(Time(s);title(Modulate Signal);Xf=abs(fft(x,N);subplot(2,2,2);stem(Xf(1:999);xlabel(Frequency(Hz);y=modulate(x,Fc,Fs,am);subplot(2,2,3);plot(t,y);xlabel(Time(s);axis(0,1,-1,1);title(Modulated Signal(AM);Xf=abs(fft(y,N);subplot(2,2,4);stem(Xf(1:999);xlabel(Frequency(Hz);title(Modulated Signal(AM);仿真结果：见代码旁结果分析：1. 取不同的RC值可以对不同频率的波进行阻断。2. 根据低通高阻原理x（t）中的x1（t）为一低频信号，故应该通过x（t）为一高频信号，应该被阻断，理想下应为一个方波，但由于带宽问题，所以有部分失真。同理x（t）cos（50t）应为一个三角波。3 调制后，频域上做相应的搬移。当由于是利用离散代替连续，似乎因为抽样，而使频谱发生了周期化。自主学习内容： Matlab部分运算语句，fft语句的使用,function的使用 阅读文献： Matlab7.0从入门到精通 人民邮电出版社 求是科技 编著 信号与系统 高等教育出版社 陈后金 编著 精通GUI图形界面编程 北京大学出版社 施晓红 周佳 编著发现问题(专题研讨或相关知识点学习中发现的问题)： 对.*这一运算