- 博客(17)
- 收藏
- 关注
原创 麦克斯韦速率分布
代码clc;clear;figure('numbertitle','off','name','麦克斯韦速率分布');global K T m;K=10;T=300;V=0:0.01:32;for m=50:50:250;subplot(1,2,1);FF1=4*pi*(m/(2*pi*K*T)).^(3/2);FF2=V.^2;FF3=exp(-m*V.^2/(2*K*T));% 坐标轴美化FF=FF1*FF2.*FF3;plot(V,FF,'linewidth',2.5
2022-11-02 01:57:19 2067
原创 2022数学建模A动态仿真
%% 2022数学建模A动态仿真clc;clear;figure('numbertitle','off','name','2022数学建模A动态仿真');%赋值global m M K C W C1 P Fm=4800;M=2400;K=80000;C=10000;W=1.8;C1=656;P=10000;%改变F=6000;%ode45求解绘图t1=[0:0.1:40];H1=[0;10;0;5];[t,H]=ode45(@dH,t1,H1);
2022-09-17 13:59:23 194
原创 有限差分法求解高阶导数
有限差分法的基本思想是把微分方程进行离散化,在运用边界条件,用一个一个不确定的差分方程进行逐级迭代去近似替换。有限差分法的实质是把特殊的偏微分方程改变成为有限差分形式迭代方程,把特殊的偏微分方程用线性代数微分方程组替换,以及把数学推导过程、改变成为数值计算,进而使连续特殊偏微分方程进行离散化。关于差分法主要讨论的问题有:差分格式的建立,把特殊的偏微分方程用差分方程进行代替,边界条件设置,以及差分矩阵相关概念理解与使用MATLAB编写相应数值分析求解的代码。在有限差分算法中,步长一般称为离散化之后的单元的
2022-07-29 18:47:23 1830
原创 三个弹簧振子理论
三个弹簧振子理论建立如图 1 所示的一维三弹簧振子的实验系 统:三个弹簧振子的质量均为m1,m2,m3,弹簧的劲度系 数分别为k1,k2和k3 ,三个谐振子在微振动过程 中相对自身的平衡位置的位移分别为x1、x2 和 x3 ,取x1、x2 和 x3 广义坐标,则系统的动力学方程可表示为:代码弹簧摆摆动图clc;clear;global k1 k2 k3 M1 M2 M3k1=145000;k2=180000;k3=121000;M1=220;M2=240;M3...
2022-07-27 23:22:38 1201 1
原创 经纬度下温度变化
clc;clear;data1=xlsread('1998tmp.xlsx');aa1=xlsread('AAA1.xlsx');aa2=xlsread('AAA2.xlsx');aa3=xlsread('AAA3.xlsx');aa4=xlsread('AAA4.xlsx');aa5=xlsread('AAA5.xlsx');aa6=xlsread('AAA6.xlsx');aa9=xlsread('AAA9.xlsx');A1=data1(:,1);a1=A1(2:end);B1=d
2022-07-21 21:58:45 233
原创 弹簧摆摆动轨迹
%% 弹簧摆摆动轨迹clc;clear;global L m k gtheta0=pi/10;m=1;k=80;g=9.8;L0=1;L=L0+m*g/k;t1=[0:0.05:15];H1=[L0 0 theta0 0];[t,H] = ode45(@dH,t1,H1);%求解动力学方程figure('numbertitle','off','name','弹簧摆摆动轨迹图');[y1,x1]=pol2cart(H(:,3),-H(:,1));plot(x1(1),y1(1)..
2022-06-23 01:20:44 402
原创 小珠子在环内的凹槽中滚动仿真图
%% 小珠子在环内的凹槽中滚动仿真图clc;clear;figure('numbertitle','off','name','小珠子在环内的凹槽中滚动仿真图');global K RK=4;w=1.3;R=0.2;r=1;t1=[0,80];H1=[0;0.3];[t,H] = ode45(@dH,t1,H1);z=r-cos(H(:,1));x=sin(H(:,1)).*cos(w*t);y=sin(H(:,1)).*sin(w*t);plot3(x,y,z,'b',.
2022-06-23 01:18:30 168
原创 Matlab相图
%% 水平动力学运动的相图clc;clear;figure('numbertitle','off','name','水平动力学运动的相图');global K RK=5;R=0.5;t1=[-10:0.1:10];for i=[-0.5 -1.5 -3 -4.5 0.5 1.5 3 4.5] H1=[0;i]; [t,H] = ode45(@dH,t1,H1); plot(H(1:end,1),H(1:end,2)); axis([-5 5 -4 4]...
2022-06-23 01:12:52 1047
原创 水平有摩擦力K的变化对运动曲线的影响
%% 水平有摩擦力K的变化对运动曲线的影响clc;clear;global K Rfigure('numbertitle','off','name','水平有摩擦力K的变化对运动曲线的影响');for K=0:0.5:3 R=0.25; t1=[0,80]; H1=[0;0.1]; [t,H] = ode45(@dH,t1,H1); plot(t,H(:,1),'LineWidth',1.5); hold on axis([0 60...
2022-06-23 01:06:35 83
原创 串联弹簧动态物理仿真
%% 弹簧动图clc;clear;global k1 k2 k3 M1 M2 M3k1=145000;k2=180000;k3=121000;M1=220;M2=240;M3=350;t1=[0,5];H1=[4;0;8;0;12;0];[t,H] = ode45(@dH,t1,H1);%求解动力学方程figure('numbertitle','off','name','弹簧摆摆动轨迹图');R=0.3;xx1=20;xx2=40;xx3=60...
2022-06-22 02:14:49 180
原创 撒克逊碗动态物理仿真
%% 撒克逊碗clc;clear;global M g p R r hd N N2 Lg=9.8;M=0.13581;R=0.05;p=1000;%千克%水的密度r=0.003;hd=0.002;N=0.0001;%碗的厚度%阻尼系数小孔的N2=1;L=0.1;%碗的阻力系数%碗的高度t1=[0,40];H1=[0;0;0;0][t,H] = ode45(@dH,t1,H1);%求解动力学方程%% 动画绘制figure('numbertitle','off','name','
2022-06-22 02:06:23 215
原创 等势线动图
clc;clear;%%等势线动图for a=2:1:14;q1=-1.602e-19;q2=1.602e-19;radius = 2; %点电荷半径x = 25; y = 25; %x和y轴显示边界值k = 9e+9; %k = 1/4πεdelta = 1;%引入距离误差以消除点电荷处点位无限大的问题phi = linspace(-pi,pi,20); %电场线发射范围[X,Y]=meshgrid(-x:0.5:x,-y:0.5:y);%获取左半部分矩阵r1=sqr...
2022-06-20 02:09:31 165
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人