帮忙用matlab作个小程序

来源：百度知道 编辑：UC知道 时间：2024/05/30 08:37:41
P1=2T1cosB+Tcos(B+A/2)+Tcos(B-A/2)+Nsin(B+A/2)-Nsin(B-A/2)
P2=2T1sinB+Tsin(B+A/2)+Tsin(B-A/2)-Ncos(B+A/2)+cos(B-A/2)
P*P=P1*P1+P2*P2
作出当T=0.577*167000*0.002/sinA*0.04
T1=0.577*167000*(0.01-0.002/tanA)0.04
N=167000*0.002/sina*0.04
时，P1与A，B
P2与A，B
P与A，B三个关系曲线图
菜鸟急需！多谢！
A的范围是0～90
B的范围是0～90
PP1P2的范围是0～5000

syms A B
T=0.577*167000*0.002/sin(A)*0.04;
T1=0.577*167000*(0.01-0.002/tan(A))*0.04;
N=167000*0.002/sin(A)*0.04;
P1=2*T1*cos(B)+T*cos(B+A/2)+T*cos(B-A/2)+N*sin(B+A/2)-N*sin(B-A/2)
P2=2*T1*sin(B)+T*sin(B+A/2)+T*sin(B-A/2)-N*cos(B+A/2)+cos(B-A/2)
P=sqrt(P1*P1+P2*P2)
figure(1)
ezsurf(P1)
figure(2)
ezsurf(P2)
figure(3)
ezsurf(P)

结果：
P1 =

(96359/1250-96359/6250/tan(A))*cos(B)+96359/12500/sin(A)*cos(B+1/2*A)+96359/12500/sin(A)*cos(-B+1/2*A)+334/25/sin(A)*sin(B+1/2*A)+334/25/sin(A)*sin(-B+1/2*A)

P2 =

(96359/1250-96359/6250/tan(A))*sin(B)+96359/12500/sin(A)*sin(B+1/2*A)-96359/12500/sin(A)*sin(-B+1/2*A)-334/25/sin(A)*cos(B+1/2*A)+cos(-B+1/2*A)

P =

(((96359/1250-96359/6250/tan(A))*cos(B)+96359/12500/sin(A)*cos(B+1/2*A)+96359/12500/sin(A)*cos(-B+1/2*A)+334/25/sin(A)*sin(B+1/2*A)+334/25/sin(A)*sin(-B+1/2*A))^2+((96359/1250-96359/6250/tan(A))*sin(B)+96359/125