那位大虾帮帮忙阿？

hanxiao · 发表于 2006-3-14 16:29

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

这是我的程序，很简单的，都是矩阵的加碱，积分，求导什么的。为什么一算就很长时间，然后自动退出呢？第一个积分命令int前的程序都能运行，到int就不行了，怎么回事呢？那位大虾帮忙看看，谢谢了先！！ syms z a=[z^2/34^2,z^3/34^3,z^4/34^4] a1=diff(a) a2=diff(a1) m=int((212160-858*z)*a'*a,0,34) kesi=int((212160-858*z)*a',0,34) p=-kesi*2.383*10^-2 r=int(2.1*10^11*(29.962-0.702*z)*a2'*a2,0,34) k2=-int(10^6*a1'*a1,0,34) k=r+k2 syms t e=k-t^2*m d=-0.1*k i=d*d+e*e yi=inv(i)*d*p yr=inv(i)*e*p x=yr(1)+yr(2)+yr(3) y=yi(1)+yi(2)+yi(3) w=x^2+y^2 b1=int(2*w,t,0,8) b2=int(2*t^2*w,t,0,8) double(b1) double(b2) z=34 a2 =[1/578,3/19652*z,3/334084*z^2] x1=a2(1)*yr(1)+a2(2)*yr(2)+a2(3)*yr(3) y1=a2(1)*yi(1)+a2(2)*yi(2)+a2(3)*yi(3) w1=x1^2+y1^2 wj1=12.8*10^11*int(2*w1,t,0,8) wj2=12.8*10^11*int(2*t^2*w1,t,0,8) double(wj1) double(wj2)

hanxiao · 发表于 2006-3-14 19:54

并求一个解决方法。不胜感激。。。

happy · 发表于 2006-3-15 09:19

<DIV class=quote>以下是引用hanxiao在2006-3-14 16:29:44的发言： kesi=int((212160-858*z)*a',0,34) </DIV>
这个积分不存在，换数值积分吧

kevinhcj · 发表于 2006-3-15 15:59

我倒是很快能算出kesi来，但是bl,b2较长时间也没有计算出来，我就没有算了。建议你把复杂的计算结果先vpa一下后在进行后续计算，是不是好一点。

hanxiao · 发表于 2006-3-15 17:27

kevinhcj： 谢谢你，这么认真地对待我的问题，好感动阿。。。谢谢！ 我上次的问题没有表达清楚，请再帮帮忙吧！见下贴

hanxiao · 发表于 2006-3-15 17:29

还是上次的源代码，当执行到b1=int(2*w,t,0,8)这个语句时，前面都可以执行，这句需要很长时间（几个小时吧），然后要么退出，要么就出现ans=NaN。改用数值积分b1=quadl(2*w,t,0,8)，就出现了下面的信息： ??? Error using ==> fcnchk If FUN is a MATLAB object, it must have an feval method.
Error in ==> quadl at 54 f = fcnchk(funfcn); 怎么回事呀？ 源程序： syms z a=[z^2/34^2,z^3/34^3,z^4/34^4] a1=diff(a) a2=diff(a1) m=int((212160-858*z)*a'*a,0,34) kesi=int((212160-858*z)*a',0,34) p=-kesi*2.383*10^-2 r=int(2.1*10^11*(29.962-0.702*z)*a2'*a2,0,34) k2=-int(10^6*a1'*a1,0,34) k=r+k2 syms t e=k-t^2*m d=-0.1*k i=d*d+e*e yi=inv(i)*d*p yr=inv(i)*e*p x=yr(1)+yr(2)+yr(3) y=yi(1)+yi(2)+yi(3) w=x^2+y^2 b1=int(2*w,t,0,8) b2=int(2*t^2*w,t,0,8) double(b1) double(b2) z=34 a2 =[1/578,3/19652*z,3/334084*z^2] x1=a2(1)*yr(1)+a2(2)*yr(2)+a2(3)*yr(3) y1=a2(1)*yi(1)+a2(2)*yi(2)+a2(3)*yi(3) w1=x1^2+y1^2 wj1=12.8*10^11*int(2*w1,t,0,8) wj2=12.8*10^11*int(2*t^2*w1,t,0,8) double(wj1) double(wj2) [em03]

kevinhcj · 发表于 2006-3-15 17:52

你在查查帮助把，quadl的格式如下 q = quadl(fun,a,b) 被积函数fun需要定义成一个函数

kevinhcj · 发表于 2006-3-15 17:54

quadl和int不一样，前者是做数值积分的，后者做符号积分

hanxiao · 发表于 2006-3-15 19:29

kevinhcj: function怎么定义函数阿？你能告诉我么？ 今天我和我同学试了一天也不行，不知道为什么？ 很着急，帮帮忙吧！真的很感谢你。。。

sffei · 发表于 2006-3-15 19:55

function定义函数随便找本书看看就行了

kevinhcj · 发表于 2006-3-15 20:18

% for example & good luck! function y = myfun(x) a=1; b=2; c=3; y = a*x^2+b*x+c;

hanxiao · 发表于 2006-3-16 09:09

我用function定义函数，为什么出来这些？ 我是新手，大家帮帮忙吧！！ function y=myfun(t) ??? function y=myfun(t) | Error: Function definitions are not permitted at the prompt or in scripts.

happy · 发表于 2006-3-16 09:33

老问题了，看帖子http://forum.vibunion.com/forum/viewthread.php?tid=5018吧

hanxiao · 发表于 2006-3-16 12:42

我用function定义函数，为什么或出现下面的情况？ function y=myfun(x) ??? function y=myfun(x) | Error: Function definitions are not permitted at the prompt or in scripts.

kevinhcj · 发表于 2006-3-16 12:56

如果你的主程序不是m函数，而只是m文件，你定义的函数就要另外保存，而且最后文件名和函数名相同。

账号		自动登录	找回密码
密码			我要加入

那位大虾帮帮忙阿？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（hanxiao）那位大虾帮帮忙阿？

回复：（hanxiao）kevinhcj:function怎么定义函数阿？...

回复：（hanxiao）我用function定义函数，为什么出来这...