声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2130|回复: 8

[编程技巧] 请帮我看看这个数值积分的程序

[复制链接]
发表于 2007-6-3 14:28 | 显示全部楼层 |阅读模式

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

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

x
我运行这个程序:
syms x y t
E= zeros(9);  
for i = 1:9,  
for j = 1:9,  
E(i,j) =int(2*sin(t)/(10000*((i-5)/200)-273-102*cos(2*t)+10000*((j-5)/200)^2-400*(j-5)/200*sin(t))^(1/2),t,0,2*pi);  
end  
end  
disp(E)

运行了很久(大概十来个小时)才得到这个结果
??? Error using ==> sym/maple
System error, ran out of memory

Error in ==> D:\toolbox\symbolic\@sym\int.m
On line 51  ==>    r = maple('map','int',f,[x.s '=' a.s '..' b.s]);

请帮我修改一下,谢谢!
回复
分享到:

使用道具 举报

发表于 2007-6-3 15:46 | 显示全部楼层
楼主为什么不先积分后循环
发表于 2007-6-3 17:17 | 显示全部楼层
楼主为什么不用数值积分.难道非用符号积分吧?
发表于 2007-6-3 18:50 | 显示全部楼层
好像不是很复杂的积分, 可以算出定积分来的。
 楼主| 发表于 2007-6-3 22:21 | 显示全部楼层
E(i,j) =int(2*sin(t)/(10000*((i-5)/200)-273-102*cos(2*t)+10000*((j-5)/200)^2-400*(j-5)/200*sin(t))^(1/2),t,0,2*pi);


这个式子就是通过int(0.02*sin(t)/sqrt(x-0.02-0.02*cos(t)^2+(y-0.02*sin(t))^2+0.0025),t,0,2*pi)得到的结果。
 楼主| 发表于 2007-6-3 22:47 | 显示全部楼层
我用quad函数,举个例子
F = inline('0.02*sin(x)/sqrt(1-0.02-0.02*cos(x).^2+(1-0.02*sin(x).^2+0.0025))');
Q = quad(F,0,2);
disp(Q)

为什么会显示一下错误啊?
???  Index exceeds matrix dimensions.

Error in ==> D:\toolbox\matlab\funfun\quad.m
On line 70  ==> if ~isfinite(y(7))
发表于 2007-6-3 23:15 | 显示全部楼层
原帖由 lenapu 于 2007-6-3 22:47 发表
我用quad函数,举个例子
F = inline('0.02*sin(x)/sqrt(1-0.02-0.02*cos(x).^2+(1-0.02*sin(x).^2+0.0025))');
Q = quad(F,0,2);
disp(Q)

为什么会显示一下错误啊?
???  Index exceeds matrix dimensi ...


用 ./ 非 /
 楼主| 发表于 2007-6-4 09:35 | 显示全部楼层
我想让quad参与循环,比如:
A=zeros(10);
for k=1:10
F = inline('0.02*sin(x)+k');
A(k)= quad(F,0,2);
end
disp(A)

可是却又显示错误,是不是积分函数里面除了积分变量以外还含有其他的变量的时候,就不能应用quad函数?
??? Error using ==> inline/feval
Not enough inputs to inline function.

Error in ==> D:\toolbox\matlab\funfun\quad.m
On line 62  ==> y = feval(f, x, varargin{:});
 楼主| 发表于 2007-6-4 15:32 | 显示全部楼层
看了‘“关于求解变参数微分方程”回复专贴’
明白了
呵呵!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 13:39 , Processed in 0.056575 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表