matlab dx dy dt,dx/dt=y,dy/dt=-sinx,求大神帮忙编一个MATL? 爱问知识人

建立m文件：

function dx=dfun(t,x) %函数名为dfun，参数为t与x

dx=[x(2);-sin(x(1))]; %以向量形式表示方程

输入：

clear

ts=-15:0。

05:15; %步长取0。

x0=[1,0]; %设定参数初值

options=odeset('reltol',1e-6,'abstol',1e-9); %提高精度

[t,x]=ode45(@dfun,ts,x0,options); %调用ode45计算

plot(x(:,1),x(:,2)),grid %作出y(x)图形

axis equal

gtext('fontsize{12}x'),gtext('fontsize{12}y') %标记字体x

但以上并非曲线y=f(x)的完整形状(调整ts的范围也无济于事)，原因是y为x的周期函数，而数值解只能求出初值附近的解

本题可以求出y=f(x)的解析表达式

由dx/dt=y，dy/dt=-sinx，得

dy/dx=(dy/dt)*1/(dx/dt)=-sinx/y

分离变量，积分得

y^2=2*cos(x) C，其中C为常数

代入初始条件y(1)=0，可求得C=-2*cos(1)

∴y^2=2*cos(x)-2*cos(1)，此式为原方程组的解析解

利用ezplot命令可绘制出完整图像

clear

syms x y

ezplot(y^2-2*cos(x) 2*cos(1),[-8,8,-3,3])

axis equal

axis([-8,8,-3,3])

grid on

另外，改变初值将得到不同的图形(为什么？请思考)，例如

初值改为：x=1,y=√[2*(cos(1) 1)]-10^(-5)

初值改为：x=1,y=√[2*(cos(1) 1)] 10^(-5)

。

全部

matlab dx dy dt,dx/dt=y,dy/dt=-sinx,求大神帮忙编一个MATL? 爱问知识人相关推荐

求大神帮忙，想要一元20次方程的解matlab代码,在线等,(解一元20次方程的解matlab代码)
求大神帮忙,想要一元20次方程的解matlab代码,在线等,(解一元20次方程的解matlab代码) 已知y是一个1*101的数组如下,函数如下,想求得x 135 130 134 128 129 13 ...
ubuntu14.04扩展屏幕后，打开matlab就死机。求大神帮忙解决！！！万分感谢
ubuntu14.04扩展屏幕后,打开matlab就死机.求大神帮忙解决!!!万分感谢
matlab里trim函数,matlab 中trim函数总出现Index exceeds matrix dimensions，求大神帮忙解决...
运行一个trim函数,可是总是错误 Index exceeds matrix dimensions. Error in trim>trim_alg (line 164) gg=[x(ix)-x0 ...
如何修改matlab中的语句,求大神帮忙看一下这个语句怎么改！！！！！
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 clc clear W=imread('e:\matlab\aaa\1.jpg'); imshow (W); I = rgb2gray(W); % 灰度处 ...
matlab数值很小出错,求大神帮忙解决一下，用MATLAB求解动力学数据总是出错~ - 计算模拟 - 小木虫 - 学术科研互动社区...
CODE: function KineticsEst5 % 动力学ODE方程模型的参数估计 % % % % The variables y here are y(1)=xB, y(2)=xoNB, y ...
matlab程序 surf算法,【求大神帮忙，surf算法源代码解析】
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 % Example, Corresponding points % Load images I1=imread('1.JPG'); I2=imread(' ...
部分最小二乘回归分析matlab,MATLAB如何用偏最小二乘回归法做预测？跪求大神帮忙 - 程序语言 - 小木虫 - 学术科研互动社区...
CODE: function question_16 %2016-12-26 clear;clc pz=xlsread('C:\Users\cirnotxm\Documents\EGDownloads ...
matlab 赋空值,未对输出参数赋值求大神帮忙解惑
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 x=imread('F:/flower.bmp'); wname='db5' [Lo_D,Hi_D] = wfilters(wname); lf=leng ...
matlab 番茄,正在做成熟番茄识别，遇到一个小问题，求大神帮忙
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 pic=imread('D:\机械视觉\新建文件夹\图片.jpg'); %载入图象 function result=avefilter(pic)%均值滤波 ...
matlab运行一直正忙,matlab一直显示正忙,没反应，求大神帮助
代码如下,运行一直处于正忙状态,求大神帮忙看看代码有什么问题:Q %% 定义 (x,t) 平面上的网格点坐标 clear;clc; dx=0.05; % x 方向的步长 dy ...

matlab dx dy dt,dx/dt=y,dy/dt=-sinx,求大神帮忙编一个MATL? 爱问知识人

matlab dx dy dt,dx/dt=y,dy/dt=-sinx,求大神帮忙编一个MATL? 爱问知识人相关推荐

最新文章

热门文章