函数中必须显式给出 dy(1)=...dy(2)=...dy(3)=...———函数没问题了，剩下的问题是 [t,y]=ode45(@weifensan,[0,3000],[10000;10000;0])；plot3（y(:,1),y(:,2),y(:,3)）分号和括号要在英文半角状态下输入，而你用的是中文 ———hold on ...
用Matlab求解偏微分方程，可以用pde工具箱来解决。例如，简单的点热源方程 求解步骤及主要函数：1、问题定义 2、创建具有单个因变量的PDE模型， createpde（）3、创建几何结构并将其追加到PDE模型中，geometryFromEdges（）4、使用边界条件，pdegplot（）5、指定PDE系数 6、指定初始条件，setInitial...
你好，求解这个方程的一般步骤用ode45。代码如下 clc;clear m = 0.00267; alpha = 4e-6;beta = 1e-6; U = 90;k = 0.5; p = 19.5216;A = 0.001256;dy=@(t, x, m, alpha, beta, U, k, p, A)[x(2);...m^(-1)*(-(k+1/beta)*x(1)+alpha*U/beta-p*...
syms a,b,c,d； （x）clc;clear all;syms a b c d;dsolve('Dy(1)=1/2*a*y(2)*j','Dy(2)=-(c+d)*y(2)*j+1/2*a*y(1)*j+1/2*b*y(3)*j','Dy(3)=-d*y(3)*j+1/2*b*y(2)*j')error('symbolic:dsolve:errmsg2','There are more ODEs than variables.')...
1、对于解析值，你可以用dsolve()函数求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解 >> syms y(x)>>Dy = diff(y)；D2y = diff(y, 2)；>>dsolve(x*D2y+x*(Dy)^2-Dy==0,'x')运行结果 2、对于数值解，你可以查看二阶常微分方程这个例子，https://zhidao.baidu.com/...
3.6 用dsolve函数可接 3.7 t=ln3
用matlab的dsolve函数计算该二阶微分方程的解析值是有一定的难度。只能用数值分析的方法来求解其数值解。当然求其数值解，必须还得已知a、θ的具体数值。该问题可以用matlab的ode45函数求得。求解方法如下：1、建立微分方程的自定义函数，odefun（x,y）2、当a=10、θ=π/6时，执行下列代码 theta=pi...
题主给出的二阶微分方程组，用matlab求解其精确解和数值解，可以使用dsolve（）函数得到其解析值，使用ode（）函数得到其数值解。一、使用dsolve（）函数求解 第一步，对x(t)、y(t)进行声明变量，即 syms x(t) y(t)第二步，使用diff（）函数分别对x(t)、y(t)进行导数计算 Dy=diff(y,t,1...
设置问题的类型。如 图，选择菜单中的options，下拉菜单中选择application。可以看到偏微分方程适用几乎所有常见数学问题类型。选择你想要求解的一类【这个 一定要选择，因为后面的方程类型和边界条件，matlab都会按照你选择的类型帮你做好初始化，你只需要动手改改参数就可以了。】8 设置方程的类型 选择菜单中...
r2*(1-N(2)/K2)-n*N(1))*N(2);dNdt = [dN1dt;dN2dt];运行结果：t = 0 0.2000 0.4000 0.6000 0.8000 1.0000 N = 0 0 0 0 0 0 0 0 0 0 0 0 再检查一下方程或初始条件吧，两种方法结果一样。如果初值不为0，N的结果不是0。