matlab解复杂方程fsolve（MATLAB用fsolve解方程）

本文目录

• MATLAB用fsolve解方程
• matlab 用fsolve解方程组的问题
• matlab中fsolve解二元方程问题
• MATLAB中fsolve如何求解n多个代数方程组
• 请问怎么用matlab的fsolve函数怎么求解这个方程呢

MATLAB用fsolve解方程

fsolve要求所有内容都是double型，你这里r1，r2，r3
是sym变量，所以出错，改为下面：
f = @(n)sqrt(n2^2-n^2)*k0*H-atan(sqrt(n^2-n1^2)*k0/(sqrt(n2^2-n^2)*k0))-atan(sqrt(n^2-n3^2)*k0/(sqrt(n2^2-n^2)*k0));
》》 fsolve(f,1.6)

matlab 用fsolve解方程组的问题

1. fsolve函数是用于求解非线性方程的，但是在你的程序里没看到非线性方程

2. 在你所给的求根初值下，你所定义的变量 fff 的值为124，而函数diff(X,n)表示求向量X的前后两项之差并重复n次，例如diff(X,2)就相当于 diff(diff(X)),在你的程序中有如下语句：

   f(1)=diff(fff,b1);

   式子右侧的fff,b1均为一个数，则diff(fff,b1)的计算结果为 ,将一个空的结果赋值给f(1)从而导致如下错误(赋值符号左右维数不匹配):

   In an assignment  A(I) = B, the number of elements in B and I must be the same.

        你先把上面的问题解决吧，如果还不可以，麻烦把你的方程贴出来 :)

matlab中fsolve解二元方程问题

function q=fxy(p)
x=p(1);y=p(2);
q(1)=2*x-y-exp(-x);
q(2)=-x+2*y-exp(-y);
end
》》 x=fsolve(@fxy,’,optimset(’Display’,’off’))
x =
0.5671
0.5671

MATLAB中fsolve如何求解n多个代数方程组

matlab中解方程组还是很方便的，例如，对于代数方程组Ax=b(A为系数矩阵，非奇异)的求解，MATLAB中有两种方法：
(1)x=inv(A)*b — 采用求逆运算解方程组；
(2)x=A\B — 采用左除运算解方程组
PS：使用左除的运算效率要比求逆矩阵的效率高很多~
例:
x1+2x2=8
2x1+3x2=13
》》A=;
》》x=inv(A)*b
x =
2.00
3.00
》》x=A\B
x =
2.00
3.00；
即二元一次方程组的解x1和x2分别是2和3。
对于同学问到的用matlab解多次的方程组，有符号解法，方法是：先解出符号解,然后用vpa(F,n)求出n位有效数字的数值解.具体步骤如下：
第一步:定义变量syms x y z ...;
第二步:求解=solve(’eqn1’,’eqn2’,...,’eqnN’,’var1’,’var2’,...’varN’);
第三步:求出n位有效数字的数值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);...。
如：解二（多）元二（高）次方程组：
x^2+3*y+1=0
y^2+4*x+1=0
解法如下：
》》syms x y;
》》=solve(’x^2+3*y+1=0’,’y^2+4*x+1=0’);
》》x=vpa(x,4);
》》y=vpa(y,4);
结果是：
x =
1.635+3.029*i
1.635-3.029*i
-.283
-2.987
y =
1.834-3.301*i
1.834+3.301*i
-.3600
-3.307。
二元二次方程组，共4个实数根；
解答如下：
基本方法是：solve(s1,s2,…,sn,v1,v2,…,vn)，即求表达式s1,s2,…,sn组成的方程组，求解变量分别v1,v2,…,vn。
具体例子如下：
x^2 + x*y + y = 3
x^2 - 4*x + 3 = 0
解法：
》》 = solve(’x^2 + x*y + y = 3’,’x^2 - 4*x + 3 = 0’)
运行结果为
x =
1 3
y =
1 -3/2
即x等于1和3；y等于1和-1.5
或
》》 = solve(’x^2 + x*y + y = 3’,’x^2 - 4*x + 3= 0’,’x’,’y’)
x =
1 3
y =
1 -3/2
结果一样，二元二方程都是4个实根。
通过这三个例子可以看出，用matlab解各类方程组都是可以的，方法也有多种，只是用到解方程组的函数，注意正确书写参数就可以了，非常方便。
2、变参数非线性方程组的求解
对于求解非线性方程组一般用fsolve命令就可以了，但是对于方程组中某一系数是变化的，该怎么求呢？
%定义方程组如下，其中k为变量
function F = myfun(x,k)
H=0.32;
Pc0=0.23;W=0.18;
F=[Pc0+H*(1+1.5*(x(1)/W-1)-0.5*(x(1)/W-1)^3)-x(2);
x(1)-k*sqrt(x(2))];
%求解过程
H=0.32;
Pc0=0.23;W=0.18;
x0 = ; % 取初值
options = optimset(’Display’,’off’);
k=0:0.01:1; % 变量取值范围
for i=1:1:length(k)
kk=k(i);
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非线性方程组
x1(i)=x(1);
x2(i)=x(2);
end
plot(k,x1,’-b’,k,x2,’-r’);
xlabel(’k’)
legend(’x1’,’x2’)

