声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1046|回复: 5

[编程技巧] 关于积分

[复制链接]
发表于 2009-11-4 11:44 | 显示全部楼层 |阅读模式

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

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

x
请教各位
一个有十个变量的复杂表达式,对其中的一个变量x10求积分,最后表达式含有其他九个变量,怎么做,int(f,x10,0,1.5)这样可以吗?会不会求不出,因为我等了很长时间,仍然busy,所以怕是变量太多解不出,各位提提意见
谢谢了

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2009-11-4 19:31 | 显示全部楼层
不好办。第10个变量和前面9个变量的耦合有没有规律。如果有规律自己先凑出点,然后再积。否则真的计算机要累死

评分

1

查看全部评分

 楼主| 发表于 2009-11-5 10:33 | 显示全部楼层

回复 沙发 VibrationMaster 的帖子

为什么这么做有错误?

syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10;
a=(-x1)+x3*cos((4*pi/3)*x10);
b=x2+x3*sin((4*pi/3)*x10);
e=2*atan((1-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))*((-a)/((a^2+b^2)^(1/2)))-((1-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))^2)^(1/2))*(b/(a^2+b^2)^(1/2)))/(((1-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))^2)^(1/2))*((-a)/((a^2+b^2)^(1/2)))-((x4^2+a^2+b^2-x5^2)/(2*x4*(a^2+b^2)^(1/2)))*(b/(a^2+b^2)^(1/2))));
c=a+x4*cos(e);
d=b+x4*sin(e);
m=a+x7*cos(e+x9);
n=b+x7*sin(e+x9);
q=asin((m+x8)/x6);
p=n+x6*cos(q);
f=abs(diff(p,x10));
>> s=int(f,x10,0,1.5)
??? 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(:),[x.s '=(' a.s ')..(' b.s ')']),size(f));
程序及错误如上

我还发现一个小问题,如果写一个简单的a^2+b^2这样的积分,本来能求出,可是运行完上面的程序,出现错误以后,clc clear 然后再运行这个简单的积分,他也会出现上面的错误,这是为什么,对了,这是用7.1解得,7.8运行一天解不出[
这是对一个表达式求导得到的,可以求导,不可以积分?
发表于 2009-11-5 11:29 | 显示全部楼层
 楼主| 发表于 2009-11-7 14:59 | 显示全部楼层

回复 地板 beyondhxf 的帖子

看了看,好像解决不了问题呀?请多指教
 楼主| 发表于 2009-11-13 09:38 | 显示全部楼层
自己顶!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-18 19:33 , Processed in 0.056151 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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