MATLAB: How to populate the next term in an If/elseif code

iif else statementif statement

Here's my code:

N=[];

for i=1:length(timeStep);

    if timeStep(i)>33;
      N(i)=1;
    elseif timeStep(i)<=33;
      N(i)=0;
    end
end
n=N';
n(1)=1;

If timeStep at i is greater than 33, I would like N of the NEXT term (i+1) to be populated with a 1. However, changing the code to N(i+1)=1 does not work.

What would make this work for me?

Thank you in advance.

Best Answer

Why should it not work to change it to N(i+1)=1?

I'm wondering what's the desired output for timeStep between 35 and 40?

In any case, you can very simply vectorize the code:

N=[1 timeStep>40];

Related Solutions

MATLAB: How to fit a common linear trend observed across multiple sensors

You can just replicate the x-values and linearize all your data:

%% Make some fake noisy measurements
timeStep = 1:100; % Time step
nSensors=10;
Temperature=zeros(nSensors,numel(timeStep));
for iSensor = 1:nSensors % Loop through sensors
    % Dimensions of Temperature: nSensors x nTime
    Temperature(iSensor,:) = (5 + rand(1,1))*timeStep + ...% Add noise to the true slope of 5
        (rand(1, length(timeStep))-0.5)*100 + 7; % Add noise to the true offset of 7
end
figure(1),clf(1)
plot(timeStep, Temperature); 
xlabel('Time'); ylabel('Temperature'); title('Noisy Temperature');
timeStep2=ones(size(Temperature)).*timeStep;%lazy repmat
p=polyfit(timeStep2(:),Temperature(:),1);
hold on
plot(timeStep,polyval(p,timeStep),'--k')
hold off

MATLAB: How to write matlab code for moments

You may need to convert your input image to double. In uint8 numbers greater than 255 became 255.

function [M]= moments(I)
I=double(I);
[r c]=size(I); 
m=zeros(r,c); 
% geometric moments 
for i=0:1 
    for j=0:1 
        for x=1:r 
            for y=1:c 
                m(i+1,j+1)=m(i+1,j+1)+(x^i*y^j*I(x,y)); 
            end
        end
    end
end
xb=m(2,1)/m(1,1); 
yb=m(1,2)/m(1,1);
% central moments 
u=[ 0 0 0 0;0 0 0 0;0 0 0 0;0 0 0 0]; 
for i=0:3 
    for j=0:3 
        for x=1:r 
            for y=1:c 
                u(i+1,j+1)=u(i+1,j+1)+(x-xb)^i*(y-yb)^j*I(x,y); 
            end
        end
    end
end
% scale invariant moments 
n=[ 0 0 0 0;0 0 0 0;0 0 0 0;0 0 0 0]; 
for i=0:3 
    for j=0:3 
        n(i+1,j+1)=u(i+1,j+1)/(u(1,1)^(1+(i+j)/2)); 
    end
end
%rotation invariant moments 
I_1= n(3,1)+ n(1,3); 
I_2=(n(3,1)- n(1,3) )^2+ (2*n(2,2))^2; 
I_3=(n(4,1)-3*n(2,3))^2+ (3*n(3,2)-n(1,4))^2; 
I_4=(n(4,1)+n(2,3))^2+ (n(3,2)+n(1,4))^2; 
I_5=(n(4,1)-3*n(2,3))*(n(4,1)+n(2,3))*((n(4,1)+n(2,3))^2-3*(n(3,2)+n(1,4))^2)...
    +(3*n(3,2)-n(1,4))*(n(3,2)+n(1,4))*(3*(n(4,1)+n(2,3))^2-(n(3,2)+n(1,4))^2); 
I_6=(n(3,1)-n(1,3))*((n(4,1)+n(2,3))^2-(n(3,2)+n(1,4))^2)+ 4*n(2,2)*(n(4,1)...
    +n(2,3))*(n(3,2)+n(1,4)); 
I_7=(3*n(3,2)-n(1,4))*(n(4,1)+n(2,3))*((n(4,1)+n(2,3))^2- 3*(n(3,2)+n(1,4))^2 )...
    - (n(1,4)-3*n(2,3))*(n(3,2)+n(1,4))*(3*(n(4,1)+n(2,3))^2-(n(3,2)+n(1,4))^2);
M= [I_1 I_2 I_3 I_4 I_5 I_6 I_7];

Best Answer

Related Solutions

MATLAB: How to fit a common linear trend observed across multiple sensors

MATLAB: How to write matlab code for moments

Related Question