Matlab基础培训XX578.docx

《Matlab基础培训XX578.docx》由会员分享，可在线阅读，更多相关《Matlab基础培训XX578.docx（7页珍藏版）》请在第一文库网上搜索。

1、Matlab基础培训XX5785.7上午:1、help elmat基本矩阵生成函数2、doc中的例子可直接按F9运行，不需要复制3、矩阵合并:C=A B & A;B4、统计有关函数help datafun矩阵的统计是对列操作，与向量的统计是完全不一致的5、数学操作函数help elfun数学操作是对向量2或者者矩阵的每个元素进行操作6、Size length numel(矩阵元素个数)7、随机数分布函数在统计工具箱8、分段操作%符的使用evluate cell,能够公布file里publication功能 cell mode9、help fileformats matalb支持的文件格式，操作

2、函数可查看help oifun textscan - Read formatted data from text file.Textread复合数据格式中的说明文字会使用cell数组存储10cell元包数组的生成3=1,匕13d加.8),1,使用的是大括号。读取的方式同样使用大括号。 b=al,l;若是 b=a (1,1)则产生新的 Cellarray10、 日期 三种格式datastr字符串datanum数字整数 datavec向量(1*6) weekday x=datenum(,l203) X =731583 x=datestr(,l203) X =02-Jan-2003 load bel

