Matlab玻尔共振数据计算-作图程序

%“自由振动”——表1的数据计算：红色部分大家填入自己记录的数据。%表示后面是解释说明语句，matlab不执行。大家可以在红色部分填入自己数据后，依次复制到matlab命令窗口去执行；也可以同时复制到matlab命令窗口一次执行得到所有结果。可执行语句中不能用中文标点%下面的ZhengFu1To10，ZhengFu11To20，Lamda1To10，T，Taverage，Beta，Omegaf，Omega0，ZhouqiT，YuanPinLv，ZhenFu，XiangWeiCha，XiangWeiChaRad，都是自己定义的“变量”名。Matlab的命令，函数，变量只支持英文字符，不支持中文字符及希腊字符ZhengFu1To10=[12345678910]%输入第1至第10个振幅数据组：θ1--θ10ZhengFu11To20=[11121314151617181920]%输入第11至第20个振幅数据组：θ11--θ20Lamda1To10=log(ZhengFu1To10./ZhengFu11To20)%一次性计算出10个Λ=ln(θi/θi+10)；Matlab中log表示自然对数。./表示涉及数组的除法LamdaAverage=(sum(Lamda1To10))/10%求Λ1至Λ10的平均值Λ_=LamdaAverageT=[1234567891011121314151617181920]%输入自己记录的20个周期数据组TAverage=mean(T)%计算周期的平均值T_=TAverageBeta=LamdaAverage/(10*TAverage)%计算阻尼系数β=BetaOmegaf=2*pi/(TAverage)%计算圆频率ω_f=OmegafOmega0=sqrt(Omegaf^2+Beta^2)%计算系统的固有频率ω_0=Omega0%受迫振动——作幅频特性曲线，相频特性曲线ZhouqiT=[1.5991.5651.5381.5221.5191.5151.5061.4971.4871.471]%输入周期数据组YuanPinLv=2*pi./(ZhouqiT)%计算ω即YuanPinLvZhenFu=[48729810811110694786044]%输入振幅数组据即ZhenFuXiangWeiCha=[233960828995116132142152]%输入相位差(度数)即XiangWeiChaXiangWeiChaRad=XiangWeiCha*pi./180%把相位差转化弧度plot(YuanPinLv,ZhenFu),xlabel('ω—'),ylabel('振幅A—'),title('幅-频特性')%以ω为横轴，振幅为纵轴作图——幅度频特性曲线。xlabel和ylabel定义坐标轴标记plot(YuanPinLv,XiangWeiCha),xlabel('ω—'),ylabel('相位差—'),title('相-频特性')%以ω为横轴，相位差为纵轴作图——相频特性曲线。xlabel和ylabel定义坐标轴标记