Matlab基础培训XX578.docx
Mat1ab基础培训XX5785.7上午：1、he1pe1mat基本矩阵生成函数2、doc中的例子可直接按F9运行，不需要复制3、矩阵合并:C=AB&A;B4、统计有关函数he1pdatafun矩阵的统计是对列操作，与向量的统计是完全不一致的5、数学操作函数he1pe1fun数学操作是对向量2或者者矩阵的每个元素进行操作6、Size1engthnume1(矩阵元素个数)7、随机数分布函数在统计工具箱8、分段操作%符的使用ev1uatece11,能够公布fi1e里PUb1iCatiOn功能ce11mode9、he1pfi1eformatsmata1b支持的文件格式，操作函数可查看he1poifuntextscan-Readformatteddatafromtextfi1e.Textread复合数据格式中的说明文字会使用ce11数组存储10.ce11元包数组的生成a=1abc'Jnt(8b1,使用的是大括号。读取的方式同样使用大括号。b=a1,1;若是b=a(1,1)则产生新的Ce1Ia1TayIOx日期三种格式datastr字符串datanum数字整数datavec向量(1*6)weekday»×=datenum('1203')X=731583»x=datestr('1203')X=02-Jan-20031oadbe1gmarriagedv=datavec(dates)5.7下午：1、在图中可现实特殊的字符SPrimf(，n(字符)')，可参见附录A2、datetick函数x1imy1imeditausmarriageaxes中能够看到效果gtext(")第七章联系最后的editboxoffice13、逻辑矩阵与逻辑索引在寻找特殊值时很有用1oadgasprices.matid×=(Germany<3.5);Year(idx)可找到油价小于3.5的年份4逻辑操作函数关键字可查看iskeyword,编写流程操纵语句等参见例子editife1seeditfor1opeditwhi1e1oop5、友好界面输入输出Dispfprintf写文件SPrintf写字符串InputmenuEditnumguess×phi,yphi=ginpu(17)t让你选用多个数据点，给你返回坐标值：»p1ot(peaks)»xphizyphi=ginput(3)6、三维空间曲面的绘制SUrfimagescimageCo1ormap:a=jet(8)jet(16)bone(8)Qnd(8,3);imagesc(a)Contour绘制二维的等高线图还有Mesh绘制的网格，但网格不上色Waterfa11ribbonSpysurf1等editmembranep1ots例子7、COVxcorrcoef计算方差协方差有关系数8、多项式拟合p1oyfit(参数1,参数2,阶数)，拟合后可用过P1oyva1计算位置点值(预测)一维插值interp1二维插值interp29、平滑与卷积ConVnfi1ter等10、线性自回归模型fftgui演示M文件加密使用pcode生成.p文件写字板打开会是乱码5.8上午：1、m文件调试，fix小工具可帮户迅速定位与修改有错误的代码。但不能完全依靠此。2、调试断点中可使用设置断点的条件(右键设置n=8)对循环做多次时才停止，否则常规的会在第一次循环时就就停止，editfor1oop3、ProfiIer里调试函数时间。可全面给出每行代码执行的时间与次数等信息，有助于快速定位需要优化的代码部分。假如只要单独显示程序运行总时间可使用tic,toco4、代码中需要对矩阵预先分配内存可有效的提高程序的效率。把标量的操作变为向量的操作也可有效的提高效率editmakeA3,meshgrid()函数生成矩阵在对矩阵行列操作，其中能够看到函数只调用了一次。»tic;A=InakeA1(500,500);toe原始循环E1apsedtimeis0.833463seconds.»tic;A=InakeA2(500,500);toe预分酉已E1apsedtimeis0.201760seconds.»tic;A=makeA3(500z500);toe向量化操作E1apsedtimeis0.011711seconds.即尽量避免使用for循环。5、数据类型：numberic(有符号与无符号各类位长，与Sing1eanddOUbIe).1ogic.ce11,structure.Functionhand1e,userc1ass.6、函数句柄的使用：h=(X,Y)XA+YA2;h(1,2).短小的计算可使用比编写m执行效率高。作为工具箱里的函数句柄使用：比如求解器里的。de45等务必以句柄的形式传递。可参见ode45he1p文档。优化工具箱里的很多指令都务必使用函数句柄。信号处理1、SoUnd听取信号的声音。StriPS可每个多久绘制图形长数据的清晰完整的显示。fs=1000;ts=0:1/fs:2;x=250+240*sin(2*pi*fs.*ts);StriPS(X,0.25,fs);sound(xzfs);%fs=10000;t=0:1/fs:1.5;x=sawtooth(2*pi*50*t,1/2);p1ot(t,x),axis(00.2-11)工具箱Sptoo1可对信号做全面的查看与分析2、editnoisyC可听噪声的效果。3、重新采样函数;downsanp1eupsamp1eresamp1einterpdecimate(抽取，类似降采样，实际上时通过一个低通滤波)，具体可参见he1p文档。4、XCOrr计算有关函数，XCOV计算协方差，conv计算卷积5、频谱分析：采样信号的频谱是对原始信号频谱的周期严拓，分析时只需要对-FS/2到FS/2或者者0到FS分析；Fftshift可对绘制的横轴频率范围进行调整。a.加窗可降低旁瓣的影响。旁瓣会平滑，但代价是主瓣会展宽。EditWindft窗函数Wintoo1工具箱，产生的窗可储存至WOrkSPace；常见的窗函数也可使用命令行使用，docsigna1/windowb.zeropadding也可用来提高频谱的细节：edittestzeropad可查看信号处理工具箱的介绍。C还可使用频谱对象来作频谱分析.Docpsd查看2、时变谱SPetCrOgram短时傅里叶分析。»1oadhande1»sound(y,Fs);»SPeCtrOgram(y,kaiser(100,5),75,512,Fs)例子specgramdemoxpsoundHA190003、线性时不变系统1T1差分方程Z变换zzp,k=Zftpk(bza)把传递函数转换成零极点形式。Zp1ane(z,p)看零极点图各类形式的转换：1inearsystemtransformation冲击响应b=120a=1O0.2impz(b,a,10,1e3)10个点，采样频率IooO例子：edittransferp1ot可查看三维零极点图。得到零极点后可在SPtoO1里涉及滤波可在此图中得到群时延,滤波器系数等。也可利用这个滤波器来实现信号通过滤波器的效果。滤波器的指令y=fi1ter(bza,x)o另外可使用离散滤波器对象。类似前面的频谱对象h。dfi1t.strucuredocdfi1tdocfvtoo1.3、HR滤波器。A使用命令行的方式实现：EditIirtypes文件里具体查看可运行五种不一致滤波器的频率响应。BSPtOOI里的fdatoo1工具箱来设计。设计好的滤波器可使用F11E下的generateMfi1e生成m文件代码。生成的都是滤波器对象而不是滤波器系数。DocfdesignC.任意响应滤波器设计yu1ewa1k（n阶数，f频率点,m幅度响应）EditdIrectstopdirectstop（10）directstop（20）4、FIR滤波器设计。相位线性好，但是阶数高，延时大要紧方法是窗函数的方法。fir1（标准的四种）fir2（任意响应滤波器）fir1s（多通带滤波器，优化的方法而非窗函数）firpm（多带滤波器，求解最大最小优化误差）firc1sfircos（升余弦滤波器）firc1s1b=0.4sinc();fvtoo1(b,1);Windemo(20)editwindemoWindemo(20)FteQX<K>>能够看出阻带通带的变化。editfir1demoeditdirectstop2;directstop2(20);PoIyStab(b)把线性相位(冲击响应关于中心店对称)滤波器转换成最小相位滤波器editMinde1ay可看出线性相位与最小相位在各方面的区别6、频域滤波长数据阶段然后fft再卷积在ifft最后再叠加各段数据。Editfi1ttimesAxedSigna11ength滤波器长度增大，运行的时间不可能有很大的增加。
