mnxzd.txt

来自「拟陀螺运动程序,该程序用MATLAB编写,如有不足之处,希望指出.」· 文本 代码 · 共 36 行

摘要：本程序模拟弦振动过程。

程序：

% 弦振动模拟
clc;clear;close all;
figure;
set(gcf,'DoubleBuffer','on');
% \copyright: zjliu
% Author's email: [e-mail] 
x=linspace(0,2,100);
u0=x*2;
u0(x>1)=4-x(x>1)*2;
Hux=plot(x,u0,'k');
axis([0,2,-4,4]);
hold on;
vx=zeros(size(x));
dx2=[x(2)-x(1)]^2;
dt=0.02;
a=1;
Tx=title('模拟弦振动, 时间: t=0');
t=0;
vxt=vx(2:end-1);
while t<20;
    t=t+dt;
    D2=[diff(u0(2:end))-diff(u0(1:end-1))]/dx2;
    vt=vxt+D2*dt*a;
    vxt=vt;
    ut0=dt*vt+u0(2:end-1);
    u0=[0,ut0,0];
    set(Hux,'YData',u0);
    set(Tx,'String',['模拟弦振动, 时间: t=',num2str(t)]);
    pause(0.1);
end