matlab计算代数方程组(Matlab如何求解代数方程)
本文目录
- Matlab如何求解代数方程
- 怎么用matlab 解代数方程组
- 如何用MATLAB计算方程组
- matlab 求代数方程组 a*x^2+b*x+c=0 x+y=0 关于x,y的解,并分别绘制x和y关于b和c的图像(a视作常数)
- matlab 求代数方程组解
Matlab如何求解代数方程
》》 syms x》》 syms y》》 solve(’x^2+y^2=2’,y)
怎么用matlab 解代数方程组
=solve(’x^2+y^2=d^2’,’(x-a)^2+(y-sqat(d^2-a^2))^2=c^2’)
如何用MATLAB计算方程组
syms x1 x2 x3 x4 x5 f1 f2 f3 f4 f5f1=(x1+x2+x3+x4+x5+22.55*x1*x4+7.63*x3*x4+148.222*x1^3*x4+5.86*10^13*x1^12*x4^7+67.03*x1*x4^2+212.17*x1*x4^6+338.08*x1*x5+4948932.908*x1^3*x5^2+3522485477*x1^4*x5^3+3.46*x4*x5+4.109*x3*x5+6.03*x3*x5^2+4.847*x3^2*x5+135498.74*x1^3*x2*x4^3+7.55*10^14*x1^11*x2*x4^7)-1;f2=0.5653*(1/3*x2+135498.74*x1^3*x2*x4^3+7.55*10^14*x1^11*x2*x4^7)-0.2946*(0.5*x1+22.55*x1*x4+444.666*x1^3*x4+7.03*10^14*x1^12*x4^7+67.03*x1*x4^2+212.17*x1*x4^6+338.08*x1*x5+14846798.72*x1^3*x5^2+1.41*10^10*x1^4*x5^3+406496.24*x1^3*x2*x4^3+8.305*10^15*x1^11*x2*x4^7);f3=0.2946*(0.5*x3+7.63*x3*x4+4.109*x3*x5+6.03*x3*x5^2+9.694*x3^2*x5)-0.05*(1/3*x2+135498.74*x1^3*x2*x4^3+7.55*10^14*x1^11*x2*x4^7);f4=0.05*(x4+22.55*x1*x4+7.63*x3*x4+148.222*x1^3*x4+4.102*10^14*x1^12*x4^7+134.06*x1*x4^2+1273.02*x1*x4^6+3.46*x4*x5+406496.24*x1^3*x2*x4^3+5.285*10^15*x1^11*x2*x4^7)-0.3456*(0.5*x3+7.63*x3*x4+4.109*x3*x5+6.03*x3*x5^2+9.694*x3^2*x5);f5=0.3456*(x5+338.08*x1*x5+9897865.816*x1^3*x5^2+1.056*10^10*x1^4*x5^3+3.46*x4*x5+4.109*x3*x5+12.06*x3*x5^2+4.847*x3^2*x5)-0.05*(x4+22.55*x1*x4+7.63*x3*x4+148.222*x1^3*x4+4.102*10^14*x1^12*x4^7+134.06*x1*x4^2+1273.02*x1*x4^6+3.46*x4*x5+406496.24*x1^3*x2*x4^3+5.285*10^15*x1^11*x2*x4^7);=solve(f1==0,f2==0,f3==0,f4==0,f5==0,x1,x2,x3,x4,x5);
matlab 求代数方程组 a*x^2+b*x+c=0 x+y=0 关于x,y的解,并分别绘制x和y关于b和c的图像(a视作常数)
1. 求代数方程组的解:
》》 =solve(’a*x^2+b*x+c=0’,’x+y=0’,’x,y’)
x =
1/2/a*(-b+(b^2-4*a*c)^(1/2))
1/2/a*(-b-(b^2-4*a*c)^(1/2))
y =
-1/2/a*(-b+(b^2-4*a*c)^(1/2))
-1/2/a*(-b-(b^2-4*a*c)^(1/2))
》》
2. 从上面的解可以看出,x,y都有两组解且x,y互为相反数。
假设a=1,这里有两种方法绘制x,y关于b,c的图像:
(1)隐函数绘图
x1=subs(x(1),’a’,1);
x2=subs(x(2),’a’,1);
y1=subs(y(1),’a’,1);
y2=subs(y(2),’a’,1);
figure
po=get(gcf,’position’);
set(gcf,’position’,);
subplot(121)
ezsurf(x1,)
hold on
ezsurf(x2,)
subplot(122)
ezsurf(y1,)
hold on
ezsurf(y2,)
(2)根据方程式直接绘图
》》 a=1;
》》 =meshgrid(-10:0.5:10);
》》 delta=b.^2-4*a*c;
》》 delta(delta《0)=NaN;
》》 x1=0.5/a*(-b+sqrt(delta));
》》 x2=0.5/a*(-b-sqrt(delta));
》》 y1=-x1;
》》 y2=-x2;
》》 figure
》》 po=get(gcf,’position’);
》》 set(gcf,’position’,);
》》 subplot(121)
》》 surf(b,c,x1)
》》 hold on
》》 surf(b,c,x2)
》》 xlabel(’b’);ylabel(’c’);zlabel(’x’)
》》 subplot(122)
》》 surf(b,c,y1)
》》 hold on
》》 surf(b,c,y2)
》》 xlabel(’b’);ylabel(’c’);zlabel(’y’)
matlab 求代数方程组解
可以用solve函数求解。(注意,该函数在高版本中使用有一点问题)
syms y1 y2 h1 h3 l1 l2
syms x1 x2
=solve(y1+y2*x2-h3*x*2^2-2*h1*x1+(l1+l2)*x2==0,y2+y1*x1-h1*x1^2-2*h3*x2+(l1+l2)*x1==0,x1,x2)
上述代码在matlab2016a通过。答案有点长。
更多文章:
美参议院批准芬兰瑞典加入北约(芬兰和瑞典正式确认加入北约的意愿,这会有哪些影响)
2024年8月2日 11:41
asp教程连接数据库(用asp.net如何连接sql数据库)
2024年4月1日 07:38
shell命令超过内存提醒(利用shell弄cpu内存磁盘的优势)
2024年7月16日 18:16
gettickcount(gettickcount 返回值为负)
2024年7月23日 21:17
支付宝小程序开发文档(支付宝小程序怎么使用,入口在哪里,和微信小程序开发环境是一样的吗)
2024年7月17日 06:09
bigdecimal累加求和(Java新手求100阶乘到200的阶乘相加的Java程序)
2024年3月19日 02:00
根据excel生成数据库脚本(怎样通过脚本文件往SignalBuilder里面导入Excel数据)
2024年5月2日 03:43
tsql语言中有什么运算(T-SQL 两个binary值作位运算)
2024年7月21日 03:14
matlab绘制极坐标函数(使用matlab绘制由极坐标描述的曲线)
2024年6月29日 11:53
postman中文意思(mybrotherisapostman中文意思是)
2024年7月24日 12:40
metaphors(metaphor 和 analogy的区别)
2024年6月26日 19:55