谐波平衡 matlab,增量谐波平衡法程序

根据网上那个范德波极限环程序编的，为什么我的幅频图是很奇怪的折线

clear all

close all

clc

tic

%=====输入基本参数(已知条件)===================================

%t=1;

syms t

syms w0

m0=0.2712e6;k0=0.3707e4;c0=0.4439e4;v21=1.6039;

mu1=0.05;mu2=0.005;

m1=m0*mu1;m2=m0*mu2;

xi1=0.1195;xi2=0.2177;%阻尼比

lambda1=0.9266;lambda2=0.8117;%频率比

w00=sqrt(k0/m0);

w1=lambda1*w00;w2=lambda2*w00;

k1=w1^2*m1;

k2=w1^2*m2;

c1=2*xi1*m1*w1;

c2=2*xi2*m2*w2;

%=======================质量矩阵m=======================================

m=diag([m0/v21 m1 m2]);

%=======================阻尼矩阵c=======================================

c=[c0/v21+c1*v21+c2*v21 -c1*v21 -c2*v21;

-c1*v21 c1*v21 0;

-c2*v21 0 c2*v21];

%========================线性刚度矩阵k====================================

k=[k0/v21+k1*v21 -k1*v21 0;

-k1*v21 k1*v21 0;

0 0 0];

%===========