I'm getting a matrix error of inconsistently on line 45, any tips how I can fix it.

Question

BAILEY MCMASTER 2023-4-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1955624-i-m-getting-a-matrix-error-of-inconsistently-on-line-45-any-tips-how-i-can-fix-it

评论： BAILEY MCMASTER 2023-5-1

在 MATLAB Online 中打开

N=10;

L=10;

D=.01;

dt=.01;

k_on=1;

k_off=.1;

F=1;

tmax=100;

x=L*rand(N,1);

y=L*rand(N,1);

state=zeros(N,1);

figure;

axis([0 L 0 L]);

set(gca,'nextplot','replacechildren');

for t = 0:dt:tmax

dx=sqrt(2*D*dt)*randn(N,1);

dy=sqrt(2*D*dt)*randn(N,1);

x=x+dx;

y=y+dy;

x(x<0)=0;

x(x>L)=L;

y(y<0)=0;

y(y>L)=L;

p_on=k_on*dt*F;

p_off=k_off*dt;

for i = 1:N

if state(i)==0

if rand < p_on

state(i)=1;

end

else

if rand<p_off

state(i)=0

end

colors= repmat([1 0 0], N,1);

colors(state == 1,:,:) = repmat([0 1 0],sum(state),1);

scatter(x,y,50,colors,'filled');

title(sprintf('Time=%.2f',t));

drawnow;

if any(state == 1)

p_bind=1-exp(-k_on*dt*F);

if rand<p_bind

[~,i]=max(state);

state(i)=2;

colors(i,:)=[0]

end

colors = 10×3

1 0 0 1 0 0 1 0 0 1 0 0 0 0 0 1 0 0 1 0 0 1 0 0 0 1 0 0 1 0

Unable to perform assignment because the size of the left side is 2-by-3 and the size of the right side is 4-by-3.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Walter Roberson 2023-5-1

What is the difference between this and https://www.mathworks.com/matlabcentral/answers/1955619-i-keep-getting-this-matrix-error-when-i-run-my-code-any-tips-how-i-can-fix-it-its-line-45?s_tid=srchtitle ?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

LeoAiE 2023-5-1

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1955624-i-m-getting-a-matrix-error-of-inconsistently-on-line-45-any-tips-how-i-can-fix-it#answer_1226369

在 MATLAB Online 中打开

It seems there is an issue with the assignment of the colors variable in the last part of your code.

N=10;
L=10;
D=.01;
dt=.01;
k_on=1;
k_off=.1;
F=1;
tmax=100;
x=L*rand(N,1);
y=L*rand(N,1);
state=zeros(N,1);
figure;
axis([0 L 0 L]);
set(gca,'nextplot','replacechildren');
for t = 0:dt:tmax
    dx=sqrt(2*D*dt)*randn(N,1);
    dy=sqrt(2*D*dt)*randn(N,1);
    x=x+dx;
    y=y+dy;
    x(x<0)=0;
    x(x>L)=L;
    y(y<0)=0;
    y(y>L)=L;
    p_on=k_on*dt*F;
    p_off=k_off*dt;
    for i = 1:N
        if state(i)==0
            if rand < p_on
                state(i)=1;
            end
        else
            if rand<p_off
                state(i)=0
            end
        end
    end
    colors = repmat([1 0 0], N,1);
    colors(state == 1,:,:) = repmat([0 1 0],sum(state == 1),1);
    scatter(x,y,50,colors,'filled');
    title(sprintf('Time=%.2f',t));
    drawnow;
    if any(state == 1)
        p_bind=1-exp(-k_on*dt*F);
        if rand<p_bind
            [~,i]=max(state);
            state(i)=2;
            colors(i,:)=[0 0 1]; % Corrected color assignment
        end
    end
end