如何对代码中大量内容进行修改
前几天发帖问过一次,现在遇到新的问题,下面是一个积分运算clear
syms E X Nrd r t Rd Rr P0 r1 r2 N1 N2 N3 Sec G d;
d=0.7;
p=8/pi;
Rr=1;
%Rd=10^(1/4)*Rr;
Rd=2*Rr;
X=(Rd^2+r^2-Rr^2)/(2*r);
Sec=(Rd^2*acos(X/Rd)+Rr^2*acos((r-X)/Rr)-r*(Rd^2-X^2)^0.5);
P1=Sec/pi/Rr^2;
r1=d*cos(G)+((d*cos(G))^2-d^2+Rd^2)^0.5;
r2=Rd+Rr;
N1=int(p*r*P1,'r',r1,r2) %N1是G的一个大函数
在这里matlab算出来 N1 是一个G为参数的大函数,大约有几万个字符,是有无数个包含cos(G)的复杂表达式,
现在需要对G积分区间是【0 ,2*pi】,
我想用quadl积分,下面这么写,错误提示很多地方需要改为点乘,可是里面的内容是在太多,如何修改呢?谢谢!
myfun = @(G) (N1大函数的字符串);
Q = quadl(myfun,0,2*pi )
[ 本帖最后由 eight 于 2008-1-22 20:02 编辑 ] 查找,替换?
回复 楼主 的帖子
复杂表达式,建议还是一开始就数值积分回复 3楼 的帖子
多谢指点,在这里找到了方法http://biosanger.spaces.live.com/blog/cns!aa5879a9e6e2616f!276.entry
页:
[1]