没有合适的资源?快使用搜索试试~ 我知道了~
MATLAB数值积分与微分.pdf
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 160 浏览量
2021-09-14
15:47:48
上传
评论
收藏 138KB PDF 举报
温馨提示
试读
3页
。。。
资源推荐
资源详情
资源评论
.
1. 数值微分
MATLAB 没有直接求数值导数的函数,diff 函数可求数
值微分,其调用格式为:
DX=diff(X):DX(i)=X(i+1)-X(i),i=1,2,…,n-1。
DX=diff(X,n):n 阶
DX=diff(A,n,dim):计算矩阵 A,dim=1 时(缺省状
态),按列;dim=2,按行。
V=vander(1:4)
DV=diff(V)
gradient 函数可按行或列求矩阵的数值微分
[dx,dy]=gradient(V) %除了第一行(列)与最末一行(列)按前向公式,其余的
dx=(y(:,i+1)-y(:,i-1))/2; dy=(y(i+1,:)-y(i-1,:))/2;
gradient(V)
编程实现以上求一阶或二阶近似的算法(diffc.m)
h=0.05;x=0:h:pi; y=sin(x)./(x+cos(x));
y1=diffc(y,h,1,2); %n=1
y2=diffc(y,h,2,2); %n=2
y3=diffc(y,h,3,2); %n=3
y4=diffc(y,h,4,2); %n=4
subplot(2,2,1);plot(x,y,x(1:length(x)-4),y1);
subplot(2,2,2);plot(x,y,x(1:length(x)-
4),y2);subplot(2,2,3);plot(x,y,x(1:length(x)-6),y3);
subplot(2,2,4);plot(x,y,x(1:length(x)-6),y4);
这种方法比循环法简单且效率高。
用不同的方法求函数 f(x)的数值导数,并在同一个坐标系
中做出 f'(x)的图像。
程序如下:
f=inline('sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+5*x+2');
g=inline('(3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-
x+12)/2+1/6./(x+5).^(5/6)+5');
x=-3:0.01:3;
p=polyfit(x,f(x),5); %用 5 次多项式 p 拟合 f(x)
.
资源评论
苦茶子12138
- 粉丝: 1w+
- 资源: 6万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功