求助！哪位高手帮我看哈这个程序
function p=mybar(u1,u2,u3,u4,u5,u6)u1=f;
u2=d;
u3=r;
u4=e;
u5=a;
u6=c;
p(1)=u1*((0-u2)*cos(u3)+0*sin(u3))*(0^2+0^2+u4^2-u5^2)-u4*0*((0-u2)^2+0^2+u1^2-u6^2)+u1*((0-u2)*sin(u3)-0*cos(u3))*(0^2+0^2+u4^2-u5^2)+
u4*0*((0-u2)^2+0^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(0*(0-u2)+0^2-u2*0*cot(u3);
p(2)=u1*((1/6-u2)*cos(u3)+1/2*sin(u3))*(1/6^2+1/2^2+u4^2-u5^2)-u4*1/6*((1/6-u2)^2+1/2^2+u1^2-u6^2)+u1*((1/6-u2)*sin(u3)-1/2*cos(u3))*(1/6^2+1/2^2+u4^2-u5^2)+
u4*1/2*((1/2-u2)^2+1/6^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1/6*(1/6-u2)+1/2^2-u2*1/2*cot(u3);
p(3)=u1*((1/4-u2)*cos(u3)+sin(pi/4)*sin(u3))*(1/4^2+sin(pi/4)^2+u4^2-u5^2)-u4*1/4*((1/4-u2)^2+sin(pi/4)^2+u1^2-u6^2)+u1*((1/4-u2)*sin(u3)-sin(pi/4)*cos(u3))*(1/4^2+sin(pi/4)^2+u4^2-u5^2)+
u4*1/4*((1/4-u2)^2+sin(pi/4)^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1/4*(1/4-u2)+sin(pi/4)^2-u2*sin(pi/4)*cot(u3);
p(4)=u1*((1/3-u2)*cos(u3)+sin(pi/3)*sin(u3))*(1/3^2+sin(pi/3)^2+u4^2-u5^2)-u4*1/3*((1/3-u2)^2+sin(pi/3)^2+u1^2-u6^2)+u1*((1/3-u2)*sin(u3)-sin(pi/3)*cos(u3))*(1/3^2+sin(pi/3)^2+u4^2-u5^2)+
u4*1/3*((1/3-u2)^2+sin(pi/3)^2+u1^2-u6^2)-2*u4*f*sin(u3)*(1/3*(1/3-u2)+sin(pi/3)^2-u2*sin(pi/3)*cot(u3);
p(5)=u1*((1/2-u2)*cos(u3)+1*sin(u3))*(1/2^2+1^2+u4^2-u5^2)-u4*1/2*((1/2-u2)^2+1^2+u1^2-u6^2)+u1*((1/2-u2)*sin(u3)-1*cos(u3))*(1/2^2+1^2+u4^2-u5^2)+
u4*1/2*((1/2-u2)^2+1^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1/2*(1/2-u2)+1^2-u2*1*cot(u3);
p(6)=u1*((1-u2)*cos(u3)+0*sin(u3))*(1^2+0^2+u4^2-u5^2)-u4*1*((1-u2)^2+0^2+u1^2-u6^2)+u1*((1-u2)*sin(u3)-0*cos(u3))*(1^2+0^2+u4^2-u5^2)+
u4*1*((1-u2)^2+0^2+u1^2-u6^2)-2*u4*u1*sin(u3)*(1*(1-u2)+0^2-u2*0*cot(u3);
p=;
本帖最后由 321forever 于 2011-5-18 01:40 编辑
回复 1 # somnus 的帖子
lz能不能说下这是什么程序有什么用处，还有给出出错的信息，我刚刚运行了下，是不是 应该把p(1),p(2)...这些变量的括号去掉，p1 p2。。。楼主这个是不是一个周期函数的程序呢
回复 1 # somnus 的帖子
楼主这个程序好像不是很复杂，只是简单的计算！
具体遇到什么问题也没说清楚！ 回复 2 # 321forever 的帖子
就是解一个6元1次方程组，
U=f[(x-d)cosr+ysinr](x^2+y^2+e^2-a^2)-ex[(x-d)^2+y^2+f^2-c^2]
V=f[(x-d)sinr-ycosr](x^2+y^2+e^2-a^2)+ex[(x-d)^2+y^2+f^2-c^2]
W=2efsinr
U+V=W
取正弦函数sin(pix)上的6个点（0,0）；（1/6,1/2）； ; ; (1/2,1);(1 ,0) 代入上式得到6元1次方程组，求解未知数a ,c , d, e, f, r,谢谢了！
回复 3 # wangyouyi 的帖子
就是解一个6元1次方程组，
U=f[(x-d)cosr+ysinr](x^2+y^2+e^2-a^2)-ex[(x-d)^2+y^2+f^2-c^2]
V=f[(x-d)sinr-ycosr](x^2+y^2+e^2-a^2)+ex[(x-d)^2+y^2+f^2-c^2]
W=2efsinr
U+V=W
取正弦函数sin(pix)上的6个点（0,0）；（1/6,1/2）； ; ; (1/2,1);(1 ,0) 代入上式得到6元1次方程组，求解未知数a ,c , d, e, f, r,谢谢了！
这个不是一次的方程吧？
如果是一次的方程组就写成*=,其中是未知数组成的向量，这样就容易了。
如果不是一次的方程就比较麻烦。
我曾经算一个八元二次的方程组，结果out of memory了
页:
[1]