求matlab解三元方程组,最好能把命令行列出来

作者&投稿:佟柳 (若有异议请与网页底部的电邮联系)
求解一个方程组,能用matlab并把命令行展示出来最好~


用vpaslove()函数可以直接求得t,x,y,z值。
syms x y z t
[t x y z]=vpasolve(y*cos(z)+x==350,y*cos(t*x)==-y*sin(z),1/tan(z)==-t*y*sin(t*x),y^4-(y^2-112.5)^2==y^2*t^2*(y^2-112.5)^2)

方程无解。可以很简单地画出两个方程的曲线来:
ezplot('-0.050+0.015+0.070-2.*(d.*0.0215).^0.5-2.*(p.*0.0504).^0.5',[-0.02 0.02])
hold on
ezplot('-0.015+d+p',[-0.02 0.02])

由图可见,两条线无交点,也就意味着方程组无解。

或者换另一种方式来分析:把方程1代入方程2:
syms d p
eq1=-0.015+d+p;
eq2=-0.050+0.015+0.070-2.*(d.*0.0215).^0.5-2.*(p.*0.0504).^0.5;
d=solve(eq1,d)
eq=subs(eq2,d,d)得到方程:
>> eq=subs(eq2,d,d)
eq =
7/200 - 2*(129/400000 - (43*p)/2000)^(1/2) - 2*((63*p)/1250)^(1/2)方程中两个包含p的项都要进行开平方,应该大于0,即:
(63*p)/1250
>=
0
129/400000
-
(43*p)/2000
>=
0
联立解得
0<=p<=0.015,画出该方程的曲线如下:

不可能与0轴相交,也即方程无解。

求matlab解此三元方程组,可以考虑用fsove()函数来求解。

1、首先自定义函数fun(x)

fun=@(x)[11*cos(x(2))+x(1)-350,cos(x(3).*x(1))+sin(x(2)),1/tan(x(2))+x(3)*11.*sin(x(3).*x(1))]; %x=x(1),y=x(2),z=x(3)

2、然后用fsove()函数求解

x0=[300 2 -0.01], [x,h]=fsolve(fun,x0)

3、求解结果

x≈350.01;y≈1.5713;z≈0.044877

4、求解界面




MATLAB 解三元二次方程
当然如果p1!=P2!=P3也可以解,并且也只有一个解,不过就不是一个圆,没有这么好理解。如果还有不明白的可以联系我。解这样的问题还用动matlab?补充说明:我试了一下,确实动手解了一下方程组,其实还是有点问题的,主要是得到的解误差太大。比如1400mm的标准解,实际得到的答案可能是1380-1420之间...

matlab 解三元二次方程组,求大神级解答
建议用数值方法解,别用符号方法。试试不动点迭代法或者牛顿迭代法。二分法可能也行。

用MATLAB求解三元非线性方程
syms x y z;[x,y,z]=solve('x^2+2*x+1=0','x+3*z=4','y*z=-1')x = -1 y = -3\/5 z = 5\/3

