matlab二元方程组求解ppt

% by dynamic of Matlab技术论坛% 2010-02-26 18:57:08%首先回答,Matlab可以能求解二元高阶微分

solve('方程一','方程二','自变量')

代码：clearallclcf=@(t,y)([y(2);2*y(2)+.15*y(4)^2-10;   y(4);.01*y(4)^2]);[t,Y]=ode45(f,

二次方程组用fsolve求解吧,用solve可能接不出来呢,因为solve是解线性方程组的.再问：fsolve也试过了，依然无法求解呀clearP0=[1e26;1e26;1e26;1e26;1e26

没错啊,解不出来很正常,solve是符号工具箱的,用来求解析解,很弱的你这种竖直问题去看fsolve函数吧

anyhow,Igotananwser,asbelow.x=(.27313763081764851197485172598348e-22-.303486256464053902190953873151

functionhhX0=[17,195,1500,4,62,200];TSPAN=[04];[t,X]=ode45(@myfun,TSPAN,X0)plot(t,X)legend('x','y','

5x+6y=15.213x-2y=-0.421+2*3得到：5x+6y+9x-6y=15.2-0.4*314x=14x=1代入1式得到：5+6y=15.26y=10.2y=1.7解得：x=1y=1.7

解题思路：考察二元一次方程组的求解问题，根据方程组的特点，求出x+y的范围，再求解解题过程：解：两个方程两边分别相加得3x+3y=6-3m所以x+y=2-m因为x+y>-3/2所以2-m>

lz需要给出y_0的初值~这个方程满足李普希兹条件,因此,解存在唯一并且可以唯一延拓到边界,应用lax等价定理可以知道,向前欧拉法具有二阶的收敛速度……所以可以尝试用向前欧拉法编写：我刚学matlab

可以用ode45求数值解可以参考我给别人回答的再问：我是想对这个方程组做图。先解这个方程组，才能做出图来是不是？请问针对这个问题，m-file文件的程序该怎么编？再答：%%%%%%%%%%%%%fun

functiony=fc(x)y(1)=-19638+x(1)*(1-2*exp(-20/x(2)))y(2)=12767+x(1)*(1-2*exp(-6.5/x(2)))y=[y(1)y(2)]

clearall;clc;load15.mat;t=celiang(:,1);u1=celiang(:,2);symsu2i5iixu5%定义u1,L1,L2,R1,C1,C3,ii等vU1=celi

[X,Y]=solve('x^2+y^2=10^2*(1.568-1.538)','x*2=arctan(y/x)')X=-0.60651959867130434289594767388317Y=1.

f=@(x)[((x(1)*cos(45/180*pi)-cos(x(2)))./(x(1)*cos(45/180*pi)+cos(x(2)))).^2+((cos(45/180*pi)-x(1).*

[x,y]=solve('x^2+y^2=1','x*y=2','x','y')x=((15^(1/2)*i)/2+1/2)^(1/2)/2-((15^(1/2)*i)/2+1/2)^(3/2)/2(

例子：求解x1-0.5sinx1-0.3cosx2=0;x2-0.5cosx1+0.3sinx2=0;程序在文件编辑区建立待求方程组文件并保存为fun.m：functiony=fun(x)y=[x(1

把第三个式子也微分.>>[x,y,z]=dsolve('Dx=(-2)*A*x+u*y','Dy=2*A*x+2*u*z-(A+u)*y','Dx+Dy+Dz=0','x(0)=1,y(0)=0,z(

用solve可以解决,helpsolve看看用法,很简单的.

程序给你稍稍改动一下,本身没有问题,调试通过,最为关键的是初值的选取.强烈建议买一本程耿东院士的《工程结构优化设计基础》一书,是国内最具权威的关于优化方法入门的书籍,大连理工大学出版社.我根据自己的理