Matlab - Graficar con datos generados en una GUI

function uibuttongroup6_SelectionChangeFcn(hObject, eventdata, handles)

if(strcmp(get(eventdata.NewValue,'Tag'),'radiobutton7'))

    set(handles.edit7,'Enable','off')

    set(handles.edit8,'Enable','off')

    z=1;

else

    set(handles.edit7,'String',num2str(0));

    set(handles.edit8,'String',num2str(0));

    set(handles.edit7,'Enable','on');

    set(handles.edit8,'Enable','on');

    z=0;

end

.

.

.

% --- Executes on button press in pushbutton1.

function pushbutton1_Callback(hObject, eventdata, handles)

% hObject    handle to pushbutton1 (see GCBO)

% eventdata  reserved - to be defined in a future version of MATLAB

% handles    structure with handles and user data (see GUIDATA)

%Par?metros fijos

To=293.15;  %Temperatura del aire estandar en grados Kelvin

f = xlsread('parametros','Hoja1','A2:A27');

CaP = xlsread('parametros','Hoja1','B2:B27');

CbP = xlsread('parametros','Hoja1','C2:C27');

CaC = xlsread('parametros','Hoja1','D2:D27');

CbC = xlsread('parametros','Hoja1','E2:E27');

Cero = xlsread('parametros','Hoja1','F2:F27');

Uno = xlsread('parametros','Hoja1','G2:G27');

dBA = xlsread('parametros','Hoja1','H2:H27');

longitud= str2double(get(handles.edit1,'String'));

velocidad = str2double(get(handles.edit2,'String'));

distancia = str2double(get(handles.edit3,'String'));

hm = str2double(get(handles.edit4,'String'));

temperatura = str2double(get(handles.edit5,'String'));

humedad = str2double(get(handles.edit6,'String'));

barrera1 = str2double(get(handles.edit7,'String'));

barrera2 = str2double(get(handles.edit8,'String'));

%Calculo para nivel de potencia sonora por roce aerodinamico

if velocidad>=0 && velocidad<=98

    La=Cero;

    xlswrite('parametros.xls', La, 'Hoja2','B2:B27');

else

    if velocidad>=99 && velocidad<=200

        aero=96*log10(velocidad/200)+67;

        La=aero*Uno;

        xlswrite('parametros.xls', La, 'Hoja2','B2:B27');

    else

        if velocidad>=201 && velocidad<=406

            aero=64*log10(velocidad/100)+60;

            La=aero*Uno;

            xlswrite('parametros.xls', La, 'Hoja2','B2:B27');

        end

    end

end

%Calculo para nivel de potencia sonora por contacto rueda riel 

x = get(handles.radiobutton1,'Value');

if x==1

    Lw1m=CaP*log10(velocidad/100)+CbP;

    Lw=Lw1m+10*log10(longitud);

    xlswrite('parametros.xls', Lw, 'Hoja2','A2:A27');

else

    Lw1m=CaC*log10(velocidad/100)+CbC;

    Lw=Lw1m+10*log10(longitud);

    xlswrite('parametros.xls', Lw, 'Hoja2','A2:A27');

end

%Calculo de atenuacion por presencia de vegetacion

y = get(handles.radiobutton3,'Value');

if y==1

    Avegetacion=Cero;

    xlswrite('parametros.xls', Avegetacion, 'Hoja2','C2:C27');

else

    y = get(handles.radiobutton4,'Value');

    if y==1

        Avegetacion=(0.18*log10(f)-0.31)*distancia;

        xlswrite('parametros.xls', Avegetacion, 'Hoja2','C2:C27');

    else

        Avegetacion=0.01*f.^(-1/3)*distancia;

        xlswrite('parametros.xls', Avegetacion, 'Hoja2','C2:C27');

    end

end

%Calculo de atenuacion por barrera

w = get(handles.radiobutton6,'Value');

if w==1

    N=(2*(barrera1+barrera2-distancia))./f;

    Abarrera=20*log10((sqrt((2*pi)./N))./(tanh(sqrt((2*pi)./N))))+5;

    xlswrite('parametros.xls', Abarrera, 'Hoja2','D2:D27');

else

    Abarrera=Cero;

    xlswrite('parametros.xls', Abarrera, 'Hoja2','D2:D27');

end

%Calculo de atenuacion por suelo

suelo=4.8-((2*hm)/distancia)*(17+300/distancia);

Asuelo=suelo*Uno;

xlswrite('parametros.xls', Asuelo, 'Hoja2','E2:E27');

%Calculo de atenuacion por edificios

Aedificio=Cero;

xlswrite('parametros.xls', Aedificio, 'Hoja2','F2:F27');

%Calculo de propagacion sonora debido a la divergencia geometrica

Adiv=10*log10(4*pi*distancia)*Uno;

xlswrite('parametros.xls', Adiv, 'Hoja2','G2:G27');

%Calculo de indice de directividad

Vson=350; %Velocidad del sonido en m/s a 20 grados centigrados

DI=10*log10((2*longitud*f)/Vson);

xlswrite('parametros.xls', DI, 'Hoja2','H2:H27');

%Calculo de coeficiente de atenuacion debido a la absorcion atmosferica

Fro=24+4.04*10^4*humedad*((0.02+humedad)/(0.391+humedad));

Frn=(temperatura/To)^(1/2)*(9+280*humedad*exp(-4.17*(temperatura/To)^(-1/3)-1));

alfa=869*f.^2.*(1.84*10^(-11)*(temperatura/To)^(1/2)+(temperatura/To)^(-5/2)*(0.01275*(exp(-2239.1/temperatura)./(Fro+(f.^2/Fro)))+0.1068*(exp(-3352/temperatura)./(Frn+(f.^2/Frn)))));

Aabs=alfa*distancia/100;

xlswrite('parametros.xls', Aabs, 'Hoja2','I2:I27');

%Nivel de presion sonora en dB 

Lp=10*log10(10.^(Lw/10)-10.^(Adiv/10)+10.^(DI/10)-10.^(11/10)-10.^(Avegetacion/10)-10.^(Asuelo/10)-10.^(Aedificio/10)-10.^(Aabs/10)-10.^(Abarrera/10));

xlswrite('parametros.xls', Lp, 'Hoja2','J2:J27');

%Nivel de presion sonora en dB con La

Lpaero=10*log10(10.^(Lw/10)+10.^(La/10)-10.^(Adiv/10)+10.^(DI/10)-10.^(11/10)-10.^(Avegetacion/10)-10.^(Abarrera/10)-10.^(Asuelo/10)-10.^(Aedificio/10)-10.^(Aabs/10));

xlswrite('parametros.xls', Lpaero, 'Hoja2','K2:K27');

%Nivel de presi?n sonora en dBA

Lpa=Lp+dBA;

xlswrite('parametros.xls', Lpa, 'Hoja2','L2:L27');

%Nivel de presion sonora en dBA con Lpaero

Lpaaero=Lpaero+dBA;

xlswrite('parametros.xls', Lpaaero, 'Hoja2','M2:M27');

%Nivel de presion son;ora ponderada en pascales

Lpp=2*10.^(-5+Lpa/20);

xlswrite('parametros.xls', Lpp, 'Hoja2','N2:N27');

%Nivel de presion sonora sin ponderar en pascales 

Lpb=2*10.^(-5+Lp/20);

xlswrite('parametros.xls', Lpb, 'Hoja2','O2:O27');

%Calculo de nivel de presion sonora general en decibelios lineales

Lpg1=10.^(Lp/10);

Lpg2=sum(Lpg1);

Lpg=10*log10(Lpg2);

%Calculo de nivel de presion sonora general en decibelios lineales

Lpga1=10.^(Lpa/10);

Lpga2=sum(Lpga1);

Lpga=10*log10(Lpga2);

winopen( 'parametros.xls' );

Y=(Lw);

plot1(Y,f);