自学MATLAB（一）

最新推荐文章于 2024-05-01 11:25:27 发布

咸咸的盐啊

最新推荐文章于 2024-05-01 11:25:27 发布

阅读量118

点赞数

分类专栏： matlab自学文章标签： matlab 开发语言

本文链接：https://blog.csdn.net/m0_58174594/article/details/124834523

版权

matlab自学专栏收录该内容

1 篇文章 0 订阅

订阅专栏

1.常微分运算

$\frac{\mathrm{d} v}{\mathrm{d} x}= g-\tfrac{_{Cd}}{m}v^{2}$ （初始条件：t=0 v=0)

%2022.5.18

%1-求解微分方程通解
syms m g cd v(t);   %符号定义
eq=diff(v,t)==g-cd/m*v^2;   %matlab中的微分方程模型
cond=v(0)==0;   %初值条件
sol=dsolve(eq,cond)  %微分方程模型求解，通解

%2-符号表达式求值
m0=68.1;
g0=10;
cd0=0.25;
vt=subs(subs(subs(sol,'m',m0),'g',g0),'cd',cd0)%变量值代入

%3-循环求解在时间0-20内，每秒的自由落体速度
i=1;
for ti=0:20
   vval(i)=double(subs(vt,'t',ti));
    i=i+1;
end

%4-对0到20s内的速度进行可视化绘图
tf=0:20;

plot(tf,vval)  %绘图函数plot
title('0-20s内自由落体速度随时间变化曲线')
xlabel('时间t(s)')
ylabel('速度v（m/s)')
grid on   %网格线
hold on   %保持当前图像状态再添加
plot(12,vval(13),'r.',"markersize",10)  %在图像中标出12s时速度对应的点

%5-求蹦极运动员的极限速度
limS=double(limit(vt,'t',inf))

输出：