2))<eps dx(2)=0;else dx(2)=450/sqrt(1+((90*t-x(1))/(120-x(2)))^2);end end matlab只能给出数值解,所谓数值解也就是指定有限点上的函数取值,一般没法用显式的函数表达式来描述.将上述代码保存为solv.m,运行[t,x]=solv([0:0.01:1])便能得到[0:1]区间上函数取值.
用数值方法吧，把t=0时的y值输入，即可解出来。function hh [t,Y]=ode45(@myfun,[0:0.1:1],[0 0 0 0])function dy=myfun(t,y)y1=y(1);y2=y(2);y3=y(3);y4=y(4);dy1=y2;dy2=0.7143*(y1*y4*y4-9.81*sin(y3));dy3=y4;dy4=(280.2913*y4*cos(y3)+420.437...
由题目可知，自变量是t，因变量是x,y,z，化为matlab时是x(1),x(2),x(3)，不会有x(4)--x(6),说明函数式有问题。
这是二阶非线性方程用dsolve（）求解无法得到其解析解。应考虑用数值方法（如ode45）来求解。求解时还得提供该微分方程的初始条件。
首先可以肯定的告诉你，matlab本身有个pde工具箱，可以求解一般的典型二阶线性偏微分方程。但是，对于非线性的偏微分方程，我估计很难直接用matlab的工具箱或者函数命令得到结果。因为这涉及到一个数学处理的问题，你首先要处理这个非线性问题，或者用数值的方法处理，然后再用MATLAB编程计算。
然后求解即可。也就是说，在solve那一句之前加上：eq2= subs(eq2,{'bb','dd','ee'},{i1(j-1), u2(j-1),fai(j-1)});eq4= subs(eq4,{'bb','dd','ee'},{i1(j-1), u2(j-1),fai(j-1)});不过，这样又带来一个新的问题：求解7次方程会得到7个解，你要保留哪一个？
用MATLAB二分法求方程2x²+4x-5的解，其方法如下：1、自定义函数 dichotomy_fun=@(x)2*x^2+4*x-5;2、用二分法函数——dichotomy（），求解在x=(0.5,1.5)之间的值x=dichotomy(dichotomy_fun,0.5,1.5,1e-8)3、运行结果 n = 27 %运算次数 x = 0.870828695595264 %x运算...
例如：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')输出计算结果为...
21372/625-21372/625*XA)/(89113/250+37477/1000*XA)+20667888735911619/9007199254740992000*exp(-16430000/1987/T)*(31801/1000*XA+13/625)/(89113/250+37477/1000*XA))^2-1288029493427961856/17149319079463875*T+5945222134290114936832/85746595397319375 给出初始条件，用ode45求解。这是动力学问题吧？
可以使用Matlab内置函数fsolve来求解非线性方程组。具体步骤如下：定义一个匿名函数，将两个方程表示为一个向量。F = @(x) [24.9697*x(2)/(sqrt((24.9697-x(1))^2+x(2)^2))-24.969; 25.0282*x(2)/(sqrt((25.0282-x(1))^2+x(2)^2))-25.029];使用fsolve函数求解方程组。