somnus 发表于 2011-5-19 11:27

跪求!哪位高手帮我看哈这个程序啊!!!

解一个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,谢谢了!
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=;

meiyongyuandeze 发表于 2011-5-19 22:51

感觉不像是一次方程吧,有平方项的!
不知道楼主的错误时是什么,最好能给个报错信息!
页: [1]
查看完整版本: 跪求!哪位高手帮我看哈这个程序啊!!!