Metody numeryczne lab 2

file='test.ckt'

file =

    'test.ckt'

readckt
loadMatrix
v=G\b

v =

    9.9751
    7.4813
    4.9875
    2.4938

dobór delta x będzie wskazywał liczbę węzłów (przyjąć to delta x) -> delta x = 1/10 -> 11 węzłów
dla prostego przypadku nie trzeba korzystać z wielu elementów aby uzyskać rozkład gładki



%%%%%
% Load the voltage sources.
for i = 1:size(Vsources,1)
    n1 = Vsources(i,2);
    n2 = Vsources(i,3);
    vol = Vsources(i,1);
    if ((n1 == 0) & (n2 ~= 0))
        sourcenodes(n2) = -vol;
    elseif (n1 ~= 0) & (n2 == 0)
        sourcenodes(n1) = vol;
    end
end
%%%%%

wynik = parabola z maksimum w środku z 2 radiatorami


%%%%%
% Load the voltage sources.
for i = 1:size(sourcenodes,2)
    n1 = i;
    vol = sourcenodes(i);

    if (sourcenodes(i) > 0)
        RHS(n1) = vol;
        Matrix(n1,:) = 0;
        Matrix(n1,n1) = 1;
    end
end
%%%%%

%%%%%%%%%%%%%%%%działa%%%%%%%%%%%%%%%%%%%
% Load the voltage sources. (dodane)
for i = 1:size(sourcenodes,2)
    volt = sourcenodes(i);
    if(sourcenodes(i)~=0)
        RHS(i) = volt;
        Matrix(i,:) = 0;
        Matrix(i,i) = 1;
    end
end