3、gmarriage dv=datavec(dates)5.7下午:1、在图中可现实特殊的字符SPrintf Cn (字符)，可参见附录A2、datetick函数xlim ylim edit ausmarriageaxes中能够看到效果gtext()第七章联系最 后的 edit boxofficel3、逻辑矩阵与逻辑索引在寻找特殊值时很有用Load gasprices.matidx=(Germany corrcoef计算方差协方差有关系数8、多项式拟合PIoyfit (参数1,参数2,阶数)，拟合后可用过PIoyVal计算位置点值(预测)一维插值interpl二维插值interp29、平滑与卷积

4、ConVn filter等10、线性自回归模型fftgui演示M文件加密使用pcode生成.p文件 写字板打开会是乱码5.8上午：1、m文件调试，fix小工具可帮户迅速定位与修改有错误的代码。但不能完全依靠此。2、调试断点中可使用设置断点的条件(右键设置n=8)对循环做多次时才停止，否则常规 的会在第一次循环时就就停止，edit forloop3、PrOfiIer里调试函数时间。可全面给出每行代码执行的时间与次数等信息，有助于快速定 位需要优化的代码部分。假如只要单独显示程序运行总时间可使用tic, tOCo 、代码中需要对矩阵预先分配内存可有效的提高程序的效率。把标量的操作变为向量的操作也可

5、有效的提高效率edit makeA3 , mShgrid ()函数生成 矩阵在对矩阵行列操作，其中能够看到函数只调用了一次。 tic; A=makAl (500,500) ; toe原始循环Elapsed tim is 0.833463 SCOnds. tic; A=makA2 (500,500) ; toe预分酉己EIaPSd tim is 0.201760 SCOnds. tic; A=makA3 (500,500) ; toe向量化操作EIaPSd tim is 0.011711 seconds .即尽量避免使用for循环。5、数据类型：numberic(有符号与无符号各类位长,与 Si

6、ngIeanddoUble). logic, cell, structure.Function handle, user class.6、函数句柄的使用：h=(X,Y)X八+Y八2; h(l,2).短小的计算可使用比编写m执行效率高。 作为工具箱里的函数句柄使用：比如求解器里的。de45等务必以句柄的形式传递。可参 见ode45help文档。优化工具箱里的很多指令都务必使用函数句柄。信号处理1、SOUnd听取信号的声音。StriPS可每个多久绘制图形长数据的清晰完整的显示。fs=1000;ts=0:l/fs:2;x=250+240*sin(2*pi*fs.*ts);Strips(x,0.25,

7、fs);sound(x,fs);%fs = 10000;t = 0:1/fs:1. 5;X = SaWtOOth(2*pi*50*t, 1/2);plot (t, x), axis (0 0. 2 -1 1)工具箱Sptool可对信号做全面的查看与分析2、edit noisyC可听噪声的效果。3、重新采样函数;downsanple upsample resample interp decimate(抽取，类似降采样, 实际上时通过一个低通滤波)，具体可参见help文档。4、XColT计算有关函数，XCoV计算协方差，corn/计算卷积5、频谱分析：采样信号的频谱是对原始信号频谱的周期严拓，分析

8、时只需要对-FS/2到FS/2或者者0 到FS分析；Fftshift可对绘制的横轴频率范围进行调整。a.加窗可降低旁瓣的影响。旁瓣会平滑，但代价是主瓣会展宽。EditwindftTruncated, Windowed SignalTruncated SignalDFT of Truncated SignalDFT of Truncated, Windowed Signal 25窗函数WintoOl工具箱，产生的窗可储存至WorkSPace；常见的窗函数也可使用命令行使 用，doc signal/windowb. zero padding也可用来提高频谱的细节：edit testzeropad可

9、查看信号处理工具箱的介绍。4096 point FFT with 1024 data points注意区分尖锋的区分，要紧原因在于点间的分辨率。Fft 的演试 demo 有： demol palyshow fftdemo phoneplayshow sunsppots5.8下午1、功率谱密度估计：A命令行方式子空间：EV (peig) MuSIC(Pmusic)分辨率更高，计算量大mF参数：周期图 perioddigram Welch multipaper (pmtm) 参数法：pyulear,pburg, pcov,B使用sptool里的工具箱。FS=IOOO;t=0:l/FS:2-l/FS

10、; x=sin(2*pi*200*t)+0.01*sin(2*pi*300*t); sptoolC还可使用频谱对象来作频谱分析.Doc psd查看2、时变谱SPetCrOgram 短时傅里叶分析。load handel sound(y,Fs); SPeCtrOgram(y,kaiser(100,5),75,512,Fs)例子 Specgramdemo xpsound HAL90003、线性时不变系统LTl差分方程Z变换z,p,k=Zftpk (b,a)把传递函数转换成零极点形式。Zplane (z,p)看零极点图各类形式的转换： Linear system transformation冲击响应

11、b=12 0a=10 0.2impz(b,a,10,le3) 10 个点，采样频率 IOoOFreqz (b,a)例子：edit transferplot可查看三维零极点图。得到零极点后可在SPtOol里涉及滤波 器。可在此图中得到群时延,滤波器系数等。也可利用这个滤波器来实现信号通过滤波器的效果。滤波器的指令y=filter(b,a,x)0另外可使用离散滤波器对象。类似前面的频谱对象h。dfilt.strucure doc dfiltdoc fvtool.3、IlR滤波器。A使用命令行的方式实现：EditlirtyPeS文件里具体查看可运行五种不一致滤波器的频率响应。B sptool里的fd

12、atool工具箱来设计。设计好的滤波器可使用FILE下的generate Mfile生成m文件代码。生成的都是滤波器对象而不是滤波器系数。Doc fdesignC.任意响应滤波器设计yulewalk（n阶数，f频率点,m幅度响应）Edit directstop directstop （10）directstop （20）Comparison of Response Magnitudes1.41.2p2c6ehpu6eiiWindemo(20)4、FIR滤波器设计。相位线性好，但是阶数高，延时大要紧方法是窗函数的方法。31（标准的四种）2（任意响应滤波器）firls（多通带滤波器，优化的方法而非

13、窗函数） firpm（多带滤波器，求解最大最小优化误差）firclsfircos（升余弦滤波器）firclslb=0.4sinc();fvtool(b,l);edit windemo Windemo(20)能够看出阻带通带的变化。edit firldemoedit directstop2;directstop2(20);PolyStab(b)把线性相位(冲击响应关于中心店对称)滤波器转换成最小相位滤波器 edit Mindelay可看出线性相位与最小相位在各方面的区别6、频域滤波长数据阶段然后fft再卷积在ifft最后再叠加各段数据。Edit filttimes(Sos)u6-1-LL。1 EP滤波器长度增大，运行的时间不可能有很大的增加。