中南大学信号与系统实验报告.docVIP

实验一 基本信号的生成 ⑴ 运行以上九个例子程序，掌握一些常用基本信号的特点及其MATLAB实现方法；改变有关参数，进一步观察信号波形的变化。 ⑵ 在 范围内产生并画出以下信号： ； ； ； 。 答：a） k=-10:10; delta=[(k-0)==0]; stem(k,delta) b) k=-10:10; delta=[(k+2)==0]; stem(k,delta) c) k=-10:10; delta=[(k-4)==0]; stem(k,delta) d) k=-10:10; delta1=[(k+2)==0]; delta2=[(k-4)==0]; stem(k,2*delta1-delta2) ⑶ 在 范围内产生并画出以下信号： ； ； 。 请问这三个信号的基波周期分别是多少？ 答：a) k=0:31; fk1=sin(pi/4*k); fk2=cos(pi/4*k); stem(k,fk1.*fk2) b) k=0:31; fk=cos(pi/4*k).*cos(pi/4*k); stem(k,fk) c) k=0:31; fk=sin(pi/4*k).*cos(pi/8*k); stem(k,fk) 实验二 信号的基本运算 ⑴ 运行以上三个例题程序，掌握信号基本运算的MATLAB实现方法；改变有关参数，考察相应信号运算结果的变化特点与规律。 ⑵ 已知信号如下图所示： 用MATLAB编程复现上图； 画出的波形； 画出的波形； 画出的波形。 答：a) t=-4:0.001:4; ft=tripuls(t,6,0); plot(t,2*ft) b) t=-4:0.001:4; ft=tripuls((2-2*t),6,0); plot(t,2*ft) c)函数名为functri（相应的.m文件名为functri.m），程序如下： function yt=functri(t) yt=2*tripuls(t,6,0); 然后利用diff函数 h=0.001;t=-4:h:4; y1=diff(functri(t))*1/h; plot(t(1:length(t)-1),y1) title(df(t)/dt) d)然后利用quad函数 t=-4:0.1:4; for x=1:length(t) y2(x)=quad(functri,-3,t(x)); end plot(t,y2) title(integral of f(t)) 实验三 系统的时域分析 ⑴ 运行以上五个例题程序，掌握求解系统响应的MATLAB分析方法；改变模型参数，考察系统响应的变化特点与规律。 ⑵ 设离散系统可由下列差分方程表示： 计算时的系统冲激响应。 答：k=-20:100 a=[1 -1 0.9]; b=[1]; h=impz(b,a,k); subplot(2,1,1) stem(k,h) ⑶设，输入，求系统输出。 （取） 答： k=-10:50;A=1;a=0.9; hk1=A*a.^k; uk1=[k=0]; hk=hk1.*uk1; uk2=[k-10=0]; fk=uk1-uk2; z=conv(fk,hk); N=length(z); stem(0:N-1,z) ⑷ 已知滤波器的传递函数： 输入信号为为随机信号。试绘出滤波器的输出信号波形。（取） 答： R=101; d=rand(1,R)-0.5; t=0:100; s=2*sin(0.05*pi*t); f=s+d; subplot(2,1,1); plot(t,d,g-,t,s,b--,t,f,r-); xlabel(Time index t); legend(d[t],s[t],f[t]); %d(t)即为随机信号ω(t) title(处理前的波形) b=[0.22 0];a=[1 -0.8]; y=filter(b,a,f); subplot(2,1,2); plot(t,s,b--,t,y,r-); xlabel(Time index t); legend(s[t],y[t]); title(滤波器输出波形) 实验四 周期信号的频域分析 1.仿照例程，实现下述周期信号的傅立叶级数分解与合成： 1 1 -3 -4 5 4 1 O 要求： （a）首先，推导出求解，，的公式，计算出前10次系数； （b）利用MATLAB求解，，的值，其中，求解前10次系数，并给出利用这些系数合成的信号波形。 答：(a) (b)源代码如下 function [A_sym,B_sym]=CTFShchsym2 % 采用符号计算求一个周期内连续时间函数f的三角级