![]() |
[求助]如何求值一表达式
2 个附件
我不知道能不能在这里显示公式,试了一下没有得行,所以只好以附件的形式传上。本来想这样已麻烦大家了,还要下附件来看就有点不好意思了:)。不过,本人初学,而有点急于求此式值,还请大家给我指点。第一,这样求对不对?有没有可能简便?;第二,结果输出也太长了,呵呵 不晓得如何化简,其中的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) |
回复: [求助]如何求值一表达式
没有人回复啊? 帮我看看啊,各位达人
|
所有时间均为北京时间。现在的时间是 14:52。 |
Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.