Saving variables in a function

Question

Paul Rogers 2021-2-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/745742-saving-variables-in-a-function

评论： Walter Roberson 2021-2-15

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

Hi everyone,

I am calling this function where I get in output the time t and the matrix y:

function [ dy ] = greitzer( t,y,P)
 
%% Energy Losses
deltahfi = P.kf.*(y(1).^2);    %Impeller Friction losses
deltahfd = P.kfd.*(y(1).^2);   %Diffuser Friction losses
%% Throttle valve's parameters
%% con lo stesso Delta da una parte e dall'altra
deltaideal = P.sigma*(P.U2^2); %Deltah0c,ideal
    %% Lossees    
    deltahii = 0.5.*((P.U1-...
        ((cot(P.beta1b).*y(1))/(P.ro1*P.Ai))).^2); %Impeller Incidence Loss
    deltahid = 0.5.*(((P.sigma*P.D2*P.U1)/P.D1)-...
        ((y(1).*cot(P.alfa2b))/(P.ro1*P.Ad))).^2; %Diffuser Incidence Loss
    deltaloss = deltahii+...
        deltahid+deltahfi+deltahfd;
%     deltah0c = deltaideal-...
%         deltaloss-deltan; %deltah0c 
    
    %% Efficiency
    eta = deltaideal./(deltaideal+deltaloss)-P.deltan;
    %% Pressure Rise in the Compressor
    psi_c =((1+((eta.*deltaideal)/(P.T01*P.cp))).^(P.k/(P.k-1)));
%% Greitzer's System
gamma_T = @(t) P.Amplitude*sin(P.w*t)+P.b; % valve's function
dy = [ P.B*(psi_c-y(2)); 
    (1/P.B.*(y(1)-gamma_T(t).*(y(2).^0.5))) ];
end 