matlab语言设计介绍第五讲(Revision).ppt

第5章 MATLAB基本应用领域;Section 1:线性代数问题;Section 1:线性代数问题;2.线性代数方程求解;3.数据拟合;T=[0:.1:2.5]; Y=[ones(size(T)) exp(-T)]*c; plot(T,Y,-,t,y,o) title(最小二乘法曲线拟合) xlabel(\itt),ylabel(\ity);4.矩阵求逆、求幂和矩阵指数;多项式表示 多项式: 系数矢量p=[1 0 -2 -5] 多项式的根 p(x)=0的解: r=roots(p) r’ = 2.0946 -1.0473 + 1.1359i -1.0473 - 1.1359i 多项式计算 p(4)的值为: y=polyval(p,4) y= 51 ;矢量的卷积和去卷积;Section 2:多项式与内插;插值所得曲线要通过节点！ 另外，还有基于FFT的内插和二维内插;max;xstd=std(x) xstd = 0.8685 y=randn(1,100); xycov=cov(x,y) xycov = 0.7543 -0.0146 -0.0146 0.8924 xycorr=corrcoef(x,y) xycorr = 1.0000 -0.0178 -0.0178 1.0000;2.数据傅里叶分析和FFT;;问题1：找到最大谱峰对应的频点？ 问题2：将20Hz 的低频分量保留，其余频率分量 的信号滤除？ 问题3: 如何将其余两个局部最大值对应的频点找 到？ ;泛函：MATLAB提供可对函数进行操作的函数;函数可直接编写成M文件，例如 function y=humps(x) y=1./ ((x-0.3).^2+0.01)+1./((x-0.9).^2+0.04)-6; 为绘制出这函数的曲线，可输入 fplot(humps,[-5,5]); xlabel(‘\itx’); ylabel(‘\ity’);grid on; xmin=fmin(humps,0.3,1) xmin = 0.6370 xzero=fzero(humps,[-1 0]) xzero = -0.1316 ;数值积分;微分方程求解过程 将常微分方程转换成一阶微分方程组,即表示成右函数形式,这是利用龙格－库塔法求解微分方程的前提；;列车以20m／s速度行驶，当制动时可获得-0.4米/平方秒的 加速度，求使列车停下的制动时间和路程。 微分方程： 由此得一阶微分方程组： ;X0=[0;20]; tspan=[0:48,49:0.01:51]; [T,X]=ode45(diffequat1,tspan,X0); N=find(X(:,2)=0); t=T(N(1)),s=X(N(1),1),v=X(N(1),2) figure(1) subplot(2,1,1),plot(T,X(:,1),r),title(列车滑动路程), grid on xlabel( t(s));ylabel( s(m)); subplot(2,1,2),plot(T,X(:,2),k),title(列车运行速度), grid on xlabel( t(s));ylabel(v(m/s));;Section 5:常微分方程求解;第5章习题