MATLAB Function 次元について

10 次查看（过去 30 天）

显示更早的评论

s 2020-1-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/499058-matlab-function

评论： stozaki 2020-1-25

Simulink内でMATLAB Function ブロック使用して演算しています。

画像のようにMATLAB Functionに入力すると

エラー内容

インデックスが行列の次元を超えています。配列 z は空であるため有効なインデックスを持ちません。

のようにエラーが発生しシミュレーションが実行できません。

次に、コメントアウトされている方(薄く表示されているA)の数値を使用すると、シミュレーションが実行されます。

この違いの原因が分かりません。どちらも出力されるデータ型はdouble型でした。

9 个评论
显示 7更早的评论隐藏 7更早的评论

s 2020-1-7

在 MATLAB Online 中打开

コードはこのようになっています。

どなたか解決策をお願いいたします。

function [aaa,aa,a,a_] = fcn(u, g, xx, time, dd, veppd)
amx=g*u;
      
    xppd = -20:0.01:0;                     
    yyt = 1:0.25:4;                         
    nny =length(yyt);                     
    
        dt_tppd=zeros(nny,1);
        dt_tcr=zeros(nny,1);
        dt_xppd=zeros(nny,1);
        dt_yppd=zeros(nny,1);
        dt_vecr=zeros(nny,1);
        dt_yt=zeros(nny,1);
        
            yppd_z1=zeros(1,length(xppd));
            t_ppd=zeros(1,length(xppd));
            vecr=zeros(1,length(xppd));
            t_cr=zeros(1,length(xppd));          
        
    for k=1:nny                           
        y2 = yyt(1,k);                      
            
        
        yppd_z1 = -y2*xppd./(xx-xppd)-dd/2;
        t_ppd= yppd_z1/veppd;
        vecr = amx*(-time+(time.^2+(-4*xppd)/(2*amx)).^0.5);
        t_cr = -xppd./vecr;                
        
        z = find(t_ppd>t_cr);               
       
        dt_tppd(k,1) = t_ppd(1,z(1,1));
        dt_tcr(k,1) = t_cr(1,z(1,1));
        dt_xppd(k,1) = xppd(1,z(1,1));
        dt_yppd(k,1) = yppd_z1(1,z(1,1));
        dt_vecr(k,1) = vecr(1,z(1,1))*3.6;      
        dt_yt(k,1) = y2-dd/2;                     
    end
    
    hn=2;
    pfuvecr = polyfit(dt_yt,dt_vecr,hn); 
    qu_vecr = polyval(pfuvecr,dt_yt);
    hn=3;
    pfuxppd = polyfit(dt_yt,dt_xppd,hn);
    qu_xppd = polyval(pfuxppd,dt_yt);
    
    aaa = pfuxppd(1);
    aa = pfuxppd(2);
    a = pfuxppd(3);
    a_ = pfuxppd(4);
 
y = u;