fdtd22.m

一个有关FDTD数值色散问题的MATLAB程序

clear all
close all
x=0:.001:2*pi;
x=rad2deg(x);
syms x;
ezpolar(sqrt(((sin(pi*1/2*cos(x))).^2)*((cos(x)).^2/(pi*1/2*cos(x)).^2)+((sin(pi*1/2*sin(x))).^2)*((sin(x)).^2/(pi*1/2*sin(x)).^2)));
hold on;
ezpolar(sqrt(((sin(pi*1/3*cos(x))).^2)*((cos(x)).^2/(pi*1/3*cos(x)).^2)+((sin(pi*1/3*sin(x))).^2)*((sin(x)).^2/(pi*1/3*sin(x)).^2)));
hold on;
ezpolar(sqrt(((sin(pi*1/4*cos(x))).^2)*((cos(x)).^2/(pi*1/4*cos(x)).^2)+((sin(pi*1/4*sin(x))).^2)*((sin(x)).^2/(pi*1/4*sin(x)).^2)));
hold on;
ezpolar(sqrt(((sin(pi*1/8*cos(x))).^2)*((cos(x)).^2/(pi*1/8*cos(x)).^2)+((sin(pi*1/8*sin(x))).^2)*((sin(x)).^2/(pi*1/8*sin(x)).^2)));
hold on;
ezpolar(sqrt(((sin(pi*1/20*cos(x))).^2)*((cos(x)).^2/(pi*1/20*cos(x)).^2)+((sin(pi*1/20*sin(x))).^2)*((sin(x)).^2/(pi*1/20*sin(x)).^2)));
grid