matlab解方程问题！！！

Question

p=linspace(6.22*10^5,12.4*10^6,100);
H=2.45*10^8;
Ax=p./H;
xx=zeros(1,100);
for n=1:100
solve('Ax(n)=0.0487*x^4-0.6591*x^3+3.3665*x^2-7.7014*x+6.6682');
x=x(1);
xx(n)=x;
end
上面的程序运行不出来。就是想运行出满足100个等式的x的值，不知道怎么解不出来。求高手帮忙解答！

Answer 1

你能保证每个方程只有一个解吗？我觉得不可能，因为都是四次方程，所以每个都有四个解。solve会返回所有的解，还包括非实数的复数解，而且你用solve的方法不对。

我给你改了下一个程序，返回的结果xx每一列含有每个方程的四个解。你试试能用不，我觉得应该有很多复数解。
clear all

p=linspace(6.22*10^5,12.4*10^6,100);

H=2.45*10^8;

Ax=p./H;

xx=zeros(4,100);

for n=1:100

syms x

F=0.0487*x^4-0.6591*x^3+3.3665*x^2-7.7014*x+6.6682-Ax(n);

xx(:,n)=solve(F);

end

追问

我知道有四组解，所以在返回的时候我只取了第一个值。
那请问有没有什么方法能解决这个问题？

追答

见改好的程序，你可以只取xx的任意一行。

Answer 2

程序改成如下：
syms x a
p=linspace(6.22*10^5,12.4*10^6,100);
H=2.45*10^8;
Ax=p./H;
X=zeros(4,100);

for n=1:100

a=solve(Ax(n)==0.0487*x^4-0.6591*x^3+3.3665*x^2-7.7014*x+6.6682,'x');
X(:,n)=vpa(a,6);

end
X%解x的值。

最后结果X
有100列。第n列对应第n个等式的解。每个等式能解出4个x值。每列有两个实数解两个虚数解