《信号与系统仿真完结版》由会员分享，可在线阅读，更多相关《信号与系统仿真完结版（26页珍藏版）》请在金锄头文库上搜索。
1、 信号与系统仿真武汉大学电气工程学院 09级7班徐碧川学号 2009302540225By Slytherin2011/12/7实验1：连续时间信号的表示及可视化；（分别取）；（分别画出不同周期个数的波形）；解：实验分析：本实验较为简单，但函数图象数目众多时不宜一个个编程实现，所以这里采用函数调用的形式一次性把所有图画出来，各函数如下：1、各函数如下：的函数编写如下：%冲击函数function y=delta(t)dt=0.01;y=(u(t)-u(t-dt)/dt;的函数编写如下：%阶跃函数function y=u(t)y=(t=0);的函数编写如下：%exp（2t）function y=y
2、31(t)y=exp(2*t);%exp(-2*t)function y=y32(t)y=exp(-2*t);的函数编写如下：%R4（t）function y=y4(t)y=u(t)-u(t-4);的函数编写如下：%Sa（2t）function y=y5(t)y=sinc(2*t/pi);的函数编写如下：%sin(2*t)function y=y6(t)y=sin(2*t);%sin(4*t)function y=y7(t)y=sin(4*t);%sin(8*t)function y=y8(t)y=sin(8*t);2、主程序如下：t=-5:0.001:5; %画图的区间plot(t,delt
3、a(t);title(单位冲击函数); %画单位冲击函数figure,plot(t,u(t);title(单位阶跃函数); %画单位阶跃函数figure,plot(t,y31(t);title(exp（2t）); %画exp（2t）figure,plot(t,y32(t);title(exp（-2t）); %画exp（-2t）figure,plot(t,y4(t);title(R（t）); %画R（t）figure,plot(t,y5(t);title(Sa（2t）); %画Sa（2t）figure,plot(t,y6(t);title(sin（2t）); %画sin（2t）figure,pl
4、ot(t,y7(t);title(sin(4*t); %画sin(4*t)figure,plot(t,y8(t);title(sin(4*t); %画sin(4*t)3、运行结果：（1）（2）（3）（4）（5）（6）实验心得：通过本次实验，使我在直观上了解了连续信号的特点。通过Matlab编程，使我了解到养成良好编程习惯的重要性，比如注释的习惯以及文件命名的注意事项。实验2：离散时间信号的表示及可视化；（分别取）；（分别取不同的N值）；（分别取不同的值）；解：实验分析：本实验是离散信号，只要对连续信号进行取样即可，所以之言控制时间轴为离散的点输出函数即可。函数图象数目众多时不宜一个个编程实现，
5、所以这里采用函数调用的形式一次性把所有图画出来，各函数如下：1、各函数如下：的函数编写如下：funcnion y=delna(n)y=(n=0);的函数编写如下：%阶跃函数funcnion y=u(n)y=(n=0);的函数编写如下：%exp（2n）funcnion y=y31(n)y=exp(2*n);%exp(-2*n)funcnion y=y32(n)y=exp(-2*n);的函数编写如下：%R4（n）funcnion y=y4(n)y=u(n)-u(n-4);的函数编写如下：%Sa（2n）funcnion y=y5(n)y=sinc(2*n/pi);的函数编写如下：%sin(2*n)f
6、uncnion y=y6(n)y=sin(2*n);%sin(4*n)funcnion y=y7(n)y=sin(4*n);%sin(8*n)funcnion y=y8(n)y=sin(8*n);2、主程序如下：n=-5:1:5; %画图的区间stem(n,delta(n);title(单位冲击函数); %画单位冲击函数figure,stem(n,u(n);title(单位阶跃函数); %画单位阶跃函数figure,stem(n,y31(n);title(exp（2n）); %画exp（2n）figure,stem(n,y32(n);title(exp（-2n）); %画exp（-2n）fig
7、ure,stem(n,y4(n);title(R（n）); %画R（n）figure,stem(n,y5(n);title(Sa（2n）); %画Sa（2n）figure,stem(n,y6(n);title(sin（2n）); %画sin（2n）figure,stem(n,y7(n);title(sin(4*n); %画sin(4*n)figure,stem(n,y8(n);title(sin(4*n); %画sin(4*n)3、运行结果：（1）（2）（3）（4）（5）（6）实验心得：通过本次实验，使我在直观上了解了离散信号的特点。离散信号可以通过对连续信号的采样得到。通过Matlab编程，
8、使我了解到养成良好编程习惯的重要性，比如注释的习惯以及文件命名的注意事项。实验3： 系统的时域求解1. 设，求：，并画出、波形。2. 求因果线性移不变系统的单位抽样响应，并绘出的幅频及相频特性曲线。解:先把两个函数编写成M函数文件待调用，对于第一问，可以直接用步长为1的conv函数进行卷积的求解，并作图。对于第二问，可以用impz函数求单位样值相应，并画出其幅频特性和相频特性。1、的函数编写如下：function y=un(n)y=(n=0);function y=xn(n)y=un(n)-un(n-10);function y=hn(n)y=power(0.9,n).*un(n);下面编写第
9、一问的主程序：n=0:10;yn=1*conv(xn(n),hn(n); %步长为1的卷积subplot(3,1,1),stem(n,xn(n),xlabel(x(n);subplot(3,1,2),stem(n,hn(n),xlabel(h(n);subplot(3,1,3),stem(yn),xlabel(y(n)=x(n)*h(n);运行结果：2、主程序编写如下：%实验3的第二问%实验3的第二问A=1,0,-0.81;B=1,0,-1;n=-10:30;%figure,impz(A,B,n);h=impz(A,B,n);figure,stem(n,h),title(单位样值响应);fig
10、ure,freqs(B,A);运行结果：实验心得：Matlab可以方便地求解卷积以及系统的响应，是求解信号与系统问题的有力工具，通过本次实验，让我对离散系统的响应有了直观的理解，同时让我进一步提高了Matlab编程的水平。实验4： 信号的DFT分析计算余弦序列的DFT。分别对N=10、16、22时计算DFT，绘出幅频特性曲线，分析是否有差别及产生差别的原因。解：实验分析：通常计算机以DFT来描述信号的频谱。而求解DFT的运算量往往很大，而FFT是DFT的快速算法，使运算大幅度减小。Matlab自带了求信号的FFT的函数，可以方便地求解信号的快速傅里叶变换。这里通过求解同一个信号的DFT和FFT
11、进行对比。1、 函数的编写：由于Matlab没有自带的求DFT的函数，所以这里根据DFT的定义编写求解信号DFT的函数如下：%求DFT的函数function WNnk=dft(xn,N)WN=exp(-j*2*pi/N); %旋转因子for k=0:N-1 for n=0:N-1 WNnk1(n+1)=xn(n+1)*WN.(n*k); end WNnk(k+1)=sum(WNnk1(1:N);endMatlab内置了求FFT的函数，即函数fft；2、 主程序的编写：clearN=10;n=0:N-1;x=cos(pi./8)*n); y=dft(x,N);yy=fft(x); %分别求DFT
12、和FFTsubplot(2,2,1),stem(n,y),title(N=10的DFT频谱);subplot(2,2,3),stem(n,yy,r),title(N=10的FFT频谱);subplot(2,2,2),stem(n,abs(y) ,title(N=10的DFT幅频特性);subplot(2,2,4),stem(n,abs(fft(x),r),title(N=10的FFT幅频特性);N=16;n=0:N-1;x=cos(pi./8)*n); y=dft(x,N);yy=fft(x); %分别求DFT和FFTfigure,subplot(2,2,1),stem(n,y),title(
13、N=16的DFT频谱);subplot(2,2,3),stem(n,yy,r),title(N=16的FFT频谱);subplot(2,2,2),stem(n,abs(y) ,title(N=16的DFT幅频特性);subplot(2,2,4),stem(n,abs(fft(x),r),title(N=16的FFT幅频特性);N=22;n=0:N-1;x=cos(pi./8)*n); y=dft(x,N);yy=fft(x); %分别求DFT和FFTfigure,subplot(2,2,1),stem(n,y),title(N=22的DFT频谱);subplot(2,2,3),stem(n,yy,r),title(N=22的FFT频谱);subplot(2,2,2),stem(n,abs(y) ,title(N=22的DFT幅频特性
