hey everyone. this is my program. I'm changing b parameter, so the plot changes, but the Shannon's entropy is always the same. the entropy of grayscale image (in rec plot) is always changing.
clear all global w0 alfa f gama teta b del=0.0001;
% w01=1; % w02=0.1; % alfa=0.1; % beta=0.01; % f1=0; % f2=1.2; % gama1=0; % gama2=0.4; % % w01=0.8; % w02=0.9; % alfa=0; % beta=0; % f1=0; % f2=0; % gama1=0; % gama2=0;
w0=0.001; %w0=0.1 alfa=1; %beta=1; f=0.0001; %f=0;
gama=1; a=1; b=0.0000064; % teta=0.005; %teta=0.09; %teta=0.06; %teta=0.009; teta=0.001;
% y10=3;
y10=0; y20=0;
t0=0; tm=50000; %tm=100; % %options=odeset('RelTol',del,'AbsTo',[del del]); %cdomileba % % [t,y]=ode45('rigid26',[0 tm],[y10 y20 y30 y40],options); % [t,y]=ode45('rigid26',[0 tm],[y10 y20 y30 y40],options); %[t,y]=ode45('rigid2b',[0:2*pi/gama1:4000*pi/gama1],[y10 y20],options);
options=odeset('RelTol',del,'AbsTo',[del del]); %cdomileba
%[t,y]=ode45('rigid2b',[0::4000*pi/gama1],[y10 y20 y30 y40],options);
[t,y]=ode45('x',[t0:1:tm],[y10 y20],options);
%y = awgn(z,0.0001,'measured'); %plot(t,[z y]) %plot(t,y(:,1),t,sin(w0*t)); %plot(t,y(:,1));
e = wentropy([t,y],'shannon')
%e=wentropy([t,y],'shannon')
Best Answer