Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
 
 
主题工具 显示模式
旧 2008-08-13, 13:54   #1
xiaoduan259
初级会员
 
注册日期: 2008-07-12
年龄: 44
帖子: 13
声望力: 17
xiaoduan259 正向着好的方向发展
默认 [求助]如何求值一表达式

我不知道能不能在这里显示公式,试了一下没有得行,所以只好以附件的形式传上。本来想这样已麻烦大家了,还要下附件来看就有点不好意思了:)。不过,本人初学,而有点急于求此式值,还请大家给我指点。第一,这样求对不对?有没有可能简便?;第二,结果输出也太长了,呵呵 不晓得如何化简,其中的Struveh怎么处理?谢谢!
function lpegrowth()
a=0;c=100;b=0.1;
C1=12;Cs=0.013158;D=0.0004;K=0.01;
t=30;
x=[a:b:c];
y=besselj(3,x);
i=find(y(1 : (c-a)/b).*y(2 : (c-a)/b+1)<0);
for j=1:size(i,2) %size(i,2)
r(j) = fzero(inline('besselj(3,x)'),[x(i(j)) x(i(j)+1)]);
end
for j=1:size(i,2)
syms x t
A(j)=(int(besselj(3,x).*x,0,r(j)))/(besselj(4,r(j))).^2;
F(j)=A(j)/r(j).^2*(besselj(4,r(j)/2)+besselj(2,r(j)/2));%取l/L=1/2
sum=symsum(F(j),1,size(i,2));
h=249.85-901.6*sum;%将各参数代入化简
end
hh=simple(h)

h=
4997/20-423257136495541022704583672537678/188467143083853026902895176025*besselj(0,1/9903520314283042199192993792*944461883059960475298757144727590783263686620542320351641600^(1/2))+306472754990119802/6507108870716640195622239493166121669178311875*944461883059960475298757144727590783263686620542320351641600^(1/2)*besselj(1,1/9903520314283042199192993792*944461883059960475298757144727590783263686620542320351641600^(1/2))+634885704743311534056875508806517/753868572335412107611580704100*pi*(besselj(0,1/9903520314283042199192993792*944461883059960475298757144727590783263686620542320351641600^(1/2))*StruveH(1,1/9903520314283042199192993792*944461883059960475298757144727590783263686620542320351641600^(1/2))-besselj(1,1/9903520314283042199192993792*944461883059960475298757144727590783263686620542320351641600^(1/2))*StruveH(0,1/9903520314283042199192993792*944461883059960475298757144727590783263686620542320351641600^(1/2)))


hh=
4997/20-423257136495541022704583672537678/188467143083853026902895176025*besselj(0,6905297936013695/70368744177664)+43132205786647276837666665005056/942335715419265134514475880125*besselj(1,6905297936013695/70368744177664)+634885704743311534056875508806517/753868572335412107611580704100*pi*besselj(0,6905297936013695/70368744177664)*StruveH(1,6905297936013695/70368744177664)-634885704743311534056875508806517/753868572335412107611580704100*pi*besselj(1,6905297936013695/70368744177664)*StruveH(0,6905297936013695/70368744177664)
上传的附件
文件类型: doc 计算.doc (19.0 KB, 3 次查看)
xiaoduan259 当前离线   回复时引用此帖
 


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 05:51


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.