Matematické Fórum

Zdravím, nebyl by tu prosím někdo tak hodný a neupravil by mi tento skript do podoby funkce, ktera je nadepsana? vubec si nevim rady:(

function Body(A,B)

p0 = 2 * 10^-5;
I0 = 1 * 10^-12;
W0 = 1 * 10^-12;
Sn = 20;

%uložení vypočtených hodnot do matice Lp1 a I1
for i=1:2
    for j=1:2
      Lp1(i,j) = 20*(log10(A(i,j)/p0));
      I1(i,j) = 10^(Lp1(i,j)/10)+log10(I0);
    end
end


%uložení vypočtených hodnot do matice Lp a I
for i=1:2
    for j=1:2
      Lp2(i,j) = 20*(log10(B(i,j)/p0));
      I2(i,j) = 10^(Lp2(i,j)/10)+log10(I0);
    end
end


I = I1+I2


%metoda mereni v bodech


for k=1:4
    for l=1:4
      Wb0(k,l) = I(k,l)*Sn;
      Lwb(k,l) = 10*(log10(Wb0(k,l)/W0));
    end
end

F=[200 250 315 400;  500 630 800 1000; 1250 1600 2000 2500; 3150 4000 5000 6300];
subplot(1,2,1);
bar(F(:),Lwb(:),2)
title('Třetinooktávový spektrogram hladiny akustického výkonu při metodě měření v bodech ')
xlabel('F[Hz]')
ylabel('Lws[dB]')

subplot(1,2,2);
image(I);
colormap hsv;
colorbar ;