somnus 发表于 2011-5-17 22:11

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

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:38

本帖最后由 321forever 于 2011-5-18 01:40 编辑

回复 1 # somnus 的帖子

lz能不能说下这是什么程序有什么用处,还有给出出错的信息,我刚刚运行了下,是不是 应该把p(1),p(2)...这些变量的括号去掉,p1 p2。。。楼主这个是不是一个周期函数的程序呢


wangyouyi 发表于 2011-5-18 08:26

回复 1 # somnus 的帖子

楼主这个程序好像不是很复杂,只是简单的计算!
具体遇到什么问题也没说清楚!

somnus 发表于 2011-5-18 12:48

回复 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,谢谢了!

somnus 发表于 2011-5-18 12:49

回复 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,谢谢了!

hebut 发表于 2011-5-18 21:19

这个不是一次的方程吧?
如果是一次的方程组就写成*=,其中是未知数组成的向量,这样就容易了。
如果不是一次的方程就比较麻烦。
我曾经算一个八元二次的方程组,结果out of memory了
页: [1]
查看完整版本: 求助!哪位高手帮我看哈这个程序