谐波平衡 matlab,请问有会增量谐波平衡法的，用matlab编程

不知道为啥连接传不上去

clear all

clc

tic

%定义各参数

syms t

w0=3;

epsR=0.001;

m=[1 0；0 1];

epsilon=0.24;

r=0.4;

delta=0.56;

k=[1+epsilon*r -epsilon*r;-epsilon*r 1+epsilon*delta];

Cs=[cos(t) cos(3*t) sin(t) sin(3*t)];

Cs1=diff(Cs，t，1);

S=[cos(t) cos(3*t) sin(t) sin(3*t) 0 0 0 0;0 0 0 0 cos(t) cos(3*t) sin(t) sin(3*t)];

A1=[1 1 1 1]';

A2=[1 1 1 1]';

A0=[A1；A2];

T1=[eye(4，4) zeros(4，4)];

T2=[zeros(4，4) eye(4，4)];

S2=diff(S，t，2);

fm=inline(S'*m*S2);

M=quadv(fm，0，2*pi);

fk=inline(S'*k*S);

K=quadv(fk，0，2*pi);

S1=diff(S，t，1);