matlab文件命名

这个m文件应该命名为什么呢？一直运行不了

clear
clc
close all

%结构参数（国际单位制）
E=206e9;
density=7800;
h=0.008;%梁的厚度
b=0.0557;%梁的宽度
l=0.609;%梁的长度

%划分网格个数
n=20;
%有限元法得到结构质量矩阵于刚度矩阵
[K,M]=finite(E,density,h,b,l,n);
%使用matlab的函数eig计算，用于验证自编程的正确性
[phi,omega]=eig(K,M);
omega=sqrt(diag(omega));
%构造比例阻尼
a=2*omega(1,1)*omega(2,1)0.01/(omega(1,1)+omega(2,1));
C=aM;

%外激励振幅
F0=10;
F=zeros(2n,1);
F(2n-1,1)=F0;
%激励力频率
omega_excitation=(omega(1,1)+omega(2,1))/2;
T_excitation=2*pi/omega_excitation;

dt=T_excitation/500;%步长
nt=500;%时间步
M_mode=phi'Mphi;%模态质量矩阵
K_mode=phi'Kphi;
C_mode=phi'Cphi;
F_mode=phi'*F;

%定义无阻尼初始条件（使伴生自由振动响应为0）
k=diag(K_mode);
v1=F_mode./k./(1-(omega_excitation./omega).^2)omega_excitation;%负的伴生自由振动模态位移幅值向量
v0=phiv1;%物理空间位移，将其作为初始位移
x0=zeros(2n,1);%初始位移
a0=M(F-Kx0);%初始加速度

rou1=1;
rou2=0;
y=CH_alpha(M,K,F0,x0,v0,a0,omega_excitation,rou2,dt,nt,n);

%模态叠加法佐证
t=0:dt:ntdt;
Hnn=0;
for k=1:2n
Hnn=Hnn+phi(2n-1,k)phi(2n-1,k)/(K_mode(k,k)-M_mode(k,k)omega_excitation^2);
end
Un=HnnF0sin(omega_excitation.*t);

figure(1)%无阻尼CH-alpha方法和模态叠加法
plot(t,y(2*n-1,:),'r',t,Un,'b--');

Undamped_CH-Alpha_Modal_superposition.m