damping matlab,Matlab绘制阻尼系数曲线

echo off

%plotting frequency responses of sdof model for different

%damping values.Calculates and plots magnitude and phase for a single

%degree of freedom system over a range of damping values.

clf;

clear all;

%assign values for mass, percentage of critical damping, and stiffnesses

%zeta is a vector of damping values from 10% to 100% in steps of 10%

m = 1;

zeta = 0.1:0.1:1;%0.1 = 10% of critical

k = 2*pi*2*pi;

wn = sqrt(k/m);

%Define a vector of frequencies to use, radians/sec.The logspace command uses

%the log10 value as limits, i.e. -1 is 10^-1 = 0.1 rad/sec, and 1 is

%10^1 = 10 rad/sec.The 400 defines 400 frequency points.

w = logspace(0,2,400);

%pre-calculate the radians to degree conversion

rad2deg = 180/pi;

%define s as the imaginary operator times t