哪位热心人能帮我用matlab解三元 三角函数方程组
ep1 = 'tan(st1)=(1-a^2)*sin(st2)\/((1+a^2)*cos(st2)-2*a)';ep2 = 'tan(st2-0.4887)=(1-a^2)*sin(st2+0.1769)\/(2*a+(a^2+1)*cos(st2+0.1769))';ep3 = 'tan(st2+0.4887)=(1-a^2)*sin(st2+1.7826)\/(2*a+(a^2+1)*cos(st2+1.7826))';[st1,st2,...

matlab解三元二次方程 当x不等于0时可以解除结果,但x=0时却解不出,这...
按下面修改:b1=-pi\/6;b2=pi\/2;b3=7*pi\/6;a1=0.790581791275871;a2=1.595624508069142;a3=0.790581791275871;x=0;y=400;z=-940;syms Rr l1 l2 ;eq1=(x-cos(b1)*Rr-l1*cos(b1)*sin(a1))^2+(y-sin(b1)*Rr-l1*sin(b1)*sin(a1))^2+(z+l1*cos(a1))^2-l2^2;eq2=(x-...

matlab解三元二次方程 为什么最后的解是矩阵样式呢? 哪里出错了呢_百度...
syms a0y a1x a1y a0x;eq1=((a1x+1)-a0x)^2+((a1y-0.5)-a0y)^2-(a1x-a0x)^2-(a1y-a0y)^2;eq2=((2-a1y)-a0x)^2+((a1x+3)-a0y)^2-(a1x-a0x)^2-(a1y-a0y)^2;eq3=(((sqrt(2)\/2)*a1x-(sqrt(2)\/2)*a1y+3)-a0x)^2+(((sqrt(2)\/2)*a1x+(sqrt(...

MatLab求解三元二次方程组
(x-x1)²+(y-y1)²+(z-z1)²=d1²(x-x2)²+(y-y2)²+(z-z2)²=d2²(x-x3)²+(y-y3)²+(z-z3)²=d3²(x-x4)²+(y-y4)²+(z-z4)²=d4²计算出位置点坐标 行列式计算也...

MATLAB 求解三元多次方程组
syms x y z;equ1=x\/2==y\/218000+(y\/2840)^(1\/0.07)equ2=x*y==7.9\/100000 equ3=x==0.89*(2*z)^(-0.07)+2.67*(2*z)^(-0.96)equs=[equ1,equ2,equ3];vars=[x,y,z];[solve_x,solve_y,solve_z]=solve(equs,vars)建议尝试用优化工具箱里的优化算法求解问题。

matlab问题,高分求matlab求三元方程的极值。
如果用xz-x+>0与xz-y-1>0约束一下,第一个约束不清晰 xz-x+>0===》xz-x+1>0 function aa y=@(x)-(x(1)-1).*(1+x(2)-x(1).*x(3)).\/((1+x(2)).*(x(1)-1)-(x(1).*x(3)).^2);[X,FVAL,EXITFLAG] =fmincon(y,[1.5 1 0.5],[],[],[],[],[1 0...

高分求Matlab高手啊: 如何求''未知''三元函数最小值。
像这种问题Matlab库函数直接可以求解。编写三个m文件,一个主程序,两个function 主程序main.m p0=[0 0 0];A=[];b=[];Aeq=[];beq=[];vlb=[-2;-2;-2];vub=[2;2;2];[x,fval]=fmincon('fun_u',p0,A,b,Aeq,beq,vlb,vub,'confun_u')fun_u.m子函数即目标函数u function y=...

芜湖县15114582219: 求高手用matlab解一个三元的方程组600+Xcos10+Ysin40 - Z=3600(Zsin65+Ysin40+Xcos10)^2+(Zcos65+Ycos40+Xsin10)^2=5200^2Xsin60+Ysin80+Z - 600=... -
朱帖福尔:[答案] 我算的你这个解都是无穷大.具体过程如下:第一步:建立M文件函数myf.mfunction f=myf(x)f(1)=600+x(1)*cos(10)+x(2)*sin(40)-x(3);f(2)=(x(3)*sin(65)+x(2)*sin(40)+x(1)*cos(10))^2+(x(3)*cos(65)+x(2)*cos(40)+x(1...

芜湖县15114582219: MATLAB如何求解三元一次方程组?比如方程组如下:x+y+z=1x+2y+z=2z+y - z=0如何使用MATLAB求解方程组的根? -
朱帖福尔:[答案] s=solve('x+y+z=1','x+2*y+z=2','x+y-z=0','x,y,z'); result=[s.x,s.y,s.z]

芜湖县15114582219: 求高手!用牛顿迭代法解三元三次方程的MATLAB程序! -
朱帖福尔: 用fsolve可解出来: 先构造函数: function output=solveproblem(X) c=X(1); m=X(2); y=X(3); output(1)=(1-c)*(1-y)*(1-m)*10.61+y*(1-c)*(1-m)*30.56+c*(1-y)*(1-m)*60.67+m*(1-c)*(1-y)*0.95+y*m*(1-c)*28.45+y*c*(1-m) *8.19+c*m*(1-y) *2.51...

芜湖县15114582219: 求大神用MATLAB解一下三元三次方程,麻烦留下步骤啊,万分感谢!!!! -
朱帖福尔: 只含有一个未知数(即“元”),并且未知数的最高次数为3(即“次”)的整式方程叫做一元三次方程(英文名:one variable cubic equation).一元三次方程的标准形式(即所有一元三次方程经整理都能得到的形式)是ax3+bx2+cx+d=0(a,b,...

芜湖县15114582219: 求高手用matlab解一个三元的方程组 -
朱帖福尔: 我算的你这个解都是无穷大.具体过程如下:第一步:建立M文件函数myf.m function f=myf(x) f(1)=600+x(1)*cos(10)+x(2)*sin(40)-x(3); f(2)=(x(3)*sin(65)+x(2)*sin(40)+x(1)*cos(10))^2+(x(3)*cos(65)+x(2)*cos(40)+x(1)*cos(10))^2-5200^2; f(3)=x(1)*...

芜湖县15114582219: matlab中如何求解三元复杂方程组 ?程序如下 -
朱帖福尔: 解方程要求把已知数全代入,把未知数的个数降到最少,否则,多一个符号未知数,就多一点解不出的可能性.所以,你这个程序中,要把所有的已知数全代入,才可能解出来.

芜湖县15114582219: 用matlab求解三元一次方程组,要求不用自带命令slove,完全按照C语言的方式来写求代码或者思路, 比如2x+3y+4z=1,5x+6y+7z=2,9x+8y+z=3 这个怎么一... -
朱帖福尔:[答案] 雅克比迭代或赛德尔迭代 %jacobi.m雅可比迭代法 function x=jacobi(a,b,x0)%x0初值 D=diag(diag(a));%取对角元素阵 U=-triu(a,1);%取上三角 L=-tril(a,-1);%取下三角 B=D\(L+U); f=D\b; x=B*x0+f; n=1; while norm(x...

芜湖县15114582219: 用MATLAB求解三元非线性方程 -
朱帖福尔: syms x y z; [x,y,z]=solve('x^2+2*x+1=0','x+3*z=4','y*z=-1')x = -1 y = -3/5 z = 5/3

芜湖县15114582219: 跪求用Matlab解带指数的三元方程组 -
朱帖福尔: 先新建一个函数 function F = fun1( x ) v1=53.56; v2=24.46; v3=10.16; T1=323.13; T2=343.13; T3=373.13; F=[log10(v1+x(1))-x(2)+x(3)*log10(T1);log10(v2+x(1))-x(2)+x(3)*log10(T2);log10(v3+x(1))-x(2)+x(3)*log10(T3)]; end 再重新建个脚本clc x0=[...

芜湖县15114582219: 如何利用excel或者matlab解比较复杂的三元一次方程组 -
朱帖福尔: syms x y z ex1 = exp(-0.4)+20*exp(-0.2)*(2*sinh(0.2*y)-sinh(0.2))-x; ex2 = 0.5+2*(0.5-x)-y; ex3 = exp(-0.4)-80*exp(-0.2)*cosh(0.2*y)-z; [x,y,z] = solve(ex1,ex2,ex3,'x,y,z'); x = double(x) y = double(y) z = double(z) % 验证 subs(ex1)<=eps subs(ex2)<=eps ...

本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网