cocoa
2008-03-11, 17:02
上周刚接触MATLAB，有个求矩阵特征值的问题怎么都搞不定，哪位高手帮帮忙，先谢啦
是这样的
K(n)=sinc((n*pi)/c+(1-1/c)*b*sinc(b*pi*n)
Q=(2*pi*n/a+k)^2*K(n-m)
其中a,b,c为已知的常数，n m为整数，均从-4取到4，k为未知数,想求矩阵Q的特征值Y关于k的图像
这是用平面波展开法计算一维光子晶体的禁带，就是两种材料周期性放置，其中a为一个周期的长度,c为介电常数，b为两种材料宽度的比值；
我现修改程序如下，特征值是算出来了，可是不会画图了。。。
syms k
n=-4:4;
m=-4:4;
[nn,mm]=meshgrid(n,m);
a=1;
b=0.8;
c=12.25;
A=sinc((nn-mm)*pi)/c+(1-1/c)*b*sinc((nn-mm)*pi*b);
B=2*pi*nn/a+k;
C=(B.*B).*A;
Y=eig(C)
谁能帮帮我~做不出来要挨骂了:lovely:
是这样的
K(n)=sinc((n*pi)/c+(1-1/c)*b*sinc(b*pi*n)
Q=(2*pi*n/a+k)^2*K(n-m)
其中a,b,c为已知的常数，n m为整数，均从-4取到4，k为未知数,想求矩阵Q的特征值Y关于k的图像
这是用平面波展开法计算一维光子晶体的禁带，就是两种材料周期性放置，其中a为一个周期的长度,c为介电常数，b为两种材料宽度的比值；
我现修改程序如下，特征值是算出来了，可是不会画图了。。。
syms k
n=-4:4;
m=-4:4;
[nn,mm]=meshgrid(n,m);
a=1;
b=0.8;
c=12.25;
A=sinc((nn-mm)*pi)/c+(1-1/c)*b*sinc((nn-mm)*pi*b);
B=2*pi*nn/a+k;
C=(B.*B).*A;
Y=eig(C)
谁能帮帮我~做不出来要挨骂了:lovely: