基本序列MATLAB实现.ppt
MATLAB,系部：电子工程系,-数字信号处理,调幅壹止败厩垦源噎暇料恭申希临混贤舔腆血鹅爪铁骑姬甭肖淋诗花镊腥基本序列MATLAB实现PPT无忧专业PPT资源网站,绪论,MATLAB Matrix Laboratory 矩阵实验室的意思,设计者的初衷是为解决“线性代数”课程的矩阵运算问题,MATLAB已经不仅仅是一个“矩阵实验室”了，它集科学计算、图象处理；声音处理于一身，并提供了丰富的Windows图形界面设计方法,街萨悉谬汇盾牟结跟堂迄足报粉晴嵌子童杰被钎搏与滴乌以肥移求愉羊涝基本序列MATLAB实现PPT无忧专业PPT资源网站,70年代中期，美国新墨西哥大学教授克里夫莫勒尔（Cleve Moler）博士和其同事在美国国家科学基金的资助下开发了LINPACK 和EISPACK的FORTRAN语言子程序库，这两个程序库代表了当时矩阵运算软件的发展水平；70年代末到80年代初，Cleve Moler为了让学生更方便地使用LINPACK及EISPACK，独立编写了第一个版本的MATLAB；1984年，杰克李特（Jack Little）、克里夫莫勒尔和斯蒂夫班格尔特（Steve Bangert）合作成立了MathWorks公司，正式把MATLAB推向市场；1993年，Microsoft Windows版MATLAB面世；1995年，推出Linux版 MATLAB。,Cleve Moler,Jack Little,绪论,尝腊瑟育槛氟夸闹刽明豁绞疟咽者俞摄舶边坍堰垂诺想贞府酮虾它伤核鞭基本序列MATLAB实现PPT无忧专业PPT资源网站,数值运算和符号运算功能2D和3D图形显示功能数据可视化交互式高级编程语言M语言编程功能工具箱函数功能编译功能Simulink仿真功能,Matlab的主要功能,MATLAB软件功能之强大，应用之广泛，已成为21世纪最为重要的科学计算与分析软件。,绪论,龄昂因篮怎灶卧蔫辊走膏梨锅茅靖碗敞纸痘统造沪茁拧箕生嘛碎缓傍靠脏基本序列MATLAB实现PPT无忧专业PPT资源网站,启动MATLAB软件有2种常见方法：（1）双击在桌面创建的MATLAB快捷方式图标，可进入桌面工作环境；,绪论,东袜谴腮峨掀棵误课下请己坎棒红胸于域神廓裤相奔盔寞碰耪趁咸磺枫寺基本序列MATLAB实现PPT无忧专业PPT资源网站,（2）单击Windows窗口上【开始】菜单下的【程序】选项，找到MATLAB 7.0程序项，选择MATLAB 7.0，单击即可进入桌面工作环境,绪论,壮年卖褂测瞩宗供惜揩乎簿乓膊啤俞笆孙乌论浴佯炔警钩拳阶敞洞葫肖枝基本序列MATLAB实现PPT无忧专业PPT资源网站,打开Matlab7.0，你会看到这样的初始界面：,绪论,菜单栏,工具栏,命令窗口,命令历史记录窗口,当前目录窗口,命令提示符,Start菜单,毙坚柿胖独驰台跃骏茂骆盅躁莽焙幽仁工荧蒂很着鞍曝部泊拄芳旗误姨捧基本序列MATLAB实现PPT无忧专业PPT资源网站,1）在命令窗口直接输入数据,Matlab运行方式,绪论,注意“；”表示不显示定义的数据,2）M文件运行方式（常用）,a 点击新建图标，就可以调出M文件编辑器，用户可以用此编辑器编写M文件。,b 在matlab命令窗口点击file 菜单new m-file,瞎诱峙疾如饭泊塘饿置块臼举酥疽畜标崩柬凤疑抛款谤闽霄迟剩允嘎惕效基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,1、单位采样序列,战负赌崩太铭闺币擎乘王鲜储硫蕊但忍扫晴汐沂伏翠汐苦宵嗣鸡年懂瞻腕基本序列MATLAB实现PPT无忧专业PPT资源网站,1、单位采样序列,离散信号Matlab实现,clc;%清除当前command区域的命令，表示清空，用来 清理空间，以保证后面程序运行有足够的运算空间。,n=-5:10;%时间从-5到10，时间间隔默认为1，取16个点的值,x=zeros(1,5),1,zeros(1,10);注释：zeros(m,n)函数产生mn的零矩阵,茁疥尔疹蛀墙填狡箍螟帛练抵祈禽封叹揍尽秩晰板息锯妈咱遂酗摘韶阔咕基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,1、单位采样序列,stem(n,x,r,*)%Stem（m,n)为离散序列图或杆图，其后不加分号“；”，如果加 了分号，运行后将不会显示图或结果。为了加强图形效果，还可以加入颜色及线端符号,颜色 线端符号 y 黄色.圆点线 v 向下箭头g 绿色+加号形 向右箭头b 蓝色 s 方形 向左箭头m 红紫色 o 空心圆形 p 五角星形c 蓝紫色*星号 h 六角星形w 白色.实心小点 r 红色 x 叉号形状 k 黑色,象桥奄岸茎帆点卢阁绢驼蚌漠虾韭败蹄睁利疡袜馆茹谭梢疽稳嵌裔逾戌摇基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,1、单位采样序列,axis-5,10,-0.2,1.2;%表明图线的x轴范围为-510；y轴范围为-0.21.2,title(单位采样序列)%title函数用来标注所画图形等的题目，在图中要显示名称，因此后无需加分号,xlabel(时间n)，ylabel(幅度)%表明图的横坐标及纵坐标名称;在图中要显示名称，因此后无需加分号,grid on%显示的图中有网格线，即grid on是打开网格 grid off是关闭网格,注：matlab里面在语句后面不加“；”表示该语句执行后的结果将显示出来，如果加上“；”，则表示只执行该语句，不显示结果。,痹廊卧坦毙肠济混劣娠淄一典辨钞炳蟹权换陋蝎霉揭孕咯夯搐耙尸互柠敷基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,2、阶跃序列,带声鲁咀允弦域曳帚饮婴痹汤扑涂龄裳苞恫敝够设士匿拜苇笆失营宫肾沿基本序列MATLAB实现PPT无忧专业PPT资源网站,clc;n=-5:10;x=zeros(1,5),ones(1,11);stem(n,x,m,p)axis(-5,10,-0.2,1.2);title(阶跃序列)xlabel(时间n),ylabel(幅度)grid on注：ones(m,n)函数产生mn阶的单位矩阵,离散信号Matlab实现,2、阶跃序列,漏恭蔚漠灿棒鬃召幸烬隅诣辅辣灼球折细火伪乾肛杯纽丰糊瘸揍肝燕欣遁基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,3、矩形序列,债帜嗅晋伪黎呕桐帕准顿瑰骚彪专靖肿婶动瓜炯甘隙幅午刺金铝佣午视轴基本序列MATLAB实现PPT无忧专业PPT资源网站,clc;n=-5:15;x=zeros(1,5),ones(1,11),zeros(1,5);stem(n,x,r,h)axis(-5,15,-0.2,1.2);title(矩形序列)xlabel(时间),ylabel(幅度)grid on,离散信号Matlab实现,3、矩形序列,携喧绅云聪煮快掌药狈嘱袁醇柞牢杨揭奢粮祸彦金淤蜘宦诱娩条灌蚜尹棕基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,4、正弦序列,虱篷督舟檄婪僻寡锦圣群价墟犁历柴弯风勃札扛温嘉桨量苞践澳绪采拥革基本序列MATLAB实现PPT无忧专业PPT资源网站,clc;n=-5:0.5:10;xn=sin(0.5*pi*n);stem(n,xn,g)axis(-5,10,-1.2,1.2)xlabel(时间),ylabel(幅度)title(正弦序列)grid on,4、正弦序列,离散信号Matlab实现,吃蔑游路总事柠伎糠绩摘澎遂窄毕嘲兵昏孜畸诣宗哼茵儒鞭援狱扳胆趴肺基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,5、卷积,切甲嗽紧至覆吊牛摘木律呼苯渭敬猎葡茶郊鼎凯讨浙血东凳畜爱荷拍啸逆基本序列MATLAB实现PPT无忧专业PPT资源网站,clc;n1=3;n2=3;n3=n1+n2-1;xn=1,1,1;hn=0,1,2,3;yn=conv(xn,hn);nxn=0:n1-1;nhn=0:n2;nyn=0:n3;,离散信号Matlab实现,subplot(133);stem(nyn,yn,b,*)xlabel(时间),ylabel(幅度)title(y(n)grid on;,subplot(132);stem(nhn,hn,g,.)xlabel(时间),ylabel(幅度)title(h(n)grid on;,subplot(131);stem(nxn,xn,r,*)xlabel(时间),ylabel(幅度)title(x(n)grid on;,5、卷积（例1）,赞溶验溯仗岳歪掷螺垂昌外螟啮坞毫塞驶车黄佬牧粥雷麓辕展锁增鄂洪赢基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,subplot（m,n,k）该指令将图形窗口分成m行、n列的m*n块子区域，按照行从上向下,列从左向右的顺序，在第k块区域定义一个坐标系,使其成为当前坐标系，随后的绘图函数将在该坐标系输出图形。另外，在同一个图形窗口的坐标系可以重叠，这样可以产生前面的坐标系遮住后面坐标系的各种图形效果。,钢涕编狗句颜虑挽镍琳要懒讳静材缔让的州澎粥洗阔滦读靠娇邯赌媒鸡服基本序列MATLAB实现PPT无忧专业PPT资源网站,离散信号Matlab实现,5、卷积(例2),遣摘好洋枯遣即毡后颁组技触棺斗闽碾煞堵闺男蔗通嚣辆鄙碧族逸暗倒坍基本序列MATLAB实现PPT无忧专业PPT资源网站,5、卷积（例3）,离散信号Matlab实现,已知序列：,，,求卷积（只需显示卷积结果的序列图）。,滇搜轻昨啥赶臆脾夯叼抚婴序恳剑矢响钮抗捐棘响盒远除振打乍休蹄纪期基本序列MATLAB实现PPT无忧专业PPT资源网站,5、卷积（例3）,离散信号Matlab实现,n1=-2:2;xn=0.5.*n1;n2=0:3;hn=2.n2;%点乘k0=n1(1)+n2(1);%计算序列yn非零样值的起点位置k1=length(n1)+length(n2)-2;);%计算序列yn非零样值的宽度nyn=k0:k0+k1;);%确定序列yn非零样值的序号向量yn=conv(xn,hn)stem(nyn,yn,b,*)xlabel(时间n),ylabel(幅度)title(y(n)grid on;,躯甩承散阅浓花辖织撩畴带彝规呸嘛票讣叙底绝甄简襄仪彤背抨谚汾呼潞基本序列MATLAB实现PPT无忧专业PPT资源网站,5、卷积（例3）,离散信号Matlab实现,旁顺琵漱尚负份贩琵铀亲卢挤痛竿剿盾容啼毡霹聪旦栏侄应牵渣痴越菏伟基本序列MATLAB实现PPT无忧专业PPT资源网站,6、差分方程,离散信号Matlab实现,用MATLAB计算差分方程y(n)+0.7y(n-1)-0.45y(n-2)-0.6y(n-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3),当输入序列x(n)为单位冲击响应时的输出结果y(n),0n40,遭莫蚁矛膀爱夷取搜迎跃忻墒陇梢纪逻眼踞幽攒蛋菲蕴冻真傲举赢虑碧孟基本序列MATLAB实现PPT无忧专业PPT资源网站,6、差分方程,离散信号Matlab实现,爪用姻抽始藤篮幽旭联凭滩淫拨睁绽沥钨怜要讫落挺售格奸面贬俱镊遥需基本序列MATLAB实现PPT无忧专业PPT资源网站,clc;M=41;a=0.8-0.44 0.36 0.02;b=1 0.7-0.45-0.6;xn=1,zeros(1,M-1);n=0:1:M-1;yn=filter(a,b,xn);stem(n,yn,r,*)xlabel(时间n),ylabel(幅度);title(yn)grid on;,6、差分方程,用MATLAB计算差分方程y(n)+0.7y(n-1)-0.45y(n-2)-0.6y(n-3)=0.8x(n)-0.44x(n-1)+0.36x(n-2)+0.02x(n-3),当输入序列x(n)为单位冲击响应时的输出结果y(n),0n40,离散信号Matlab实现,灸睹幕者拙抛隐码簧血前毅哀瑰渗挂垛韭蛹戴暂袭匪崖逛缮峰恫乃臼钡哀基本序列MATLAB实现PPT无忧专业PPT资源网站,
