微分方程5_如何理解$e^i*pi

在之前的微积分学习中,通过实例了解到了 e t e^t et的性质

  1. e t e^t et的导数等于它本身,那么 d / d t ( e a t ) = a ∗ e a t d/dt(e^{at})=a*e^{at} d/dt(eat)=aeat

  2. 任何数的0次方等于1

在这里插入图片描述

假如函数 e t e^t et描述了数轴上的位置,它随着事件而变化y= e t e^t et,t代表时间,输出值y等于位置
在这里插入图片描述
那么位置将从1开始,因为t=0时y=1
在这里插入图片描述
因为 e t e^t et的导数等于本身,任何时候速度的值等于位置本身

在这里插入图片描述
通过位置与速度的关系,就可以只管的理解这个函数关系
在这里插入图片描述
当指数是常数倍的时候,速度就等于位置的常数倍
在这里插入图片描述
在这里插入图片描述
如果指数的倍数是负数,表示先翻转向量,再缩放
在这里插入图片描述
在这里插入图片描述

当倍数常数是代表-1开方的i,
在这里插入图片描述

在这里插入图片描述

在这里插入图片描述
旋转90°,已经超出了数轴,我们用复平面来表示

在这里插入图片描述

在这里插入图片描述

这个向量,就是 a + b ∗ i a+b*i a+bi

它的初始状态 e i ∗ 0.00 = 1.00 + 0.00. ∗ i e^{i*0.00}=1.00+0.00.*i ei0.00=1.00+0.00.i

在这里插入图片描述

1个单位时间后,速度为1个单位,走过了一个圆
在这里插入图片描述

π \pi π秒之后,到达了x轴负半轴,所以 e i π = − 1 + 0.00 i = − 1 e^{i\pi}=-1+0.00i=-1 eiπ=1+0.00i=1
在这里插入图片描述
τ \tau τ秒之后,( τ \tau τ=2 π \pi π),走了一个圆周的距离,
e i τ = 1 + 0.00 i = 1 e^{i\tau}=1+0.00i=1 eiτ=1+0.00i=1

在这里插入图片描述
由此可以看出,复指数可以看着做单位弧度
在这里插入图片描述

