主页 / coding / 问题 / 78427505
Accepted
Jonathan Frutschy
Asked: 2024-05-04 09:16:01 +0800 CST

实现非线性函数的可变步长 MATLAB

  • 772

我有以下函数，它返回作为时间函数的非线性信号：

function y = NARMA_optimized(u,mu,delta,phi,psi,m,N,k)
y = zeros(N,1);
% y(0)=0, and also y(-1)=0 etc.
t = 1;
y(t) = mu.* u(t-m-k).* u(t-k) + delta;
for t = 2:N
    v = 1;
    for i = 1:m
        if t-i > 0  % else we add 0
            v = v + y(t-i);
        end
    end
    v = phi.*y(t-1) + v.*psi.*y(t-1) + mu.*u(t-m-k).*u(t-k) + delta;
    y(t) = v;
end  
end

要调用该函数，我执行以下操作：

% NARMA Model Comparison
alpha = 0.2;
beta = 0.15;
f_0_i = 1; 
delta = 0.1;
Mu = 1.5;
m = 10;
phi = 0.3;
psi = 0.05;
k_shift = 0;
Sigma = 1;
M = 100; % Upper simulation time limit [seconds]
f = @(t) Sigma.*f_0_i.*sin(alpha.*t).*cos(beta.*t);
dt = 10^-2; % timescale step size [seconds]
t_span = [0:dt:M];
T = NARMA_optimized(f,Mu,delta,phi,psi,m,length(t_span),k_shift);
figure; 
plot(T);

我想更改当前函数的时间刻度以匹配我表示为的时间的离散表示t_span。

我知道这需要对函数进行两个关键更改：

  1. 开始于t=0, 不t=1.
  2. 将函数输出的更新增加y相当于 的步长dt，而不是1秒。我不确定当前的 for 循环结构是否可以修改以适应这些必要的更改，或者我是否需要从头开始。我知道 MATLAB 从索引 1 开始，而不是索引 0。循环中从索引 0 开始的解决方法如下：
dt = 10^-2; % [seconds]
M = 100; % [seconds]
t_span = [0:dt:M];
for t=1:numel(t_span)
y(t) = [insert value here]
end

我不确定这是否是他们dt在新函数中从索引 0 开始计算的最佳方式。

function

1 个回答

  1. Best Answer

    您可以t_span直接传入该函数。如果您愿意，实时时间将用于u、 或f。然后定义一个indexforfor循环；这index将用于y.

    function y = NARMA_optimized(u,mu,delta,phi,psi,m,t_span,k)
y = zeros(length(t_span),1);
% y(0)=0, and also y(-1)=0 etc.
%t = 1;
t = t_span(1);
index = 1;
y(index) = mu.* u(t-m-k).* u(t-k) + delta;
for t = t_span(2:end)
    index = index + 1;
    v = 1;
    for i = 1:m
        if index-i > 0  % else we add 0
            v = v + y(index-i);
        end
    end
    v = phi.*y(index-1) + v.*psi.*y(index-1) + mu.*u(t-m-k).*u(t-k) + delta;
    y(index) = v;
end  
end

    在主代码中更改以下内容：

    M = 1; % Upper simulation time limit [seconds]
T = NARMA_optimized(f,Mu,delta,phi,psi,m,t_span,k_shift);
plot(t_span,T);

    这是输出： 在此输入图像描述

    • 1

相关问题