clear;
clc;
syms f1 f2 f3 f4 f5 f6 f7 %定义未知参数频率
%syms m
S=0.6e-3;
G=0.1e-3;
P=0.8e-3;
a=2.0e-3;
D=0.57e-3;
d=0.4e-3;
c=3e+08;
k1=zeros(1,7);%以上是慢波结构几何参数
beta=[1 121 241 361 481 601 721]; %相位常数的定义值
beta=beta*pi/180;
m=0:1:7; %函数里面循环
f=[f1 f2 f3 f4 f5 f6 f7];
aa = [1 zeros(1,6)]; %冲激响应,a(1,1)=1,其他都等于0
H =sym(zeros(7,7)) ;
%H=0;
I=eye(7,7);
for jj=1:7
w(1,jj)=2*pi.*f(1,jj);%角频率w=2*pi*f
k0(1,jj)=w(1,jj)/c;% k0=w/c
for ii=1:7
for nn=1:7
for mm=1:7
%这里把公式拆成各个单元模块,并构建矩阵A,B,C,线性代数中齐次方程组((A-I)*a-B*b=0 Ba-(C+I)b=0)
%再有线性代数中齐次线性方程组有解
%则系数行列式为0 构建E