mesh_electro.m

% Faculty of Engineering, University of Fukui. 2004/10/10
% coded by Sokthai Chan ([e-mail])
% Electric Distribution     

clear all;

OMEGA=2;
N=67;
L=pi;
t=linspace(-L,L,N);
[X,Y]=meshgrid(t);

% coordinates setting
v=[-L L -L L -20 0];
% sin wave function with the variation by OMEGA
Z1 = abs(sin(OMEGA*X*0.7).*sin(OMEGA*0.8*Y));
% cos wave function with the variation by OMEGA
Z2 = abs(cos(OMEGA*X*0.18).*cos(OMEGA*0.19*Y));

% Z3 : sin and cos wave are multipled
Z3 = max(-20,10*log10(max(0.01,abs(sin(OMEGA*X*0.7).*sin(OMEGA*0.8*Y)))));

% change it to logarithm
Z = max(-20,10*log10(max(0.01,Z1.*Z2)));

% plot in mesh 
mesh(X,Y,Z);
axis(v);