请问怎么用matlab的fsolve函数怎么求解这个方程呢

满满的字，看着就怕。。。
（1）先将你的算式中的r用x(1)代替，a用x(2)代替。然后在前面加上f=@(x)并去掉后面的=0
（2）根据你的几何意义，找出一个满足你要求的可能点，比如r=1,a=0（或你的初始条件）（在全域可能具有无穷多解。给出初始点后，fslove将在这个点附近找满足方程的最近的一个解）
（3）运行fsolve(f,格式的都可以
全部程序如下
f=@(x)(2.*x(1) - (8228485965250453.*cos((pi.*(x(2) + 72))./180))./2251799813685248)./(4.*(x(1).^2 - (8228485965250453.*cos((pi.*(x(2) + 72))./180).*x(1))./2251799813685248 + 7517095976830443./2251799813685248).^(1./2)) - (2.*x(1) - (8228485965250453.*cos((pi.*(x(2) - 72))./180))./2251799813685248)./(4.*(x(1).^2 - (8228485965250453.*cos((pi.*(x(2) - 72))./180).*x(1))./2251799813685248 + 7517095976830443./2251799813685248).^(1./2)) - (2.*sin((pi.*(x(2) - 6))./180) - 2.*x(1) + (8228485965250453.*sin((pi.*(x(2) + 6))./180))./2251799813685248)./(2.*(x(1).^2 - (8228485965250453.*x(1).*sin((pi.*(x(2) + 6))./180))./2251799813685248 - 2.*x(1).*sin((pi.*(x(2) - 6))./180) + 17817569595097215./2251799813685248).^(1./2)) - (2.*x(1) - (8228485965250453.*cos((pi.*(x(2) - 144))./180))./2251799813685248)./(4.*(x(1).^2 - (8228485965250453.*cos((pi.*(x(2) - 144))./180).*x(1))./2251799813685248 + 7517095976830443./2251799813685248).^(1./2)) - (2.*x(1) - (8228485965250453.*cos((pi.*(x(2) + 144))./180))./2251799813685248)./(4.*(x(1).^2 - (8228485965250453.*cos((pi.*(x(2) + 144))./180).*x(1))./2251799813685248 + 7517095976830443./2251799813685248).^(1./2)) - (2.*x(1) - (8228485965250453.*cos((pi.*x(2))./180))./2251799813685248)./(4.*(x(1).^2 - (8228485965250453.*cos((pi.*x(2))./180).*x(1))./2251799813685248 + 7517095976830443./2251799813685248).^(1./2));
x0=fsolve(f,);x0为解