[问题] 帮忙求积分 - MATLAB爱好者论坛-LabFans.com

tannengzhi0504 · 2010-12-15, 14:37

clear;clc
syms ds;
p0=6*10^6;
v0=50*10^-3;
p2=16*10^6;
m=30;
n=1.4;
A=pi*0.125^2/4;
t=sqrt([2*ds*m*[p0-(3*ds*A/2*v0)^n]]/(p0*p2*A));
s=int(t,ds,0,0.53);
ss=vpa(s);
0.0044965893291238346267263438504854
不知道对不对，你参考一下

justin_mao · 2010-12-15, 21:51

数值积分可以直接用quadl的呀···

代码:

p0=6*10^6;
v0=50*10^-3;
p2=16*10^6;
m=30;
n=1.4;
A=pi*0.125^2/4;
t=@(ds) sqrt([2.*ds.*m.*[p0-(3.*ds.*A./2.*v0).^n]]/(p0.*p2.*A));
quadl(t,0,0.53)

应该可以的

2010-12-15, 14:37	#1
tannengzhi0504 初级会员注册日期: 2010-05-31 帖子: 21 声望力: 15	回复: 帮忙求积分 clear;clc syms ds; p0=610^6; v0=5010^-3; p2=1610^6; m=30; n=1.4; A=pi0.125^2/4; t=sqrt([2dsm[p0-(3dsA/2v0)^n]]/(p0p2A)); s=int(t,ds,0,0.53); ss=vpa(s); 0.0044965893291238346267263438504854 不知道对不对，你参考一下

2010-12-15, 21:51	#2
justin_mao 初级会员注册日期: 2010-12-07 帖子: 2 声望力: 0	回复: 帮忙求积分数值积分可以直接用quadl的呀··· 代码: p0=610^6; v0=5010^-3; p2=1610^6; m=30; n=1.4; A=pi0.125^2/4; t=@(ds) sqrt([2.ds.m.[p0-(3.ds.A./2.v0).^n]]/(p0.p2.A)); quadl(t,0,0.53) 应该可以的