如何用matlab求解一个二阶常系数微分方程组
一、理论解，用dsolve函数（在command window 中输入doc dsolve可以查看帮助）
示例：
代码：
[x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)')
二、数值解，用ode45，或ode23, ode15s其他函数。（在command window 中输入doc ode45可以查看帮助）
示例：
1. 建立M文件
function dy=vdp(t,y)
dy=[y(2);
1000*(1-y(1)^2)*y(2)-y(1);];
2. 调用MATLAB 函数ODE15S
[T,Y]=ode15s('vdp',[0 1000],[2 0]);
plot(T,Y(:,1));
图像为：
方程组求解要降阶才可以求解
如何用Matlab求解级数的和??
在Matlab中，我们可以使用内置函数来求解级数的和。这里我们以一个无穷级数为例，例如：求和 ∑(i=1 to ∞) (1\/i^2)。步骤如下：首先，我们需要定义一个函数来计算每一项的值。在这个例子中，我们的函数是 f(i) = 1\/i^2。然后，我们需要定义一个变量来存储级数的和。接下来，我们需要使用一...
怎么用matlab求方程的根?
用matlab迭代法求方程的根，其求解思路是这样的：1、创建迭代公式，即 x(k+1)=sqrt(10\/x(k)-4*x(k))2、确定初值，x(1)=1.5 3、使用while循环语句，进行迭代 4、当x(k+1)-x(k)<ε=1e-5，则x(k)为方程的根。实现代码：x(1)=1.5 k=1;while k<10000 x(k+1)=sqrt(10\/x...
matlab求逆
MATLAB是一款强大的数学计算软件，提供了丰富的数学函数库。在MATLAB中，可以使用“inv”函数来求矩阵的逆。具体方法是，直接在MATLAB的命令窗口中输入“inv”，其中A是要计算逆的矩阵。MATLAB会自动返回矩阵A的逆矩阵。需要注意的是，只有方阵才能求逆，对于非方阵，MATLAB会报错。三、使用步骤和注意事项 ...
如何用MATLAB求π的近似值?
如下：y=0;for n=1:1000000 xx=(-1)^(n-1)\/(2*n-1);if abs(xx)>10^(-6)y=y+xx;else break end end fprintf('π的近似值为%.6f ，共需迭代%d步\\n',y*4,n)。介绍 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。行矩阵运算、...
如何使用MATLAB求方程的解?
(a)如果希望利用由Solve得到的解计算表达式的值,可以利用取代运算符\/. ,这样Mathematica就会自动带入相应的值.(b)由于解就是列表,因此可以用Part或[[]]从列表中“提取”解.在下面两个例子中演示这些方法的使用.假设要计算下述方程所有根的平方和:x^6-21x^5+175x^4-735x^3+1624x^2-1764x+720...
用matlab计算在[-3π,3π]区间,求y=cos(x)*In(x^2)的值。
【计算代码】打开matlab软件，可以直接在命令窗口中输入下列语句：>>x=[-3*pi:pi:3*pi] ←创建x等差数列数组，公差为pi（π）>>y=cos(x).*log(x.^2) ←使用点运算符，计算其对应于x的y值 运行结果 【本题相关知识点】1、等差数列的建立。方法一：用x=[x0：d：xn]创建数组，x0是...
matlab中求n次方根的n该用什么函数?
简单来说就是，使用power(x,1\/3),或者简单的使用 x^(1\/3)即可。matlab中求平方根（二次方根）用sqrt，但没有专门表示三次方根的函数。计算三次方根等价于计算某个数的三分之一次方。所以可以使用power(x,1\/n)的方法(表示开n次方根)。举例：>> x=4 x =4 >> power(x,1\/3)ans =1....
在matlab中如何求积分,用什么函数?
首先要将 m,x,a,b 这四个变量定义为符号变量 syms m x a b;Fx = a*x^2;int(Fx,x,m,n)3、通过上面这个方法，就能够求得任意一个函数在给定区间的积分，如果想看到书写的格式，可以用pretty命令，这样显示更接近平常的表示方法。1、在matlab中，积分运算有多种方式，为了便于查看不同方式...
如何在MATLAB中使用power函数求平方?
1、第一步打开matlab，在命令行窗口输入a=3^2，如下图所示：2、第二步按回车键之后，可以看到结果为9，表示3的平方，如下图所示：3、第三步命令行窗口输入 b=power(4,2)，表示4的平方，如下图所示：4、第四步按回车键，可以看到结果为16，power函数可以表示是数字的几次方，后面是2就是平方...
如何使用matlab,求AX=B?
a\\b = inv(a)*b；b\/a = b*inv(a)；2。数组除法：A\/B表示A中元素与B中元素对应相除。MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。MATLAB是matrix&laboratory两个词的组合，意为...
运家安内： 一、理论解,用dsolve函数(在command window 中输入doc dsolve可以查看帮助) 示例:代码: [x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)')二、数值解,用ode45,或ode23, ode15s其他函数.(在command window 中输入...
平武县15820754927： matlab里ode45解决二阶常系数微分方程用法 - ？
运家安内： ^1、相关常数都没给,这里随便设置一些数值: m=1;c=2;k=3; y0=4;yp0=5; tf=10; dy=@(t,y)[y(2);-(c*y(2)+k*y(1))/m]; [t,y]=ode45(dy,[0 tf],[y0 yp0]); plot(t,y(:,1)) 2、function test() [t,y]=ode45(@func1, [0,1], [0;0;1;2;2;2]); figure(1); clf; plot(t, y); legend('x','...
平武县15820754927： 如何用matlab求解二阶常微分方程组 - ？
运家安内： ^mathematica的解法要不要…… sol = NDSolve[{D[x[t], t, t] == x[t]/(x[t] + y[t])^2,D[y[t], t, t] == y[t]/(x[t] + y[t])^2, x[0] == 1, x'[0] == 1,y[0] == 2, y'[0] == 1}, {x[t], y[t]}, {t, 0, 100}]; Plot[{x[t], y[t]} /. Flatten@sol, {t, 0, 100}]解的图形:
平武县15820754927： matlab怎么解二阶非齐次常系数微分方程 - ？
运家安内： >> y=dsolve('D2y+100*Dy+200*y=0.6','x') y = C2*exp(x*(10*23^(1/2) - 50)) + C3*exp(-x*(10*23^(1/2) + 50)) + 3/1000
平武县15820754927： 一类二阶常微分方程的几种解法 - ？
运家安内：[答案] 1、引言常微分方程有着深刻而生动的实际背景,它从生产实践与科学技术中产生,而又称为现代科学技术中分析问题与解决问题的一个强有力的工具.人们对二阶及以上微分方程(包括线性、常系数、隐性)的研究,产生了许多理论成果.如胡爱莲[1]...
平武县15820754927： 用matlab建立一个二阶系统模型,测试在不同输入下的响应 - ？
运家安内： 参考代码: % 用传递函数的形式创建二阶系统 1/(s^2+s+1) G=tf(1,[1 1 1]);% 求阶跃响应 subplot 221 step(G)% 求脉冲响应 subplot 222 impulse(G)% 求余弦输入响应 subplot 223 t = 0:0.1:10; u = cos(t); lsim(G,u,t);% 求正弦输入响应 subplot 224 t = 0:0.1:10; u = sin(t); lsim(G,u,t);都是很基础的东西,直接调函数就行,似乎没什么好解释的了,如果实在看不明白再问吧.
平武县15820754927： 如何用matlab求解二阶微分方程,以及程序实例 - ？
运家安内： 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/question/1112476448513751339.html(龙格-库塔方法)
平武县15820754927： matlab 中二阶常微分方程的数值解法？
运家安内： odefun=@(t,x)[x(2);3*x(2)-2*x(1)+1];[t,y]=ode45(odefun,[0:0.01:2],[1 0]);plot(t,y)[t y] 结果 y(0.5000)=0.7896 y= dsolve('D2y-3*Dy+2*y=1','Dy(0)=0','y(0)=1');>> yy =exp(t) - exp(2*t)/2 + 1/2>> feval(@(t)exp(t) - exp(2*t)/2 + 1/2,0.5)ans = 0.7896
平武县15820754927： matlab求解二阶常微分方程 - ？
运家安内： d2θ/ds2=cosθ-sinθ两边同时乘dθdθ/dsd(dθ/ds)=cosθdθ-sinθdθ(dθ/ds)^2=2sinθ+2cosθ+cds=dθ/√(2sinθ+2cosθ+c)这个积分不能用初等函数表示:在https://baike.baidu.com/item/超越积分/5292075?fr=aladdin第8个,所以方程的解也不能用初等函数表示.
平武县15820754927： 如何用Matlab求二阶导数 - ？
运家安内： 用MATLAB求导函数:diff如求函数f关于x的n阶导数可以写为:diff(f,x,n); 例如:求f=x^3=2*x^2+5的二阶导数,代码为: >> syms x; >> f = x^3 + 2*x^2 + 5; >> diff(f,x,2) 运行结果: 得到答案:6*x+4
