Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

hanxiao 发表于 2006-3-15 15:40

happy大侠，请进!上次没说清，请再帮帮忙阿，急死我了！

大侠，上次说错了！事情的经过是这样子的：还是上次的源代码，当执行到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= 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 = 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-15 19:31

happy大侠： 你跑拿去了？！

happy 发表于 2006-3-16 09:24

回复：（hanxiao）happy大侠，请进!上次没说清，请再帮...

好奇怪，今天重新运行了一下居然没问题了 结果是 double(wj1)=4.3267e-002 double(wj2)=1.8272e+000

hanxiao 发表于 2006-3-16 11:39

大侠: 这究竟是怎么回事呢?你算了很长时间么？ 我和我同学昨天研究一天也没研究明白，觉得这个程序中了邪了！！ 我曾经用它算过a=,a为1*2阶的，很好用。而且你的3阶 和我2阶的结果很接近，这也是理论上应该得到的结果。为什么我就算不了？ 你的matlab是什么版本？

kevinhcj 发表于 2006-3-16 12:53

你的式子太复杂，中间变量也很大，我算的时候也是自动退出，这是内存的问题，如果内存足够大，也没有关系。当然，如果内存一般，运气好的时候也可能算出来，呵呵。

hanxiao 发表于 2006-3-16 21:39

kevinhcj: 真是太感谢你了！处处支持我！。。。谢谢阿！！

hanxiao 发表于 2006-3-16 21:39

同时也要感谢happy大侠！

happy 发表于 2006-3-17 07:24

回复：（kevinhcj）你的式子太复杂，中间变量也很大，...

<DIV class=quote>以下是引用kevinhcj在2006-3-16 12:53:57的发言： 你的式子太复杂，中间变量也很大，我算的时候也是自动退出，这是内存的问题，如果内存足够大，也没有关系。当然，如果内存一般，运气好的时候也可能算出来，呵呵。</DIV>
 很有道理，我第一次运行的时候，计算机已经运行了很久，估计内存被吃掉了很多 第二次运行的时候计算机刚刚重新启动

hanxiao 发表于 2006-3-17 13:58

如果这样我就对自己的程序比较有信心了，哪天找个快点的机器再算一下！ 谢谢大家的支持！！

hanxiao 发表于 2006-3-19 10:05

happy大侠： 我还是觉得很奇怪，我昨天特意找了一个很不错的机器，一启动就运行了程序，还是不行。 出现了下面的错误信息： Warning: Explicit integral could not be found. > In sym.int at 58
??? Error using ==> reshape To RESHAPE the number of elements must not change.
Error in ==> sym.maple at 94 result = reshape(result,size(varargin{3}));
Error in ==> sym.int at 51 r = reshape(maple('map','int',f(:),),size(f)); 这些都是什么意思呀？reshape命令我昨天也没研究明白，是元素维数不对么？ 还有大侠，你还能告诉我double(b1)和double(b2)的值么？麻烦你了。。。 谢谢。。。

happy 发表于 2006-3-19 10:16

回复：（hanxiao）happy大侠：我还是觉得很奇怪，我...

<DIV class=quote>以下是引用hanxiao在2006-3-19 10:05:41的发言： 
happy大侠： 我还是觉得很奇怪，我昨天特意找了一个很不错的机器，一启动就运行了程序，还是不行。 出现了下面的错误信息： Warning: Explicit integral could not be found. > In sym.int at 58
??? Error using ==> reshape To RESHAPE the number of elements must not change.
Error in ==> sym.maple at 94 result = reshape(result,size(varargin{3}));
Error in ==> sym.int at 51 r = reshape(maple('map','int',f(:),),size(f)); 这些都是什么意思呀？reshape命令我昨天也没研究明白，是元素维数不对么？ 还有大侠，你还能告诉我double(b1)和double(b2)的值么？麻烦你了。。。 谢谢。。。</DIV>
 reshape的问题一直搞不太清楚，不好意思 另外double(b1)和double(b2)上面的回帖中我已经给出了

hanxiao 发表于 2006-3-19 21:45

不！上次给的是double(wj1)和double(wj2). 其实我要的结果就是上面的四个数：double(wj1)double(wj2) 还有double(b1)和double(b2). 拜托大侠了！！

hanxiao 发表于 2006-3-20 14:26

ding

happy 发表于 2006-3-20 15:17

回复：（hanxiao）不！上次给的是double(wj1)和double...

<DIV class=quote>以下是引用hanxiao在2006-3-19 21:45:19的发言： 不！上次给的是double(wj1)和double(wj2). 其实我要的结果就是上面的四个数：double(wj1)double(wj2) 还有double(b1)和double(b2). 拜托大侠了！！ </DIV>
4.2283e-007 1.0471e-005

hanxiao 发表于 2006-3-21 09:08

happy大侠： 你可真是大好人啊！！。太感谢你了！！ 不过，昨天我已经用数值积分得到了其他相关的数据。你给我的数据和数值积分得到的数据正好吻合！ 不过出于严谨的科研态度，我还是不知道为什么你可以算，而我却算不了？我机器内存是512的，版本matlab7。其实我昨天用原来的程序又算了一遍，还是自动退出。我同学1g的内存，matlab6.5，也白搭！ 大侠，你太好了，小弟谢谢大哥了！

页: [1] 2

声振论坛's Archiver

happy大侠，请进!上次没说清，请再帮帮忙阿，急死我了！

回复：（hanxiao）happy大侠，请进!上次没说清，请再帮...

回复：（kevinhcj）你的式子太复杂，中间变量也很大，...

回复：（hanxiao）happy大侠： 我还是觉得很奇怪，我...

回复：（hanxiao）不！上次给的是double(wj1)和double...

回复：（hanxiao）happy大侠：我还是觉得很奇怪，我...