Question

matlab求解

matlab求解
tg(y/2)=[300-350*cosy+250*x*ctg(y)]/250.
a=350*cosy
b=250+350*siny-250x
求a和b的关系式
需要作图，Y取值0到90度

Answer 1

>> clear

syms a b x y z

A=zeros(1,2);

A=solve('a=350*cos(y)',y);

for i=1:2

y=A(i);

z=tan(y/2)-((300-350*cos(y)+250*x*atan(y))/250);

x=solve('z',x);

disp('得到b和a的关系为：')

b=250+350*sin(y)-250*(x)

end

disp('画出a和b的图像：')

clf

y=1:0.01:pi/2;

a=min(350*cos(y)):0.01:max(350*cos(y));

subs(b);

plot(a,subs(b),'r')

title('b关于a的函数图像');

xlabel('a轴');

ylabel('b轴');

grid on

Warning: Explicit solution could not be found. 

> In solve at 83

得到b和a的关系为：

 

b =

 

[ empty sym ]

 

得到b和a的关系为：

 

b =

 

250 - 350*(1 - a^2/122500)^(1/2)

 

画出a和b的图像：

追问

如果改为
tg(y/2)=[300-250*x*cot(y)]/250.
a=350*cosy
b=250+350*siny-250*x
求a和b的关系式
需要作图，Y取值0到90度呢

怎么和我人工作图的结果不一样？

追答

做法一样的哇，你自己试试呗。。
他们做出来是什么样的图哇？

画出来的结果还是一样的哇。、。、
>> clear
>> syms a b x y z
A=zeros(1,2);
A=solve('a=350*cos(y)',y);
for i=1:2
y=A(i);
z=tan(y/2)-((300-250*x*cot(y))/250);

x=solve('z',x);
disp('得到b和a的关系为：')
b=250+350*sin(y)-250*(x)
end
disp('画出a和b的图像：')
clf
y=1:0.01:pi/2;
a=min(350*cos(y)):0.01:max(350*cos(y));
subs(b);
plot(a,subs(b),'r')
title('b关于a的函数图像');
xlabel('a轴');
ylabel('b轴');
grid on

Warning: Explicit solution could not be found.
> In solve at 83
得到b和a的关系为：

b =

[ empty sym ]

得到b和a的关系为：

b =

250 - 350*(1 - a^2/122500)^(1/2)

和上面的a和b的表达式一样，图像也肯定是一样的哇、、、

追问

Y的取值没有为（0,0.5*pi）不含端点。

追答

大哥，or大姐？我真晕你哇，你又没说清楚要作a和b的图还是x和y的图、、、你们作的是x和y的图，我作的是a和b的图，当然不一样哇、、唉