Simulink Embedded MATLAB Function

Question

I have tried to introduce the next function in MATLAB embedded function:

function [R2, A2, Fi2, sigma2,fk2,w2] = proof(N, fhz, t, K,tao, sE)
%#codegen

[ll,kk]=size(sE);

A1=zeros(ll,kk);
R1=zeros(ll,kk);

fk1=zeros(ll,kk);
Fi1=zeros(ll,kk);
sigma1=zeros(ll,kk);

g=0;c=1;
sE1=zeros(ll,kk);
     for s=N+1:length(sE)
        d=1;
       
        for r=c:N+g
            
            sE1(d)=sE(r);
            d=d+1;
        end
        c=c+1;

        [R, sigma, w, fk]=pronyF2(sE1,tao);
        
	R1 = R(1)*exp(-1i*2*pi*fhz*t(s));               % Guardar vector de 							cada fasor rotado
        A1(s,:) = abs(R1);                             % Guardar vector de 							Amplitud
        Fi1(s,:) = angle(R1);                          % Guardar vector cada 							¡ngulo
        sigma1(s)=sigma;                             % Guardar vector de cada 							Amortiguamiento
        fk1(s,:)=fk(1);                                   % Guardar vector de 							cada frecuencia
        w1(s,:)=w(1);
        
        g=g+1;
     end

R2=R1;
A2=A1;
Fi2=Fi1;
sigma2=sigma1;
fk2=fk1;
w2=w1;


        
%C·lculo de los coeficientes del Modelo de PredicciÛn Lineal (LPM)
______________________________________________________________________________

function [R,sigma,w,fk] =pronyF2(sE1, tao)
[ll,kk]=size(sE1');
K=1;
n=2*K;
Y=zeros(ll,kk);
q=zeros(ll);
N=5;
        for m=1:N-2
            g=1;
            K=1;
            o=1;
            
            for o=1:2*K 
                Y(m,o)= sE1(n+1-g);
                g=g+1;
            end
            n=n+1;
        end
        K=1;
        n=2*K;
        h=1; 
        for k=n+1:N
            q(h)=sE1(k);
            h=h+1;
        end
       a1=pinv(Y)*q';
% Segundo paso: C·lculo de RaÌces del polinomio caracterÌstico
        P=[1 -a1'];
        z=roots(P);
 %% Tercer paso: C·lculo de Magnitudes y ·ngulos de la seÒal estimada
        Z=complex(zeros(ll,kk));
        for u=1:N
             for v=1:length(z)
                Z(u,v)= z(v,1)^(u-1);
             end  
        end
     
R= 2*pinv(Z)*(sE1');
sigma=(real(log(z)))/tao;
w=angle(z)/tao;
fk=(w/(2*pi));

        %% Tercer paso: C·lculo de Magnitudes y ·ngulos de la seÒal estimada
        Z=complex(zeros(ll,kk));
        for u=1:N
             for v=1:length(z)
                Z(u,v)= z(v,1)^(u-1);
             end  
        end
     
R= 2*pinv(Z)*(sE1');
sigma=(real(log(z)))/tao;
w=angle(z)/tao;
fk=(w/(2*pi));
______________________________________________________________________________
When i press to Run the simulink stop and begin the debugger. This MATLAB function must take whatever signal and perform an estimation in the function pronyF2, then, give these results to R, sigma,w and fk. With this function, we use  the function proof to calculate [R2, A2, Fi2, sigma2,fk2,w2]. I don't know what is the problem?, i'm beginner in simulink using matlab function.

Thank you for your help!!

Simulink Embedded MATLAB Function

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（0 个）

类别

产品

标签

Community Treasure Hunt

Simulink Embedded MATLAB Function

1 个评论 显示 -1更早的评论 隐藏 -1更早的评论

回答（0 个）

类别

产品

标签

另请参阅

Community Treasure Hunt

1 个评论
显示 -1更早的评论隐藏 -1更早的评论