MATLAB语言及应用3_zeng.pptVIP

（2） 常微分方程的数值解 matlab为解常微分方程初值问题提供了一组配套齐全,结构严整的指令, 包括: ode45, ode23, ode113, ode23t, ode15s, ode23s, ode23tb。最常用的是ode45，介绍其基本使用方法. ode45使用方法: [x, y]=ode45(‘file_name’, xspan, y0) 输出量：x---微分区间的点系列 y---原函数在微分区间点系列上的函数值 输入量：file_name---文档名,保存的M文件名（该函数文件的输出必须是待解函数的一阶导数，即y’=f(y,x)） xspan---微分区间（一般为二元向量[x0,xf]） y0---初始值 【例】 解微分方程：y’=3x2, y(0)=1,微分区间[0,5] % M文件,保存为dy_file.m function dy=f(x,y) dy=3*x^2; y0=1; xspan=[0,5]; [x, y]=ode45(dy_file,xspan,y0); plot(x,y,r), xlabel(x), ylabel(y) （3）求函数极值点 许多科学研究和工程计算问题都可以归结为一个极值问题,如能量最小，时间最短,最佳拟合,最短路径等.又因为f(x)的极小值问题等价于-f(x)的极大值问题,所以matlab只有处理极小值的指令． 确切地说,这里讨论的只是“局域极值”问题.“全域最小”问题要复杂得多.至今没有一个系统性的方法可求解一般的“全域最小”问题.对于一元二元函数,作图观察对确定全域最小值有很好的应用价值.但更多元的函数,就很难使用作图法. matlab求函数极值的三条指令 [x, fval, exitflag, output]=fminbnd(fun,x1,x2,opt)---求一元函数在区间(x1,x2)中极小值的指令． [x, fval, exitflag, output]=fminsearch(fun,x0,opt)---单纯形法求多元函数极值点的指令．* [x,fval,exitflag,output,grad,hessian]=fminunc(fun,x0,opt)---拟牛顿法求多元函数极值点的指令．* 【例】 求函数y=2x-10x的最小值 (1) 函数采用字符串表达式 func=2^x-10*x; x1=0;x2=8; % 通过观察, 发现最小值在此区间 [Xmin,Ymin]=fminbnd(func,x1,x2) Xmin = 3.8507 Ymin = -24.0800 (2) 函数采用M文件形式 % M文件内容，保存为func.m function y=f(x) y=2^x-10*x; % M文件内容 x1=0; x2=8; [Xmin,Ymin]=fminbnd(func,x1,x2) xx=0:0.01:8; yy=2.^xx-10.*xx; plot(xx,yy) 4 信号处理和系统分析举例——数字滤波 y=filter(B, A, x)---运用数字滤波对输入信号x滤波 滤波器的作用:通过对采样数据信号进行数学运算处理来达到频域滤波的目的. 例 设计一个低通滤波器,从受噪声干扰的多频率混合信号x(t) 中获取10HZ的信号. x(t)=sin(2π*10*t)+cos(2π*100*t) +n(t). randn(state,1); ws=1000; % 采样频率 t=0:1/ws:0.4; % 产生带噪声的多频率信号 x=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t)); wn=ws/2;fc=30; % Nyquest频率wn,绝对截止频率fc [B,A]=butter(10,30/wn); % 产生截止频率为fc/wn的10阶Butterworth低通滤波器分子 、分母多项式B和A，作为下面滤波处理的输入 y=filter(B, A, x); % 滤波处理 plot(t,x,b-,t,y,r.,MarkerSize,10); Poly2sym P306 把向量多项式转换为符号多项式 Polyval P171 Poly P171 ：“根据多项式根求多项式各项系数”的操作，输入向量的各元素表示多项式的根，输出量是多项式的系数向量 Class P298 拟合问题 P173 插值问题 “精通6.5”P158 [interp1:一维插值 参数见表 或见本ppt P14] 欧姆定理U=I*R：斜率为R的直线 插值方式：(数字图像插值算法