Question

matlab解非线性方程组

接下来是：
c*x*(h^2 + x^2 + y^2 + h^2*cotd(30)^2 - 2*h*x*cotd(30))^(1/2)*(h^2 + x^2 + y^2 + h^2*cotd(45)^2 - 2*h*x*cotd(45))^(1/2)*(h^2 + x^2 + y^2 + h^2*cotd(60)^2 - 2*h*x*cotd(60))^(1/2) - c*h*cotd(45)*(h^2 + x^2 + y^2 + h^2*cotd(30)^2 - 2*h*x*cotd(30))^(1/2)*(h^2 + x^2 + y^2 + h^2*cotd(60)^2 - 2*h*x*cotd(60)) + c*h*cotd(60)*(h^2 + x^2 + y^2 + h^2*cotd(30)^2 - 2*h*x*cotd(30))^(1/2)*(h^2 + x^2 + y^2 + h^2*cotd(45)^2 - 2*h*x*cotd(45))^(1/2)*(h^2 + x^2 + y^2 + h^2*cotd(60)^2 - 2*h*x*cotd(60))^(1/2);
S=solve(f1,f2,'x','y')
方程2很蛋疼啊

Answer 1

可以使用solve函数。如下参考：

例如：

x^2+y^3=10

x^3-y^2=1

其中x,y为方程组的未知量，在Matlab的命名窗口中输入：

y =

- 37^(1/2)/2 - 1/2

37^(1/2)/2 - 1/2

37^(1/2)/2 - 1/2

- 37^(1/2)/2 - 1/2

syms x y

[x y]=solve('x^2+y=10','x^2-y^2=1','x','y')输出计算结果为：x =

(37^(1/2)/2 + 21/2)^(1/2)

(21/2 - 37^(1/2)/2)^(1/2)

-(21/2 - 1/2*37^(1/2))^(1/2)

-(1/2*37^(1/2) + 21/2)^(1/2)

扩展资料：

注意事项

求解是一个基本的内置函数，用于方程的符号解，它返回一个符号变量矩阵(m×nm×n sym)。当没有可用的符号解时，将抛出一个警告并输出一个数值解。

基本形式为:solve(eqn, var, Name, Val);

% eqn是符号表达式/符号变量/符号表达式的函数句柄，var是未知数;Name是附加的需求，Val是它的值。解一维方程。对于多项式，solve返回它的所有值。

函数1等于@(x)x^3-20 x^2-25 x + 500;创建一个函数句柄。句柄中的变量不是符号变量，不需要定义。

syms x exp1;    % 定义符号变量 x, exp1；

exp1 = x^3 - 20*x^2 - 25*x + 500;    % 符号表达式，包含符号变量. 符号变量必须先有上一行定义。

solve(exp1 == 0, x)    % 命令行输入a，传入一个包含符号表达式的等式，x为所要求的变量

solve(exp1, x)    % 命令行输入b，传入一个符号表达式，函数默认求其零点

solve(func1(x), x)    % 命令行输入c，传入参数func1(x)等价于传入了符号表达式，和输入b完全一样

solve(func1(x) == 0, x)    % 命令行输入d，这句话和a完全一样

solve(func1, x)    % 命令行输入e，传入参数func1，这是一个函数句柄，函数默认求其零

ans =    % 命令行输出，三个解，为3*1的符号向量。对以上五种输入输出都完全一样

-5

5

20

对于不可符号求解的函数零点/方程解，solve抛出警告并返回一个数值解：

exp1 = atan(x) - x - 1;    % 不可符号求零点的表达式

solve(exp1 == 0, x)    % 命令行输入

% 命令行输出：

警告: Cannot solve symbolically. Returning a numeric approximation instead.

ans =-2.132267725272885131625420696936

Answer 2

matlab中有专门的函数来解决方程组的，我给你举一个例子好了，你一看就会了。这个我帮别人写过的一道题，分享给你
(a-x)^2+(b-y)^2=e^2
(C-x)^2+(D-y)^2=v^2
已知a,b,c,d,e,v 值
求解 X,Y 请问用 matlab 如何写，就是求2个园的交点问题。
仿真程序为：
global a b c d e v;
>> a=1;b=0;c=-1;d=0;e=1.5;v=1.5;%%%%%%%设定你这几个未知数的值
>> syms x y;%%%%%%x,y是变量
>> [x,y]=solve('x^2+y^2-2*a*x-2*b*y=e^2-a^2-b^2','x^2+y^2-2*c*x-2*d*y=v^2-c^2-d^2');%%%%我把平方展开了
>> x=vpa(x,4);y=vpa(y,4);%%%%%%%%%%%取4位有效数字
%%%%%%%%%%
格式就是用solve(方程1，方程2，…求解变量1,变量2，…）；
希望对您有帮助

追问

这是原始公式，其中已知量的值：v=800、c=300000、f=120 000 (就是120MHz，带入公式时单位换算成KHz)、h=10。谢谢啊。

Answer 3

用solve函数，格式solve(方程1，方程2，…求解变量1,变量2，…）；
其中方程“代数式=0”只写成代数式，不写变量的话，默认有findsym自动确定。
g=solve(eq1,eq2,…,val1,val2…）所得g为构架数组，显示每个变量的结果用指令g.val1, g.val2
例子：
syms x y
g=solve(x^2+y-8,x-y^2+y-10,x,y);
g.x
g.y