MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   MATLAB论坛 (https://www.labfans.com/bbs/forumdisplay.php?f=6)
-   -   [MATLAB数学相关] 一个貌似很简单但是我不会的问题 (https://www.labfans.com/bbs/showthread.php?t=8339)

hly-neu 2009-05-14 10:03

一个貌似很简单但是我不会的问题
 
key1=[10000 10000];
i=1;
syms s1=0 s2=0 D=0;
e=4-sqrt(3);
n=1;
R=1;
c=pi./6;
h=4;
for r=[0:0.05:10]
for b=[0:pi./20:2*pi]
for a=[0:pi./20:2*pi]
y=(1./(2.*pi.*pi.*10.*10.*(((n+1)./(4.*pi.*pi)).*((((e-0.*sin(b)).*(sin(c))+h.*cos(c)).^n*((h-R.*sin(a).*sin(c))))./(((0.*cos(b)-R.*cos(a)).^2+((0.*sin(b)-e).*cos(c)+(h.*sin(c))-R.*sin(a)).^2+((e-0.*sin(b)).*sin(c)+h.*cos(c)).^2).^2))).^((n+3)./2))).*(((((n+1)./(4.*pi.*pi)).*((((e-r.*sin(b)).*(sin(c))+h.*cos(c)).^n*((h-R.*sin(a).*sin(c))))./(((r.*cos(b)-R.*cos(a)).^2+((r.*sin(b)-e).*cos(c)+(h.*sin(c))-R.*sin(a)).^2+((e-r.*sin(b)).*sin(c)+h.*cos(c)).^2).^2))).^((n+3)./2)));
s=(pi./20).*sum(y);
end
k=(pi./20).*sum(s);
end
D=(k.^2).*r.*0.05;
key1(i)=D;
i=i+1;
end


??? Error using ==> syms at 58
Not a valid variable name.
为什么运行程序是总是提示上面的错误啊?
哪位高手帮忙解决一下呗,谢谢!!!


所有时间均为北京时间。现在的时间是 12:46

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.