常微分方程一般可以用ode函数来求其数值解,然后根据z和ρ的数值绘制其z(ρ)的函数图形。
实现方法:
1、根据常微分方程,建立其自定义函数
fun=@(rho,z)sqrt((0.5+0.5/sqrt(1-rho^2))^2-1)
2、利用matlab的ode函数求解。
[rho,z]=ode45(fun,tspan,z0);
3、根据z和ρ的数值绘制其z(ρ)的函数图
plot(rho,z);
xlabel('ρ'),ylabel('z(ρ)')
编程后运行,可得到如下z(ρ)的函数图
一阶的比较简单,设置一下函数就可以绘制出来