随便填一下.docVIP

MATLAB数值微分自学报告 数值微分，也称数值求导，是以函数的离散值求某点函数的导数值，通常应用于函数以离散数值形式给出的情况。应用MATLAB,可进行导数、偏导数、梯度、矩阵等计算。 MATLAB可以以符号形式计算一元函数的微分、多元函数偏微分、全微分等，调用函数diff.m，调用函数如下： 符号求导的命令功能 yx=diff(f(x),x) 求函数 y = f (x)对x的一阶导函数 y= f (x) dy=diff (f(x),x)*dx 求函数 y = f (x)的一阶微分dy = f (x)dx yxx=diff(f(x),x,2) 或 yxx=diff(yx,x) 求函数 y = f (x)对x的二阶导函数 y’’ = f x yxxx=diff(f(x),x,3) 或 yxxx=diff(yxx,x) 求函数 y = f (x)对x的三阶导函数y’’’ = f x yxn=diff (f(x),x,n) 求函数 y = f (x)对x的n阶导函数 y(n) = f (n) (x) yn=simple(yxn) 将n阶导函数 y(n) = f (n) (x)化简，并记作yn pretty(diff(f(x),x) 输出一个符合日常书写习惯的一阶导函数的表达式 (偏微分函数调用与上表近似，易查知，略。) 函数以离散数值形式给出时，可以： 一、一阶导数的数值计算及其MATLAB程序 利用差商求导、中心差商公式求导 程序如下： function [n,xi,yx,wuc]=sandian(h,xi,fi,M) n=length(fi); yx=zeros(1,n); wuc=zeros(1,n); x1= xi(1); x2= xi(2); x3= xi(3); y1=fi(1); y2=fi(2); y3=fi(3); xn= xi(n); xn1= xi(n-1); xn2= xi(n-2); yn=fi(n); yn1=fi(n-1); yn2=fi(n-2); for k=2:n-1 yx(1)=(-3*y1+4*y2-y3)/(2*h); yx(n)=(yn2-4*yn1+3*yn)/(2*h); yx(2)=( fi(3)- fi(1))/(2*h); yx(k)=( fi(k+1)- fi(k-1))./(2*h); wuc(1)=abs(h.^2.*M./3); wuc(n)=abs(h.^2.*M./3); wuc(2:n-1)=abs(-h.^2.*M./6); 理查森（Richardson）外推法求导 此方法利用矩阵对导函数值进行近似、公式不断推进矩阵直至满足精度要求，输出满足条件导函数及其对应矩阵的坐标，MATLAB程序如下： function [Dy,dy,jdw,n]=rendiffext1(fun,x0,jdwc,max1) h=1;j=1; n=1;jdW=1;xdW=1; x1=x0+h;x2=x0-h; Dy(1,1)=(feval(fun,x1)- feval(fun,x2))/(2*h); while((jdWjdwc)(jmax1)) j;x1=x0+2^(-j)*h;x2=x0-2^(-j)*h; Dy(j+1,1)=(feval(fun,x1)-feval(fun,x2))/(2^(1-j)*h); for k=1:j k;Dy(j+1,k+1)= Dy(j+1,k)+( Dy(j+1,k)- Dy(j,k))/(4^k-1); end jdW=abs(Dy(j+1,j+1)-Dy(j+1,j)); j=j+1; end [n,n]=size(Dy);jdw=abs(Dy(n,n)-Dy(n,n-1)); dy= Dy(n,n); 用精度为O(h4 )的中心差商公式计算f (x)的近似值和误差估计的MATLAB程序 function [x0,yx,wuc]=zxcs4(h,x0,fi,M) xi=[x0-2*h,x0-h,x0,x0+h,x0+2*h]; x1= xi(1); x2= xi(2); x3= xi(3); x4= xi(4); x5= xi(5);y1=fi(1); y2=fi(2); y3=fi(3); y4=fi(4); y5=fi(5); yx=(8*y4-8*y2-y5+y1)/(12*h); wuc=abs(h.^4*M/30); 变步长的中心差商公式及其MATLAB程序 输入的量：fun是函数 y = f (x)，x0是导数点, wu是 f (x)的近似值的精度, max1是计算次数的最大值，h0是初始步长。 输出的量：H表示步长h /10k 组成的数组，Dy 是导数值组成的数组，W 是误差限