-1

% Lathi MS2P4 pg 232 の Signals and Systems からわずかに変更されたコード。

このコードを関数にするにはどうすればよいですか? そのため、「呼び出す」だけで使用できます。この全体を何度も貼り付ける代わりに。これは、2 つの関数の畳み込みアニメーションです。

figure(1) % Create figure window and make visible on screen 
x = @(t) cos(2*pi*3*t).* u(t);
h = @(t) u(t)-u(t-1);
dtau = 0.005;
tau = -1:dtau:4;
ti = 0;
tvec = -1:0.1:5;
y = NaN*zeros(1,length(tvec)); % Pre-allocate memory
for t = tvec,
    ti = ti + 1; % Time index
    xh = x(t-tau).*h(tau);
    lxh = length(xh);
    y(ti) = sum(xh.*dtau); % Trapezoidal approximation of integral
    subplot(2,1,1)
    plot(tau,h(tau),'k-',tau,x(t-tau),'k--',t,0,'ok')
    axis([tau(1) tau(end) -2.0 2.5])
    patch([tau(1:end-1);tau(1:end-1);tau(2:end);tau(2:end)],...
        [zeros(1,lxh-1);xh(1:end-1);xh(2:end);zeros(1,lxh-1)],...
        [0.8 0.8 0.8],'edgecolor','none')
    xlabel('\tau')
    legend('h(\tau)','x(t-\tau)','t','h(\tau)x(t-\tau)')
    c = get(gca,'children');
    set(gca,'children',[c(2);c(3);c(4);c(1)]);
    subplot(2,1,2)
    plot(tvec,y,'k',tvec(ti),y(ti),'ok')
    xlabel('t')
    ylabel('y(t) = \int h(\tau)x(t-\tau) d\tau')
    axis([tau(1) tau(end) -1.0 2.0])
    grid
    drawnow
%   pause
end
4

2 に答える 2

2

ここでfunction説明されているように、Matlab の を使用できます。

そのため、一番上に「function yourFunctionName」と書く必要がありyourFunctionName、この場合はその名前を使用して呼び出すことができます。

于 2013-02-10T19:03:12.103 に答える
0

正確なコードで機能したもの、簡単な関数定義、およびu(t)OPから欠落していた追加が次のとおりです。

figure(1) % figure window outside function
u = @(t)heaviside(t); % define step function
x = @(t) cos(2*pi*3*t).* u(t);
h = @(t) u(t)-u(t-1);
dtau = 0.005;

convol_anime(x, h, dtau); % pass function handles + parameters

すでに定義している場所：

function convol_anime(x, h, dtau)

tau = -1:dtau:4; ti = 0; tvec = -1:0.1:5;
y = nan*zeros(1,length(tvec)); 
...
于 2013-02-10T20:53:38.103 に答える