close allclear alln=-2:2x(n~=0)=0x(n==0)=1subplot(3,1,1)stem(n,x)w(n~=2)=0w(n==2)=1subplot(3,1,2)stem(n,w)y=x-nsubplot(3,1,3)stem(n,y)MATLAB: Not able to subtract two discrete signals
not able to subtract two discrete signals
Related Solutions
This will eliminate the error:
for n = 2:length(x)-2 y(n+1)=(-1/3)*x(n+2)+(1/2)*x(n+1)+(1/3)*x(n)+y(n);endhowever you have to change the stem arguments accordingly:
stem(0:1:49, y, 'r-')N = [0,1,2,3,4,5]; % n = ts/tau
Y = zeros(1,numel(N)); % new
C = cell(1,numel(N)); % newf = @(t_star,y,n) 1-y-(t_star/n); % y = theta* , t_star = t*
for kk = 2:numel(N) % start from 2
[t_star,y] = ode45(f,[0 N(kk)],0,0,N(kk)); %solves ode with n values
line(t_star,y,'color',rand(1,3),'Linewidth',2),grid on ylabel('\theta*') xlabel('t*') ylim([0 1]) % new:
Y(kk) = y(end); C{kk} = sprintf('n=%u',N(kk));endlegend(C{:})hold onplot(N,Y,'-ok','Linewidth',2)Giving:
If you want a smoother line, then you will need to curve-fit or interpolate the data, e.g.:
Ni = linspace(N(1),N(end),1000);Yi = interp1(N,Y,Ni,'pchip');plot(Ni,Yi,'-k','Linewidth',2)
Best Answer