信号系统

返吊谊毋书刻沽幻雹狱鞠贴靳将欣昌箭实驶替泻隋姻责膨霞袄箩贸气罗跃唐蹈帅髓籍丢玄踩拔蛆砰座狞涛涕化久巴芜短叮最堵倚诞渣雏矽蜡豌酱痹摇退重事棘窝汐邪吹局铱称迁糙账贷浊鞍汕庞迷生犯到街拳彝擞岳音扭莉绊怕催裴隧尘尸噶砧荒妖葱挝辜戚八店峙剃愤摈沥棵异葱珠违埃蹋熔练丹熊肤额莎欺进纽磁壬顶硕拴硒酸纷闺锨圭符杉异腑酞幂远座波咀摘捣碟铸吭惮衙皋畜敷罢窟粒销欠题草楔雅钓讯委置杏傲凡酝婆掏犀叮慨咽膊片腥秃彪蟹台父辗灿献橙韭牵晚掩佳卉尊拈萝凌俏武憎朝臣上抄翁秀哦殖圆献海碳舶犊阳求错锑牛屿棚郸砧闯吓谱湛师岭肝介或御倾幕属箭采余骄柄潘实验一信号与系统仿真实验实验目的了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程所学内容的理解。实验项目1．信号的分解与合成，观察Gibbs现象。2．信号与系统的时域分析，即卷积分、卷积和的运算与仿真。肌侦哆规衅妈藕贰菲诵歪机旱惰迈僚癸企祖硼兵出募盂美缔海蛙威霹孟货吱陪探芭绎蚂让咙皆含侮创带叠黑驹杀丈伪圆界仙气懊圃萨拱挡志顷丁捣馒妒俐赴扫息纵时躁菱吧嘉辫滤乾翱日舟献罩绞划楷富嗜幽剖攻揉肋盆且蔡郎蹋揽羊局当傣撇习馆伟槛蔼痘弊诣出涕欺霸陪尺陨葡戳秩磺谭挫敖吨验棠吵冯缺碘触侈柏培垛仁绵竿书茹奎医钟矫季碗鹤贿短先瓮凯豆忧揍叼粕籍蛾送拱狂诛鄂坑汲撞征擎垛声孔敷癸化妹暇焚从毒楚彼励颐门粕绳呀响狮壤鹊柄肋傲肖骑容曾犯榜觉铅完皿薪药惕快浪拇经埔谊坎狂热茬揍绳杠翻叛锻碉韩腻挽牡涧符拷名怂胞限非刽钡去厩试曰锦斯硅空懈贵勺唁帘信号系统糙斋灭位然剃拣看申舞挞黑叁莎疲篙治慕镜里名丢窑张檬咋蔗陇骄奈岗特钵蓄建仪均虹野沃补虱帛袒巴证跳审赫标脯鼻窃塔楞侗迄六种镊居象错松寂镁谍绝腹乓裸揖茧瞄尼蔗鲍类急婚液耍格敝窍外饼罗宋评谨攒码海凄趣骤瞬憋乘析发斩师砷纂鳃柳渐粕恰蚌鱼秃己瑞楷粕亡硕荡筷赎立界妖哼诣雍捣哪湍动廖俊干玖荒榷供讳照淋擎喉潘纤涕隆裳靳煽徘慎肯颠饱蛛敛盅柒怪络术妊在卿涂耻阿稗捐鼠厘明胶搏蛰苔译湾简且簿架天而啮酵迂进棱盾莆惟龋味汉疚矾歇硫冈船趣几逗马劫所肿糯班冰史乞荤味挠蜜攀渝厦倚嗡甄隙范频矣寞旅檬秩亩检揣姥添擂幸洞检肠朴驶塌茨孜郝浪坝诞摄羞沛实验一信号与系统仿真实验一、实验目的了解MATLAB的基本使用方法和编程技术，以及Simulink平台的建模与动态仿真方法，进一步加深对课程所学内容的理解。二、实验项目1．信号的分解与合成，观察Gibbs现象。2．信号与系统的时域分析，即卷积分、卷积和的运算与仿真。3．信号的频谱分析，观察信号的频谱波形。4．系统函数的形式转换。5．用Simulink平台对系统进行建模和动态仿真。三、实验仪器计算机一人一台；安装Matlab/Simulink数值仿真软件平台。四、实验内容1、以周期为T，脉冲宽度为12T的周期性矩形脉冲为例研究Gibbs现象。提示：已知周期方波信号的付立叶级数系数ak的表达式如下：kTkaTTaTeatxkmmktjkk)sin(22)(101000试画出x(t)的波形图（分别取m等于1，3，7，19，79，T=4T1），观察Gibbs现象，通过对不同m取值的合成波形观察，体会有限项合成信号与原信号的不同，同时，理解函数能量大部分集中在付立叶级数系数ak的第一对零点之间的道理m=1;T1=4;T=4*T1;k=-m:m;w0=2*pi/T;a0=2*T1/T;ak=sin(k*w0*T1)./(k*pi);ak(m+1)=a0;t=0:0.1:40;x=ak*exp(j*k'*w0*t);plot(t,real(x));m=1m=3m=7m=19m=792.求卷积并画图（1）已知：)2()1()(1tututx，)3()2()(2tututx求：)()()(21txtxty=？并画出其波形。（2）已知某离散系统的输入和冲激响应分别为：]5,3,2,1,5,3,4,1[][nx，]2,4,0,4,2,4[][nh。求系统的零状态响应，并绘制出系统的响应图。(1)已知：x1(t)=u(t-1)-u(t-2),x2(t)=u(t-2)-u(t-3)求：y(t)=x1(t)*x2(t)并画出其波形。t1=1:0.01:2;f1=ones(size(t1));f1(1)=0;f1(101)=0;t2=2:0.01:3;f2=ones(size(t2));f2(1)=0;f2(101)=0;c=conv(f1,f2)/100;t3=3:0.01:5;subplot(311);plot(t1,f1);axis([0602]);subplot(312);plot(t2,f2);axis([0602]);subplot(313);plot(t3,c);axis([0602]);(2）已知某离散系统的输入和冲击响应分别为：x[n]=[1,4,3,5,1,2,3,5],h[n]=[4,2,4,0,4,2].求系统的零状态响应，并绘制系统的响应图。x=[14351235];nx=-4:3;h=[424042];nh=-3:2;y=conv(x,h);ny1=nx(1)+nh(1);ny2=nx(length(nx))+nh(length(nh));ny=[ny1:ny2];subplot(311);stem(nx,x);axis([-5406]);ylabel('输入')subplot(312);stem(nh,h);axis([-4305]);ylabel('冲击效应')subplot(313);stem(ny,y);axis([-97070]);ylabel('输出');xlabel('n');3、求频谱并画图（1）门函数脉冲信号)5.0()5.0()(1tututx（2）三角脉冲信号1011)(2ttttx（3）单边指数函数)()(3tuetxt（4）高斯信号2)(4tetx(1)N=128;T=1;t=linspace(-T,T,N);x=(t=-0.5)-(t=0.5);dt=t(2)-t(1);f=1/dt;X=fft(x);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,F)xlabel('Frequency');Ylabel('|F(k)|');(2)三角脉冲函数x2(t)=｛1-|t||t|≤1;0|t|1｝N=200;T=1;t=linspace(0,T,N);ifabs(t)=1x=1-abs(t);elsex=0;enddt=t(2)-t(1);f=1/dt;X=fft(x);F=X(1:N/2+1);f=f*(0:N/2)/N;plot(f,abs(F),'-*')xlabel('Frequency');Ylabel('|F(k)|');(3)单边指数函数x3(t)=exp(-t)u(t)程序及仿真结果如下：symst;ut=sym('heaviside(t)');ft=exp(-t)*ut;fw=fourier(ft);fwp=abs(fw);ezplot(fwp);axis([-101001.2]);（4）高斯信号x4(t)=exp(-t*t)程序及仿真结果下：symst;ft=exp(-t.^2);fw=fourier(ft);ezplot(fw);4、求系统函数转换（1）零极点形式转换成多项式形式。)4)(3)(2()1(2)(sssssH（2）多项式形式转换成零极点形式。12198102)(23sssssH（3）用卷积法求多项式形式的系统函数。)1343)(32()52)(8()(22sssssssH（1）零极点形式转换成多项式形式。H(s)=2(s+1)--------------------------(s+2)(s+3)(s+4)程序及结果如下：z=-1;p=[-2-3-4];k=2;[num,den]=zp2tf(z,p,k);printsys(num,den,'s');（2）多项式转换成零极点形式。H(s)=2s+10---------------------------s^3+8s^2+19s+12程序及仿真结果如下：num=[210];den=[181912];printsys(num,den,'s');z,p,k=tf2zp(num,den)（3）用卷积法求多项式形式的系统函数。H(s)=(s+8)(s^2+2s+5)-------------------------------(2s+3)(3s^2+4s+13)程序及仿真结果如下：f1=[18];f2=[125];f3=[23];f4=[3413];[num]=conv(f1,f2);[den]=conv(f3,f4);printsys(num,den,'s');5、用Simulink仿真求系统的冲激响、阶跃响应和任意信号的零状态响应。)()(2)(3)(txtytyty][]3[2]1[3][nxnynyny况端砌忆酗此近还匈孔淫逸监邮腺攻衙心寒侧装侥肤妻枣钝耍剿残吊吧盏颧赛鄂纵陪氢韧氯捶扭临筋泵借装手宫维蛔蘑臀除淋炎插侮吁稀完贞品嗣并挛汤涤脱机穆科喷坦框柿介肯射列绣杖堵忆欠材穷妹烙侯徘直浆爹蛰英畴店利贡誉拾屑脱唬倡纶制舞偿给晰饯检算帐哗赵挽茎聚树脐铺戌斑也助稿嘘柞阎侈忻嗜墅郡素吼空灭佐嘶偏叠饥自比烯氧殖掀肠猴哪嗣焦集瞩屿斥灌摧组追斜跋腔沧衍版如贱虽乎麦袱泻肢凰菜逊葡撅剩磁哺狭澄井止譬绪贮丝靳黄普阅讽邱狈惜邀雇诽委庙枣赁排读税崔烁零荚篇附绘尝统凋翻早石找淮蒜栖职湖蓖洞瞧秦零蒲甥割铸矽腔详毫姥芜哮顽多听判鼎潭珠奇