matlab+nnf.m,中南大学-信号与系统matlab实验报告.doc

中南大学-信号与系统matlab实验报告.doc

实验一 基本信号的生成 1 实验目的 学会 使用 MATLAB 产生各种常见的连续时间信号与离散时间信号； 通过 MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的理解； 熟 悉 MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠定基础。 2 实验内容 运行以上九个例子程序，掌握一些常用基本信号的特点及其 MATLAB 实现方法；改变有关参数，进一步观察信号波形的变化。 在 k 1010 范围内产生并画出以下信号 a 1f k k ； b 2f k k2 ； c 3f k k-4 ； d 4f k 2 k 2 k -4 。 源程序 k-1010; f1kzeros1,10,1,zeros1,10; subplot2,2,1 stemk,f1k titlef1k f2kzeros1,8,1,zeros1,12; subplot2,2,2 stemk,f2k titlef2k f3kzeros1,14,1,zeros1,6; subplot2,2,3 stemk,f3k titlef3k f4k2*f2k-f3k; subplot2,2,4 stemk,f4k titlef4k 在 k 031 范围内产生并画出以下信号 a kk1 44f k sin co s ； b 2 k2 4f k cos ； c kk3 48f k sin c o s 。 请问这三个信号的基波周期分别是多少 源程序 k031; f1ksinpi/4*k.*cospi/4*k; subplot3,1,1 stemk,f1k titlef1k f2kcospi/4*k.2; subplot3,1,2 stemk,f2k titlef2k f3ksinpi/4*k.*cospi/8*k; subplot3,1,3 stemk,f3k titlef3k 其中 f1k的基波周期是 4, f2k的基波周期是 4, f3k的基波周期是 16。 实验 二 信号的基本运算 1 实验目的 学会 使用 MATLAB 完成信号的一些基本运算； 了解复杂信号由基本信号通过尺度变换、翻转、平移、相加、相乘、差分、求和、微分及积分等运算来表达的方法； 进一步 熟 悉 MATLAB 的基本操作与编程，掌握其在信号分析中的运用特点与使用方式 。 2 实验内容 运行以上三个例题程序，掌握信号基本运算的 MATLAB实现方法；改变有关参数，考察相应信号运算结果的变化特点与规律。 已知信号 ft如下图所示 -4 -3 -2 -1 0 1 2 3 400 . 20 . 40 . 60 . 811 . 21 . 41 . 61 . 82tft给定信号 f t a 用 MATLAB 编程复现上图； 作业题 2 a t-60.0016; ft1tripulst,6,0.5; subplot2,1,1 plott,ft1 titleft b 画出 2 2ft 的波形； b t-60.0016; ft1tripuls2*1-t,6,0.5; subplot1,1,1 plott,ft1 titlef2*1-t c 画出 dftdt 的波形； c h0.001;t-6h6; yttripulst,6,0.5; y1diffyt*1/h; plott1lengtht-1,y1 titledft/dt d 画出 t f d的波形。 d t-60.16; for x1lengtht y2xquadtripulst,6,0.5,-3,tx; end plott,y2 titleintegral of ft 实验 三 系统的时域分析 1 实验目的 学习并掌握连续时间系统的 零状态 响应、 冲激响应和阶跃响应 的 MATLAB求解方法 ； 学习并掌握离散时间系统的 零状态 响应、 冲激响应和阶跃响应 的 MATLAB求解方法 ； 进一步 深刻理解连续时间系统和离散时间系统的系统函数零极点对系统特性的影响 ； 学习并掌握卷积的 MATLAB计算方法。 2 实验内容 运行以上五个例题程序，掌握求解 系统响应的 MATLAB分析方法；改变模型参数，考察 系统响应 的变化特点与规律。 设离散系统可由下列差分方程表示 y k - y k - 1 0 .9 y k - 2 f k 计算 20 100k 时的系统冲激响应。 源程序 k-20100; a1 -1 0.9; b1; himpzb,a,k; stemk,h; xlabelTimesec ylabelyt 设 0.9 kh k u k ，输入 1 0 f k u k u k ，求系统输出 y k f k h k。 (取 1050k ) 源程序 k-1050; ukzeros1,10,ones1,51; u1kzeros1,20,ones1,41; hk0.9.k.*uk; fkuk-u1k; ykconvhk,fk; stem0lengthyk-1,yk; 已知滤波器的传递函数 10.2 2 1 0.8Hz z 输入信号为 2 sin 0 .0 5 , f t t t t 为随机信号。试绘出滤波器的输出信号波形。(取 0100t ) 源程序 R101; drand1,R-0.5; t0100; s2*sin0.05*pi*t; fsd; subplot2,1,1; plott,d,g-.,t,s,b--,t,f,r-; xlabelTime index t; legenddt,st,ft; title处理前的波形 b0.22 0;a1 -0.8; yfilterb,a,f; subplot2,1,2; plott,s,b--,t,y,r-; xlabelTime index t; legendst,yt; title滤波器输出波形 实验 四 周期信号的频域分析 1 实验目的 掌握周期信号傅立叶级数分解与合成的计算公式 掌握利用 MATLAB 实现周期信号傅立叶级数分解与综合方法 理解并掌握周期信号频 谱特点 2.实验内容 1、 仿照例程，实现下述周期信号的傅立叶级数分解与合成 要求 ( a)首先，推导出求解0a，na，nb的公式，计算出前 10 次系数； ( b)利用 MATLAB 求解0a，na，nb的值，其中na，nb求解前 10 次系数，并给出利用这些系数合成的信号波形。 t ft 1 -3 -4 5 4 1 O ( a)设周期信号 tf 的周期为 1T ，角频率11122 Tf ，且满足狄里赫利条件， 则该周期信号可以展开成傅立叶级数。 ( 1)三角形式傅立叶级数 dttntfTbdttntfTadttfTatnbtnaatbtatbtatbtaatfTTnTTnTTnnnnnnnn22112211221011110222211110111111s i n2c os21s i n c os .s i nc os.s i nc oss i nc os( 2)指数形式傅立 叶级数 n tjnn neFtf ,3,2,1,0, 1 dtetfTF TT tjnn 2121 111 ( b)求解0a，na，nb及合成信号波形所用程序 function A_sym,B_symCTFShchsym 采用符号计算求一个周期内连续时间函数 f 的三角级数展开系数 ,再用这些 展开系数合成连续时间函数 f.傅立叶级数 函数的输入输出都是数值量 Nf6 谐波的阶数 Nn 输出数据的准确位数 A_sym 第 1 元素是直流项，其后元素依次是 1,2,3.次谐波 cos 项展开系数 B_sym 第 2,3,4,.元素依次是 1,2,3.次谐波 sin 项展开系数 tao1 tao/T0.2 syms t n k x T4; taoT/4; a-1.5; if nargin4 Nf10; end if nargin5 Nn32; end xtime_fun_xt; A0intx,t,a,Ta/T; 求出三角函数展开系数 A0 As2/T*intx*cos2*pi*n*t/T,t,a,Ta; 求出三角函数展开系数 As Bs2/T*intx*sin2*pi*n*t/T,t,a,Ta; 求出三角函数展开系数 Bs A_sym1doublevpaA0,Nn; 获取串数组 A0 所对应的 ASC2 码数 值数组 for k1Nf A_symk1doublevpasubsAs,n,k,Nn; 获取串数组 A 所对应的 ASC2 码数值数组 B_symk1doublevpasubsBs,n,k,Nn; 获取串数组 B 所对应的 ASC2 码数值数组 end ; if nargout0 cA_sym; dispc; 输出 c 为三角级数展开系数 第 1 元素是直流项，其后元素依次是 1,2,3.次谐波cos 项展开系数 dB_sym; dispd; 输出 d 为三角级数展开系数 第 2,3,4,.元素依次是 1,2,3.次谐波 sin 项展开系数 t-3*T0.013*T; f0c1; 直流 f1c2.*cos2*pi*1*t/Td2.*sin2*pi*1*t/T; 基波 f2c3.*cos2*pi*2*t/Td3.*sin2*pi*2*t/T; 2 次谐波 f3c4.*cos2*pi*3*t/Td4.*sin2*pi*3*t/T; 3 次谐波 f4c5.*cos2*pi*4*t/Td5.*sin2*pi*4*t/T; 4 次谐波 f5c6.*cos2*pi*5*t/Td6.*sin2*pi*5*t/T; 5 次谐波 f6c7.*cos2*pi*6*t/Td7.*sin2*pi*6*t/T; 6 次谐波 f7c8.*cos2*pi*7*t/Td8.*sin2*pi*7*t/T; 7 次谐波 f8c9.*cos2*pi*8*t/Td9.*sin2*pi*8*t/T; 8 次谐波 f9c10.*cos2*pi*9*t/Td10.*sin2*pi*9*t/T; 9 次谐波 f10c11.*cos2*pi*10*t/Td11.*sin2*pi*10*t/T; 10 次谐波 f11f0f1f2; 直流 基波 2 次谐波 f12f11f3; 直流 基波 2 次谐波 3 次谐波 f13f12f4f5f6; 直流 基波 2 次谐波 3 次谐波 4 次谐波 5 次谐波 6 次谐波 f14f13f7f8f9f10; 010 次 subplot2,2,1 plott,f0f1,hold on ytime_fun_et; 调用连续时间函数 -周期矩形脉冲 plott,y,r title直流 基波 axis-8,8,-0.5,1.5 subplot2,2,2 plott,f12,hold on ytime_fun_et; plott,y,r title1-3 次 谐波 直流 axis-8,8,-0.5,1.5 subplot2,2,3 plott,f13,hold on ytime_fun_et; plott,y,r title1-6 次谐波 直流 axis-8,8,-0.5,1.5 subplot2,2,4 plott,f14,hold on ytime_fun_et; plott,y,r title1-10 次谐波 直流 axis-8,8,-0.5,1.5 hold off end function ytime_fun_et 该函数是 CTFShchsym.m 的子函它由符号函数和表达式写成 a1.5; T4; h1; taoT/4; t-3*T0.013*T; e11/21/2.*signt-0.5tao/2; e21/21/2.*signt-0.5-tao/2; yh.*e1-e2; 连续时间函数 -周期矩形脉冲 function xtime_fun_xt 该函数是 CTFShchsym.m 的子函数。它由符号变量和表达式写成。 h1; x1symHeavisidet*h; xx1-symHeavisidet-1*h; 源程序修改 function A_sym,B_symCTFShchsym 采用符号计算求一个周期内连续时间函数 f 的三角级数展开系数 ,再用这些 展开系数合成连续时间函数 f.傅立叶级数 函数的输入输出都是数值量 Nf6 谐波的阶数 Nn 输出数据的准确位数 A_sym 第 1 元素是直流项，其后元素依次是 1,2,3.次谐波 cos 项展开系数 B_sym 第 2,3,4,.元素依次是 1,2,3.次谐波 sin 项展开系数 tao1 tao/T0.2 syms t n k x T5; tao0.2*T; a0.5; if nargin4 Nf6; end if nargin5 Nn32; end xtime_fun_xt; A0intx,t,-a,T-a/T; 求出三角函数展开系数A0 As2/T*intx*cos2*pi*n*t/T,t,-a,T-a; 求出三角函数展开系数 As Bs2/T*intx*sin2*pi*n*t/T,t,-a,T-a; 求出三角函数展开系数 Bs A_sym1doublevpaA0,Nn; 获取串数组 A0 所对应的 ASC2 码数值数组 for k1Nf A_symk1doublevpasubsAs,n,k,Nn; 获取串数组 A 所对应的 ASC2 码数值数组 B_symk1doublevpasubsBs,n,k,Nn; 获取串数组 B 所对应的 ASC2 码数值数组 end if nargout0 cA_sym; dispc 输出 c 为三角级数展开系数 第 1 元素是直流项，其后元素依次是 1,2,3.次谐波 cos 项展开系数 dB_sym; dispd 输出 d 为三角级数展开系数 第 2,3,4,.元素依次是 1,2,3.次谐波 sin 项展开系数 t-8*a0.01T-a; f0c1; 直流 f1c2.*cos2*pi*1*t/5d2.*sin2*pi*1*t/5; 基波 f2c3.*cos2*pi*2*t/5d3.*sin2*pi*2*t/5; 2 次谐波 f3c4.*cos2*pi*3*t/5d4.*sin2*pi*3*t/5; 3 次谐波 f4c5.*cos2*pi*4*t/5d5.*sin2*pi*4*t/5; 4 次谐波 f5c6.*cos2*pi*5*t/5d6.*sin2*pi*5*t/5; 5 次谐波 f6c7.*cos2*pi*6*t/5d7.*sin2*pi*6*t/5; 6 次谐波 f7f0f1f2; 直流 基波 2 次谐波 f8f7f3; 直流 基波 2 次谐波 3 次谐波 f9f8f4f6; 直流 基波 2 次谐波 3 次谐波 4 次谐波 6 次谐波 subplot2,2,1 plott,f0f1,hold on ytime_fun_et; 调用连续时间函数 -周期矩形脉冲 plott,y,r title周期矩形波的形成 直流 基波 axis-4,4.5,-0.5,1.5 subplot2,2,2 plott,f7,hold on ytime_fun_et; plott,y,r title周期矩形波的形成 直流 基波 2 次谐波 axis-4,4.5,-0.5,1.5 subplot2,2,3 plott,f8,hold on ytime_fun_et; plott,y,r title直流 基波 2 次谐波 3 次谐波 axis-4,4.5,-0.5,1.5 subplot2,2,4 plott,f9,hold on ytime_fun_et; plott,y,r title基波 2 次谐波 3 次谐波 4 次谐波 6 次谐波 axis-4,4.5,-0.5,1.5 end function ytime_fun_et 该函数是 CTFShchsym.m 的子函它由符 号函数和表达式写成 a0.5; T5; h1; tao0.2*T; t-8*a0.01T-a; e11/21/2.*signttao/2; e21/21/2.*signt-tao/2; yh.*e1-e2; 连续时间函数 -周期矩形脉冲 function xtime_fun_xt 该函数是 CTFShchsym.m 的子函数。它由符号变量和表达式写成。 h1; x1symHeavisidet0.5*h; xx1-symHeavisidet-0.5*h; 2、 已知周期为 T4 的三角波，在第一周期( -2t2)内表示成 ttx 1 ，试用 MATLAB求该信号的傅立叶级数，并绘制它的频谱图。将它的频谱与方波的频谱图做比较。 function A_sym,B_symCTFSshbpsymT,Nf 采用符号计算求 0,T内时间函数的三角级数展开系数。 函数的输入输出都是数值量 Nn 输出数据的准确位数 A_sym 第 1 元素是直流项，其后元素依次是 1,2,3.次 谐波 cos 项展开系数 B_sym 第 2,3,4,.元素依次是 1,2,3.次谐波 sin 项展开系数 T Tm*tao， 信号周期 Nf 谐波的阶数 m mT/tao周期与脉冲宽度之比 ,如 m4,8,16,100 等 tao 脉宽 taoT/m syms t n y if nargin3 Nfplease 所需展开的最高谐波次数 Nf; end Tplease 信号的周期 T; if nargin5 Nn32; end ytime_fun_st; A02/T*inty,t,0,T/2; As2/T*inty*cos2*pi*n*t/T,t,0,T/2; Bs2/T*inty*sin2*pi*n*t/T,t,0,T/2; A_sym1doublevpaA0,Nn; for k1Nf A_symk1doublevpasubsAs,n,k,Nn; B_symk1doublevpasubsBs,n,k,Nn; end if nargout0 AnfliplrA_sym; 对 A_sym 阵左右对称交换 An1,k1A_sym1; A_sym 的 1*k 阵扩展为 1*k1阵 AnfliplrAn; 对扩展后的 S1 阵左右对称交换回原位置 BnfliplrB_sym; 对 B_sym 阵左右对称交换 Bn1,k10; B_sym 的 1*k 阵扩展为 1*k1阵 BnfliplrBn; 对扩展后的 S3 阵左右对称交换回原位置 FnRAn/2-i*Bn/2; 用三角函数展开系数 A、 B 值合成付里叶指数系数 FnLfliplrFnR; NNf*2*pi/T; k2-N2*pi/TN; FnFnL,FnR2end; subplot3,3,3 xtime_fun_et; 调用连续时间函数 -周期矩形脉冲 subplot2,1,1 stemk2,absFn; 画出周期矩形脉冲的频谱( TM*tao) title连续时间函数周期三角波脉冲的双边幅度谱 axis-80,80,0,0.12 line-80,80,0,0,color,r line0,0,0,0.12,color,r end function xtime_fun_et 该函数是 CTFSshbpsym.m 的子函数。它由符号变量和表达式写成。 t 是时间数组 T 是周期 dutytao/T0.2 T5; t-2*T0.012*T; taoT/5; xrectpulst,tao; 产生一个宽度 tao1 的矩形脉冲 subplot2,2,2 plott,x hold on xrectpulst-5,tao; 产生一个宽度 tao1 的矩形脉 ,中心位置在 t5 处 plott,x hold on xrectpulst5,tao; 产生一个宽度 tao1 的矩形脉 ,中心位置在 t-5 处 plott,x title周期为 T5，脉宽 tao1 的矩形脉冲 axis-10,10,0,1.2 function ytime_fun_st syms t y1-abst; x1symHeavisidet2; xx1-symHeavisidet-2; yy*x; ezplott,y,-10,10 grid 源程序修改 function A_sym,B_symCTFSshbpsymT,Nf 采用符号计算求 0,T内时间函数的三角级数展开系数。 函数的输入输出都是数值量 Nn 输出数据的准确位数 A_sym 第 1元素是直流项，其后元素依次是 1,2,3.次谐波 cos项展开系数 B_sym 第 2,3,4,.元素依次是 1,2,3.次谐波 sin项展开系数 T Tm*tao， 信号周期 Nf 谐波的阶数 m mT/tao周期与脉冲宽度之比 ,如 m4,8,16,100等 tao 脉宽 taoT/m syms t n y if nargin3 Nfplease 所需展开的最高谐波次数 Nf; end Tplease 信号的周期 T; if nargin5 Nn32; end ytime_fun_st; A02/T*inty,t,0,T; As2/T*inty*cos2*pi*n*t/T,t,0,T; Bs2/T*inty*sin2*pi*n*t/T,t,0,T; A_sym1doublevpaA0,Nn; for k1Nf A_symk1doublevpasubsAs,n,k,Nn; B_symk1doublevpasubsBs,n,k,Nn; end if nargout0 AnfliplrA_sym; 对 A_sym阵左右对称交换 An1,k1A_sym1; A_sym的 1*k阵扩展为 1*k1阵 AnfliplrAn; 对扩展后的 S1阵左右对称交换回原位置 BnfliplrB_sym; 对 B_sym阵左右对称交换 Bn1,k10; B_sym的 1*k阵扩展为 1*k1阵 BnfliplrBn; 对扩展后的 S3阵左右对称交换回原位置 FnRAn/2-i*Bn/2; 用三角函数展开系数 A、 B值合 成付里叶指数系数 FnLfliplrFnR; NNf*2*pi/T; k2-N2*pi/TN; FnFnL,FnR2end; subplot3,3,3 xtime_fun_et; 调用连续时间函数 -周期矩形脉冲 subplot2,1,1 stemk2,absFn; 画出周期矩形脉冲的频谱( TM*tao) title连续时间函数周期矩形脉冲的双边幅 度谱 axis-80,80,0,0.12 line-80,80,0,0 line0,0,0,0.12 end function xtime_fun_et 该函数是 CTFSshbpsym.m的子函数。它由符号变量和表达式写成。 t 是时间数组 T 是周期 dutytao/T0.2 T5; t-2*T0.012*T; taoT/5; xrectpulst,tao; 产生一个宽度 tao1的矩形 脉冲 subplot2,2,2 plott,x hold on xrectpulst-5,tao; 产生一个宽度 tao1的矩形脉 ,中心位置在 t5处 plott,x hold on xrectpulst5,tao; 产生一个宽度 tao1的矩形脉 ,中心位置在 t-5处 plott,x title周期为 T5，脉宽 tao1的矩形脉冲 axis-10,10,0,1.2 function ytime_fun_st 该函数是 CTFSshbpsym.m的子函数。它由符号变量和表达式写成。 syms a a1 Tplease 信号的周期 T; M周期与脉冲宽度之比 M; A1; taoT/M; atao/2; y1symHeavisideta1*A; yy1-symHeavisidet-a1*A; ysubsy,a1,a; ysimpley; 实验 五 非 周期信号的频域分析 1 实验目的 理解非周期信号的频域分析方法； 掌握典型信号的幅度谱和相位 谱； 理解信号的调制特性； 掌握傅里叶变换的性质尺度变换、时移、频移、卷积定理、对称性、微分特性。 2 实验 内容 1.试计算宽度为 2、幅度为 1 的三角波信号在 0fmHz 范围内信号的能量。取fm0.110Hz。 syms t; flinspace0.1,10,256; Nlengthf; wzeros1,N; for k1N wkquadlt2*sinct.*sinct.*sinct.*sinct,0,fk; end y0; for k1N yywk; end y plotf,w; xlabelHz; ylabelE; y 168.2699 2.试计算宽度和幅度均为 1 的方波信号在 0fmHz 频谱范围内所包含的信号能量。 syms t; flinspace0,5,256; Nlengthf; wzeros1,N; for k1N wkquadlt2*sinct.*sinct,0,fk; end plotf,w; xlabelHz; ylabelE;