作业帮matlab里循环求积分作图的问题.
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/15 06:28:45
matlab里循环求积分作图的问题.
我现在用matlab作出了一张这样的图~（这个图是用excel画的.比较好看.matlab的那个一样的）
曲线的函数已知,例如ii档的函数是
1/a=1/( -0,284300594+0,30105565*u-0,019711833*u^2+0,000567339*u^3-6,22615*10^(-6)*u^4)
就是图中的x轴是u,y轴是1/a
现在要求曲线下面积的积分,t=∫1/( -0,284300594+0,30105565*u-0,019711833*u^2+0,000567339*u^3-6,22615*10^(-6)*u^4) du , 但是u要从5到34并以0.001为增量,u每一次改变求出相应的t然后作t-u图,在matlab里面要怎么编啊.现在编出来的就是直接求出面积的.做不出图来哇.
我现在用matlab作出了一张这样的图~（这个图是用excel画的.比较好看.matlab的那个一样的）
曲线的函数已知,例如ii档的函数是
1/a=1/( -0,284300594+0,30105565*u-0,019711833*u^2+0,000567339*u^3-6,22615*10^(-6)*u^4)
就是图中的x轴是u,y轴是1/a
现在要求曲线下面积的积分,t=∫1/( -0,284300594+0,30105565*u-0,019711833*u^2+0,000567339*u^3-6,22615*10^(-6)*u^4) du , 但是u要从5到34并以0.001为增量,u每一次改变求出相应的t然后作t-u图,在matlab里面要怎么编啊.现在编出来的就是直接求出面积的.做不出图来哇.
clear;clc
u=5:0.001:34;
inva=1./( -0.284300594+0.30105565*u-0.019711833*u.^2+0.000567339*u.^3-6.22615*10^(-6)*u.^4);
t=cumtrapz(u,inva);
[a,h1,h2]=plotyy(u,inva,u,t);
legend([h1,h2],'II档','II档积分')
u=5:0.001:34;
inva=1./( -0.284300594+0.30105565*u-0.019711833*u.^2+0.000567339*u.^3-6.22615*10^(-6)*u.^4);
t=cumtrapz(u,inva);
[a,h1,h2]=plotyy(u,inva,u,t);
legend([h1,h2],'II档','II档积分')