信号上机实验三.docVIP

信号上机实验三 实验目的 掌握信号时域运算的MATLAB实现方法。 掌握信号时域变换的MATLAB实现方法 实验原理 信号的时域运算包括信号的相加、相乘，信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换 实验内容 将实验原理中提到的例子在计算机上全部练习一遍； 解: 例3-2. 求下列两序列的和序列： f1(k)={-2,-1,0,1,2}, f2(k)={1,1,1,} 其MATLAB命令为 function[f,k]=isxj(f1,f2,k1,k2) k=min(min(k1),min(k2)):max(max(k1),max(k2)); s1=zeros(1,length(k));s2=s1; s1(find((k=min(k1))(k=max(k1))==1))=f1; s2(find((k=min(k2))(k=max(k2))==1))=f2; f=s1+s2; stem(k,f,filled) axis([(min(min(k1),min(k2))-1),(max(max(k1),max(k2))+1),(min(f)-0.5),(max(f)+0.5)]) 例3-3：己知f(k)=2k -3≤k≤3，画出f(-k)波形。 function[f,k]=isfz(f1,k1) f=fliplr(f1);k=-fliplr(k1); stem(k,f,filled) axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]) 例3-4. 己知f(k)=k2 -4≤k≤4，画出f(k-2)波形。 function[f,k]=isyw(ff,kk,k0) k=kk+k0;f=ff; stem(k,f,filled) axis([min(k)-1,max(k)+1,min(f)-0.5,max(f)+0.5]) 已知信号，画出的波形； 解; 3、已知两个连续信号，用MATLAB绘出下列信号的波形： 1） 2） 3） 4） 解; 4、离散序列，用MATLAB绘出下列序列的波形。 1） 2） 3） 4） 解：（1） （2） （3） (4) 实验思考题 在对信号进行平移、反折和尺度变换时，运算顺序对结果是否有影响？在运算中应该注意什么？ 解：运算顺序对结果没影响。为了画f（-t-t0）或f（-k-k0），最好是先平移，将f(t)平移为f(t-t0)或 f(k)平移为f(k-k0),然后再反转，将变量t和k相应地换为-t或-k。如果反转后再进行平移，由于这时自变量变为-t或-k，故平移方向与前述相反。 信号上机实验五 实验目的 1.用MATLAB实现周期信号傅里叶级数分解与综合。 2.利用MATLAB对周期信号进行频谱分析 实验原理 任何满足狄里赫里条件的周期信号，都可以表示成三角形式或指数形式的傅里叶级数展开式。 （5.1） （5.2） 一般来说，傅里叶级数有无限个非零值，即任何具有有限个间断点的周期信号都一定有一个无限非零系数的傅里叶级数表示。但对数值计算来说，这是无法实现的。在实际应用中，我们可以用有限项的傅里叶级数求和来逼近。即对有限项和： （5.3） 当N值取得较大时，上式就是原周期信号的一个很好的近似。上式常称作截断傅里叶级数表示。 实验内容 将例5-1在计算机上练习一遍，观察运行结果是否正确； 解; function [A_sym,B_sym]=CTFShchsym syms t n k x T=5;tao=0.2*T;a=0.5; if nargin4;Nf=6;end if nargin5;Nn=32;end x=time_fun_x(t); A0=2*int(x,t,-a,T-a)/T; As=int(2*x*cos(2*pi*n*t/T)/T,t,-a,T-a); Bs=int(2*x*sin(2*pi*n*t/T)/T,t,-a,T-a); A_sym(1)=double(vpa(A0,Nn)); for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn)); B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); end if nargout==0 c=