Non-linear simulink function block help.

Question

Mark Dawson 2019-9-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/480317-non-linear-simulink-function-block-help

回答： Priyanshu Mishra 2019-9-16

So I have my actual matlab file working perfectly, however, I need to transfer this to a simulink file that takes input x(t) and returns output phi(t). Just in standard input, process, output format. Only problem is I've never used simulink, and ive got no Idea where to start with this. The code I used is shown bellow:

    %System Constants
    L1 = 20;
    L1A = 10;
    L1B = 10;
    L2 = 50;
    H= 0.47;
    %System Initial Physical Parameters
    lambda1p = 2;
    lambda2p = 2;
    lambda3p = 1/(lambda1p*lambda2p);
    %System Physical Constants
    M = 17.5;                                 %Mass of the mass (g)
    c = 0.696;                                %Damping Co-efficient
    mu = 562.01;                              %Small Strain shear modulus (kPa)
    epsilon = 2.39e-11;                       %Di-electric Permitivity (F/m)
    roe = 1522;
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%Mathematical Data%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %Input Waveform Charecteristics
    A = 0.011;                                          %Desired Peak Amplitude (mm)
    T = 20;                                             %Waveform Period (S)
    f = 1/T;                                            %Frequency (Hz) of Waveform
    %Symbolic Equations
    syms t
    %x = @(t)A*sin(2*pi*f*t);                                                                   %Sinusoidal Input Trajectory 
    %x = @(t)(2*A/pi)*asin((sin(2*pi*f*t)));                                                    %Triangular Input Trajectory
    x = @(t)((((heaviside(t-1)+heaviside(9-t))/2)+((heaviside(t-3)+heaviside(6-t))/2))-1)*A;    %Stair Input Trajectory
    lambda1A = @(t)((L1A*lambda1p)+x(t))/L1A;                                                   %Stretch Factor 1A eqn
    lambda1B = @(t)((2*lambda2p)-lambda1A(t));                                                  %Stretch Factor 1B eqn
    lambda1Adot=symfun(diff(lambda1A(t)),t);                                                    %1st Differential of Stretch Factor 1A (array)
    lambda1Adot_t=@(t) eval(lambda1Adot);                                                       %1st Differential of Stretch Factor 1A (function)
    lambda1Bdot=symfun(diff(lambda1B(t)),t);                                                    %1st Differential of Stretch Factor 1B (array)
    lambda1Bdot_t=@(t) eval(lambda1Bdot);                                                       %1st Differential of Stretch Factor 1B (function)
    lambda1A2dot=symfun(diff(diff(lambda1Adot(t)),t),t);                                        %2nd Differential of Stretch Factor 1A (array)
    lambda1A2dot_t=@(t) eval(lambda1A2dot);                                                     %2nd Differential of Stretch Factor 1A (function)
    lambda1B2dot=symfun(diff(diff(lambda1Bdot(t)),t),t);                                        %2nd Differential of Stretch Factor 1B (array)
    lambda1B2dot_t=@(t) eval(lambda1B2dot);                                                     %2nd Differential of Stretch Factor 1B (function)
    %Timebase
    t = [0.01:0.1:9.99];
    %Numerical Arrays
    x = x(t);                                                                                   %Store Numerical results of function for given t
    lambda1A = lambda1A(t);                                                                     %Store Numerical results of function for given t
    lambda1B = lambda1B(t);                                                                     %Store Numerical results of function for given t
    lambda1Adot_t=lambda1Adot_t(t);                                                             %Store Numerical results of function for given t
    lambda1Bdot_t=lambda1Bdot_t(t);                                                             %Store Numerical results of function for given t
    lambda1A2dot_t=lambda1A2dot_t(t);                                                           %Store Numerical results of function for given t
    lambda1B2dot_t=lambda1B2dot_t(t);                                                           %Store Numerical results of function for given t
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%Charecteristic Equation%%%%%%%%%%%%%%%%%%%%%%%%%
    Phi = (((L1A*H)./(L2*epsilon*(lambda2p^2)*lambda1A)).*((M.*lambda1A2dot_t)+(c.*lambda1Adot_t)-(((mu*L2*H)./L1A).*((2*lambda1p)-(2*lambda1A)-((((2*lambda1p)-(lambda1A)).^-3)*(lambda2p^-2))+((lambda1A.^-3)*(lambda2p^-2)))))).^0.5;