Is_More_ 2021-06-13 14:57 采纳率: 100%
浏览 38
已采纳

matlab解常微分方程

请问这是matlab解常微分方程的哪一种方法

clc

clear all

g=9.8;

m=1;

C=0.5;

r=0.05;

A=pi*r^2;

rou=1.2;

D=(rou*C*A)/2;

detat=0.01;

theta=(35/180*pi);

x(1)=0;

y(1)=0;

V(1)=50;

Vx(1)=V(1)*cos(theta);

Vy(1)=V(1)*sin(theta);

N=700;

for n=1:N

    V(n)=sqrt(Vx(n)^2+Vy(n)^2)

    ax(n)=-(D/m)*V(n)*Vx(n);

    ay(n)=-g-(D/m)*V(n)*Vy(n);

    Vx(n+1)=Vx(n)+ax(n)*detat;

    Vy(n+1)=Vy(n)+ay(n)*detat;

    x(n+1)=x(n)+Vx(n)*detat+0.5*ax(n)*detat^2;

    y(n+1)=y(n)+Vy(n)*detat+0.5*ay(n)*detat^2;

    px(n)=Vx(1)*n*detat;

    py(n)=Vy(1)*n*detat-0.5*g*(n*detat)^2;

end

plot(x,y,'r',px,py,'g');

grid

xlabel('x'),ylabel('y');

title('有空气阻力的抛射体运动')

  • 写回答

1条回答 默认 最新

  • CSDN专家-Matlab_Fans 2021-06-13 17:22
    关注

    这是向前欧拉法求解常微分方程的数值解法,具体算法可以参考我的博客:

    https://blog.csdn.net/MatlabFans_Mfun/article/details/114198233

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论

报告相同问题?

悬赏问题

  • ¥15 求micropython modbus-RTU 从机的代码或库?
  • ¥15 铜与钢双金属板叠加在一起每种材料300mm长,18mm宽,4mm厚一端固定并加热至80℃,当加热端温度保持不变时另一端的稳态温度。ansys
  • ¥15 django5安装失败
  • ¥15 Java与Hbase相关问题
  • ¥15 后缀 crn 游戏文件提取资源
  • ¥15 ANSYS分析简单钎焊问题
  • ¥20 bash代码推送不上去 git fetch origin master #失败了
  • ¥15 LOL外服加入了反作弊系统,现在游戏录像rofl文件离线都无法打开
  • ¥15 在centos7安装conda
  • ¥20 WPF 如何实现多语言,label 和cs(live Charts)中是否都能翻译