The formula is in the picture.
MATLAB: Could you help me to write this formula in matlap ,and to get it’s code ?
MATLAB
Related Solutions
It would be something like (provide your own functions for f, and set up the values of M and N),
clc;clear;close; % set up, clean up..
% start of Pij part: Pij = (f(i,j))/(∑_(i=1)^M∑_(j=1)^N〖f(i,j)〗)
M=...; N=...; % give valuesfunction f=functionF(x,y)% ... how to find f for x (x is imported i) and y (it is j )
end % function
Pij_single=zeros(M,N); P_of_J=zeros(1,M); % initialize P's
for i=1:M for j=1:N Pij_single(i,j)=functionF(i,j); end P_of_J(i)=sum(Pij_single(i,:));end Pij_sumOnly=sum(P_of_J); % it must be a single value
for i=1:M for j=1:M P(i,j)=functionF(i,j)./Pij_sumOnly; end end% End of Pij part
% Start of H part: H = ─∑_(i=1)^M∑_(j=1)^NPijlog2Pij
clear P_single, P_of_J; % in case you need to save memory
H_of_J=zeros(1,M);Hij=zeros(N,M);% Mark 1
% Note: if you need H = ─∑_(i=1)^M∑_(j=1)^NPijlog2Pij use the code in mark2.
% If you mean H = ─∑_(i=1)^(M-1)∑_(j=1)^(N-1)Pij*log2(Pij) use below:
H=entropy(Pij);disp('H is: '); H% Mark 2. Delete the lines below if not needed
for i=1:M; for j=1:N Hij(i,j)=P(i,j).*log2(P(i,j)); end H_of_J(i)=sum(Hij(:,1)); end H=-sum(H_of_J); disp('H is: '); H % End of the code. Double check initial indices of M, N please.
close allclearvars%x=input('Enter x: ')
x=sin(2*pi*0.1.*(1:1:11)); %h=input('Enter h: ')
h=[1 2 3 4 5 3 1 -1];% convolution
m=length(x);n=length(h);X=[x,zeros(1,n)];H=[h,zeros(1,m)];for i=1:n+m-1 Y(i)=0; for j=1:m if(i-j+1>0) Y(i)=Y(i)+X(j)*H(i-j+1); else end endend% plot results
figure;subplot(3,1,1); stem(x, '-b^'); xlabel('n');ylabel('x[n]'); grid on;subplot(3,1,2); stem(h, '-ms');xlabel('n'); ylabel('h[n]'); grid on;subplot(3,1,3); stem(Y, '-ro');ylabel('Y[n]'); xlabel('----->n'); grid on;title('Convolution of Two Signals without conv function');
Best Answer