fdfd.m

Voltage Distribution of a Triangle (FEM & FDFD Method)

%%%%%%%%%%%%%%%%%%%%%%% In The Name of God %%%%%%%%%%%%%%%%%%%%%%
% 
% Numerical Method in Electromagnetic
% Dr. M.S Abrishamian
% 
% Mohsen Bahrami Panah        8600374                14/2/87
% 
% Problem 20 -- *** Voltage Distribution of a Triangle (FDFD) ***

% clear all
% close all
% clc
% clf
function [Q] = FDFD(g)
k = 0;
% g = 10;
f = g;
for i = 1:g
    for j = 1:f
        x(i,j) = j-1;
        y(i,j) = k;
%         plot(x(i,j),y(i,j),'*b')
%         getframe;
    end
    k = k+1;
    f = f-1;
end
v(1,1:g-1) = 0;
v(2:g-1,1) = 0;
v(g,1) = 50;
v(1,g) = 50;
a = g;
for d = 2:(g-1)
    v(a-1,d) = 100;
    a = a-1;
end
N = 0;
for L = 3:g-1;
    N = N+g-L;
end
hh = waitbar(0,'Please wait...');
for w = 1:N
    waitbar(w/N,hh)
    A(w,w) = 4;
end
k = -1;
for m = (g-5):-1:0
    k = k+1;
    kk = k;
    for n = m:-1:0
        A(1+n+kk,2+n+kk) = -1;
        A(2+n+kk,1+n+kk) = -1;
        k = k+1;
    end
end
k = -1;
p = g+1;
for mm = (g-5):-1:0
    p = p-1;
    k = k+1;
    kk = k;
    for nn = mm:-1:0
        A(1+nn+kk,2+nn+kk+p-4) = -1;
        A(2+nn+kk+p-4,1+nn+kk) = -1;
        k = k+1;
    end
end
close(hh)
k = -1;
B = 200*ones(N,1);
for q = g-5:-1:0
    k = k+1;
    kk = k;
    for s = q:-1:0
        B(1+s+kk,1) = 0;
        k = k+1;
    end
end
V = A\B;
dd = 0;
aa = g+1;
for r = 2:(g-2)
    aa = aa-1;
    for rr = 2:(aa-2)
        dd = dd+1;
        v(r,rr) = V(dd);
    end
end
BB=1;
AA=1;
Delta_X=AA/(g-1);
Delta_Y=BB/(g-1);
a = 0;
for i = 1:g
    vv(i,1:g)=v(g+a,1:g);
    a = a-1;
end
bb = 0;
cc = 0;
for I = g:-1:1
    for J = 1:g-cc
    bb = bb+1;
    Q(bb) = vv(I,J);
    end
    cc = cc+1;
end

% figure(1)
% set(gcf,'Color',[1,0.9,0.9])
% surf(0:Delta_X:AA,0:Delta_X:AA,v)
% colormap(hsv(100))
% title('\bf\it Voltage Distribution of a Triangular with Conductive Sides (Numerical Methode)',...
%     'color','b')
% grid
% % print -djpeg -r300 Voltage_Distribution_Triangular_1.jpg

% figure(1)
% set(gcf,'Color',[1,0.9,0.9])
% surf(0:Delta_X:AA,0:Delta_X:AA,vv')
% colormap(hsv(100))
% title('\bf\it Voltage Distribution of a Triangular with Conductive Sides (Numerical Methode)',...
%     'color','b')
% grid
% print -djpeg -r300 Voltage_Distribution_Triangular_2.jpg

% figure(3)
% set(gcf,'Color',[1,0.9,0.9])
% contour(0:Delta_X:AA,0:Delta_X:AA,v)
% colormap(hsv(100))
% title('\bf\it Voltage Distribution of a Triangular with Conductive Sides (Numerical Methode)',...
%     'color','b')
% grid
% % print -djpeg -r300 Voltage_Distribution_Triangular_3.jpg
% pause(3)
% % close all