在这里插入图片描述

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
我们可以使用 MATLAB 中的数值求解函数 ode45 来求解这个微分方程组。首先,需要将微分方程组转化为函数形式,具体实现代码如下: ```matlab function dydt = water_evaporation(t, y, k, RH0) r = (3*y(1)/(4*pi))^(1/3); % 求解当前时刻水滴半径 S = 4*pi*r^2; % 求解当前时刻水滴表面积 V = 4/3*pi*r^3; % 求解当前时刻水滴体积 RHt = y(2); % 当前时刻相对湿度 T = 25+273.15; % 温度设定为 25 度 esat = 611.2*exp(17.67*(T-273.15)/(T-29.65)); % 计算当前时刻的饱和水汽分压 et = esat*RHt; % 计算当前时刻水汽的分压 dydt = zeros(2,1); % 初始化返回的微分方程组值 dydt(1) = -k*S*sqrt(RHt/RH0)*et/V; % 计算体积的变化率 dydt(2) = -2*k*S*sqrt(RHt/RH0)*et/esat; % 计算相对湿度的变化率 end ``` 在上述函数中,t 和 y 分别表示当前时刻和对应的状态变量,k 和 RH0 分别表示水的蒸发速率常数和空气的初始相对湿度。在函数中,我们首先根据当前时刻水滴半径计算出水滴的表面积和体积,然后根据当前时刻的相对湿度和温度计算出饱和水汽分压和水汽分压。接着,根据微分方程组的定义,计算出体积和相对湿度的变化率,并将其保存在 dydt 中返回。 接下来,我们可以使用 ode45 函数来求解微分方程组。具体实现代码如下: ```matlab k = 1e-10; % 设置水的蒸发速率常数 RH0 = 0.8; % 设置初始相对湿度 r0 = 1e-3; % 设置初始水滴半径 tspan = [0 3600]; % 设置求解时间范围 y0 = [r0^3*pi*4/3 RH0]; % 设置初始状态变量 [t, y] = ode45(@(t,y)water_evaporation(t,y,k,RH0), tspan, y0); % 求解微分方程组 ``` 在上述代码中,我们首先设置了水的蒸发速率常数和初始相对湿度,并设定了初始时刻的水滴半径、求解时间范围和初始状态变量。接着,使用 ode45 函数来求解微分方程组,其中第一个参数 @(t,y)water_evaporation(t,y,k,RH0) 表示微分方程组的函数句柄,tspan 表示求解时间范围,y0 表示初始状态变量。求解完成后,t 和 y 分别表示求解出的时间和对应的状态变量。 最后,我们可以使用 plot 函数将求解结果可视化。具体实现代码如下: ```matlab r = (3*y(:,1)/(4*pi)).^(1/3); % 计算对应时刻的水滴半径 V = y(:,1); % 计算对应时刻的水滴体积 RH = y(:,2); % 计算对应时刻的相对湿度 figure; % 新建一个图形窗口 subplot(2,2,1); plot(t, r*1e6); xlabel('Time (s)'); ylabel('Radius (\mum)'); % 绘制水滴半径随时间变化的曲线 subplot(2,2,2); plot(t, V*1e9); xlabel('Time (s)'); ylabel('Volume (nL)'); % 绘制水滴体积随时间变化的曲线 subplot(2,2,3); plot(t, RH); xlabel('Time (s)'); ylabel('Relative Humidity'); % 绘制相对湿度随时间变化的曲线 subplot(2,2,4); plot(V*1e9, RH); xlabel('Volume (nL)'); ylabel('Relative Humidity'); % 绘制相对湿度与水滴体积的关系 ``` 在上述代码中,我们首先根据求解结果计算出对应时刻的水滴半径、体积和相对湿度,并将其分别保存在 r、V 和 RH 中。接着,使用 subplot 函数在同一个图形窗口中绘制出水滴半径、体积和相对湿度随时间变化的曲线,同时绘制出相对湿度与水滴体积的关系图。 完整代码如下: ```matlab function dydt = water_evaporation(t, y, k, RH0) r = (3*y(1)/(4*pi))^(1/3); % 求解当前时刻水滴半径 S = 4*pi*r^2; % 求解当前时刻水滴表面积 V = 4/3*pi*r^3; % 求解当前时刻水滴体积 RHt = y(2); % 当前时刻相对湿度 T = 25+273.15; % 温度设定为 25 度 esat = 611.2*exp(17.67*(T-273.15)/(T-29.65)); % 计算当前时刻的饱和水汽分压 et = esat*RHt; % 计算当前时刻水汽的分压 dydt = zeros(2,1); % 初始化返回的微分方程组值 dydt(1) = -k*S*sqrt(RHt/RH0)*et/V; % 计算体积的变化率 dydt(2) = -2*k*S*sqrt(RHt/RH0)*et/esat; % 计算相对湿度的变化率 end k = 1e-10; % 设置水的蒸发速率常数 RH0 = 0.8; % 设置初始相对湿度 r0 = 1e-3; % 设置初始水滴半径 tspan = [0 3600]; % 设置求解时间范围 y0 = [r0^3*pi*4/3 RH0]; % 设置初始状态变量 [t, y] = ode45(@(t,y)water_evaporation(t,y,k,RH0), tspan, y0); % 求解微分方程组 r = (3*y(:,1)/(4*pi)).^(1/3); % 计算对应时刻的水滴半径 V = y(:,1); % 计算对应时刻的水滴体积 RH = y(:,2); % 计算对应时刻的相对湿度 figure; % 新建一个图形窗口 subplot(2,2,1); plot(t, r*1e6); xlabel('Time (s)'); ylabel('Radius (\mum)'); % 绘制水滴半径随时间变化的曲线 subplot(2,2,2); plot(t, V*1e9); xlabel('Time (s)'); ylabel('Volume (nL)'); % 绘制水滴体积随时间变化的曲线 subplot(2,2,3); plot(t, RH); xlabel('Time (s)'); ylabel('Relative Humidity'); % 绘制相对湿度随时间变化的曲线 subplot(2,2,4); plot(V*1e9, RH); xlabel('Volume (nL)'); ylabel('Relative Humidity'); % 绘制相对湿度与水滴体积的关系 ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值