Question

如何用matlab中的fsolve或solve函数解这个方程

方程：u.*cos(pi/9)^2*(120+2.125*cot(0.0173*pi)-u.*cos(pi/9)*cos(x))+cos(pi/9)*(2.125*x+u.*sin(pi/9))*(2.125*sin(x)-u.*cos(x)*sin(pi/9))-120*cos(pi/9)*cot(0.0173*pi)*(2.125*cos(x)+u.*sin(pi/9)*sin(x))

参数：u=linspace(64,79.813,1000)
从而得到u与x的一 一对应关系

Answer 1

求解此类问题，可以用循环语句和fsolve函数一起来解。求解过程如下：

u0=linspace(64,79.813,10);

for i=1:length(u0)

u=u0(i);

func=@(x,u)u.*cos(pi/9).^2.*(120+2.125*cot(0.0173*pi)-u*cos(pi/9).*cos(x))+cos(pi/9).*(2.125*x+u*sin(pi/9)).*(2.125*sin(x)-u.*cos(x)*sin(pi/9))-120*cos(pi/9).*cot(0.0173*pi).*(2.125*cos(x)+u*sin(pi/9).*sin(x));

x0=fsolve(@(x)func(x,u),u);

x(i)=x0;

end

u=u0;

[x'  u']

运行结果