matlab教学ppt第4章.ppt

第四章、数值计算; 内容提要; 4.1矩阵的相关计算;例题开讲; 4.2 线性方程组的解;4.2 线性方程组的解;求解下列线性方程组：（矩阵求逆法）; 4.3 多项式计算; 4.3多项式计算;例题开讲; 4.3 多项式计算;例题开讲; 4.3 多项式计算;4.3 多项式计算; 4.3 多项式计算;1、矩阵的相关计算 2、线性方程组的求解 3、多项式和卷积; 4.3 多项式计算; 4.4 拟合与插值;4.4 拟合与插值;4.4 拟合与插值; 4.4 拟合与插值;4.4 拟合与插值;4.4 拟合与插值; 4.4 拟合与插值; 4.4 拟合与插值;给出概率积分 的数据如下表所示：;4.5 matlab的泛函指令;4.5.1 数值积分;4.5.2 常微分方程求解;%(2)根据一阶微分方程组编写M函数文件 [dydt.m] function ydot=dydt(t,y) mu=2; ydot=[y(2);mu*(1-y(1)^2)*y(2)-y(1)];;例题开讲; 4.6 符号计算;4.6 符号计算;4.6 符号计算;例题开讲;例二： pi1=sym(‘pi’); k1=sym(‘8’); k2=sym(‘3’); pi2=pi; r1=8; r2=3; sin(pi1/3) sin(pi2/3) sqrt(k1+sqrt(k2)) sqrt(r1+sqrt(r2)) ;例题开讲;4.6 符号计算;4.6 符号计算; 4.6 符号计算; 4.6 符号计算;4.5.1 求函数的零点;P1=0.1; P2=0.5; y_C=sin(x).^2.*exp(-P1*x)-P2*abs(x); x=-10:0.01:10; Y=eval(y_C); plot(x,Y,r); hold on plot(x,zeros(size(x)),k); xlabel(t); ylabel(y(t)) hold off zoom on [tt,yy]=ginput(5); zoom off tt [t4,y4,exitflag]=fzero(y_C,tt(4),[],P1,P2);%（2）作图法观察函数零点分布 x=-10:0.01:10; %对自变量采样，采样步长不宜太大。 Y=eval(y_C); %在采样点上计算函数 plot(x,Y,r);hold on,plot(x,zeros(size(x)),k); xlabel(t);ylabel(y(t)),hold off