Breaking News: Grepper is joining You.com. Read the official announcement!

octave simulator

Itchy Iguana answered on May 14, 2024 Popularity 1/10 Helpfulness 1/10

answer octave simulator

octave simulator

Comment

function [y,n] = sigadd(n1, x1, n2, x2)
% This implements y(n) = x(n) + h(n).
% In order to add two signals, we need to
% make their length equal by padding zeros
% otherwise dimension mismatch will occur.
%n1 = index vector of signal x1
%n2 = index vector of signal x2
%x1 = signal 1
%x2 = signal 2
n = min(min(n1),min(n2)):max(max(n1),max(max(n2))); % Defining the index
%vector for the output signal. For example if n1 = −3:2 and n2 = 0:4
% then n will run from −3 to 4. So we need to define 'n' from the minimum
% of individual minimum of n1 and n2. The same goes for the maximum.
y1 = zeros(1, length(n)); %initializing two vectors to store
y2 = y1; % the zero padded signals
y1(find( (n>=min(n1)) & (n<=max(n1))==1) )= x1; % this places x1 within its
%range/indexes inside the corresponding part of n while the rest stay zeros
%(zero padding).
y2(find( (n>=min(n2)) & (n<=max(n2))==1) )=x2; % same way as x1
% Alternative way as follows
%y1(((n>=min(n1)) & (n<=max(n1)))) = x1;
%y2(((n>=min(n2)) & (n<=max(n2)))) = x2;
y = y1+y2; % the sum of the two signals
end
Now	we	can	verify	the	distributive	property	as	follows.
% verifying distributive law
% x*(h1+h2) = x*h1+x*h2
nx=−2:0;%indexforx
x = [1,2,3]; % signal x
nh1 = −3:2; % index for h1
h1 = nh1.^2; % signal h1 = n^2
nh2 = 0:3; % index for h2
h2 = nh2+2; % signal h2 = n+2
% first calculate y = h1+h2 using sigadd function
[y, ny] = sigadd(nh1, h1, nh2, h2);
%now calculate z= x*y = x*(h1+h2)
 [z,k] = convolution_func(nx, x, ny, y);
 % next, calculate x*h1+x*h2
 %p1=x*h1
 [p1, np1] =convolution_func(nx, x, nh1, h1);
 %p2=x*h2
 [p2, np2] =convolution_func(nx, x, nh2, h2);
 % finally p = p1+p2 = x*h1+x*h2
 [p, np] = sigadd(np1, p1, np2, p2);
 %displaying the results to show z = p
 subplot(211)
 stem(k,z,'Linewidth', 2);
 xlabel('k');ylabel('z[k]'); title('z= x*y = x*(h1+h2)')
 subplot(212)
 stem(np, p,'Linewidth', 2);
 xlabel('np');ylabel('p[np]'); title('p = p1+p2 = x*h1+x*h2')

Popularity 1/10 Helpfulness 1/10 Language whatever

Source: octave.org

Tags: whatever

Link to this answer
Share Copy Link

Contributed on May 14 2024

Itchy Iguana

0 Answers Avg Quality 2/10

octave simulator

Contents

More Related Answers

octave simulator

Grepper

Documentation

Social

Legal

Contact

Oops, You will need to install Grepper and log-in to perform this action.