0

你好再次合乎逻辑的朋友！我知道这是一个相当复杂的问题，所以请多多包涵！我想我已经设法把它归结为两个细节： - 我需要两个我似乎无法工作的循环......首先；变量rollers(1).ink 是一个包含墨水值的(12x1) 向量。该程序在每个连接处的辊子之间平均分配墨水。我试图让rollers(1).ink 仅在特定时间步与rollers(2) 交互。每整转一次，墨水应转移到系统中，即 nTimesSteps = nBins_max 的每个倍数。当系统旋转时，墨水不应转移回滚轮 (1).ink - 每转一次，它应该只将墨水引入系统一次，而不会带回任何墨水。目前我已经设置了 rollings(1).ink = one 但仅用于测试。我真的被困在这里了！第二; 它需要这样做的原因是因为在模拟结束时我还希望以打印图像的形式去除墨水。图像应该是我系统中最后一个滚筒上墨水的反射，并且该值的一半应该从最后一个滚筒中移除，并在每次旋转时从系统中取出。最后一个滚筒上剩余的墨水应回收并在系统中“重新分离”，为下一次旋转做好准备。所以……我认为它是在循环开始的第 86 行，我需要做所有这些事情。在伪中，对于间歇性的馈送，我一直在尝试类似的东西：最后一个滚筒上剩余的墨水应回收并在系统中“重新分离”，为下一次旋转做好准备。所以……我认为它是在循环开始的第 86 行，我需要做所有这些事情。在伪中，对于间歇性的馈送，我一直在尝试类似的东西：最后一个滚筒上剩余的墨水应回收并在系统中“重新分离”，为下一次旋转做好准备。所以……我认为它是在循环开始的第 86 行，我需要做所有这些事情。在伪中，对于间歇性的馈送，我一直在尝试类似的东西：

For k = 1:nTimeSteps
While nTimesSteps = mod(nTimeSteps, nBins_max) == 0 % This should only output when nTimeSteps is a whole multiple of nBins_max i.e. one full revolution
‘Give me the ink on each segment at each time step in a matrix’
End

averageAmountOfInk 的输出是我想要返回此数据的确切格式，除了我真的不需要平均值，只需要每个时刻的实际值。当我尝试使用以下方法重新创建它时，我不断收到尺寸不匹配的错误：

For m = 1:nTimeSteps
    For n = 1:N
Rollers(m,n) = rollers(n).ink’;
    End
End

如果有人有兴趣了解它目前的功能，我将在下面发布完整的代码。最后还有一个功能，当然需要将其保存到单独的文件中。我已经多次发布了这个问题的变体，但我完全意识到这是一个相当棘手的问题，而且我发现很难通过互联网传达我的意图！如果有人对我缺乏编程技能有任何想法/建议/一般侮辱，请随时回复！

%% Simple roller train
% # Single forme roller
% # Ink film thickness = 1 micron

clc
clear all
clf
% # Initial state
C = [0,70;   % # Roller centres (x, y)
     10,70;
     21,61;
     11,48;
     21,34;
     27,16;
     0,0
     ];
R = [5.6,4.42,9.8,6.65,10.59,8.4,23];    % # Roller radii (r)
% # Direction of rotation (clockwise = -1, anticlockwise = 1)
rotDir = [1,-1,1,-1,1,-1,1]';
 N = numel(R);           % # Amount of rollers

% # Find connected rollers
isconn = @(m, n)(sum(([1, -1] * C([m, n], :)).^2)...
    -sum(R([m, n])).^2 < eps);
[Y, X] = meshgrid(1:N, 1:N);
conn = reshape(arrayfun(isconn, X(:), Y(:)), N, N) - eye(N);

% # Number of bins for biggest roller
nBins_max = 50;
nBins = round(nBins_max*R/max(R))';

% # Initialize roller struct
rollers = struct('position',{}','ink',{}','connections',{}',...
    'rotDirection',{}');

% # Initialise matrices for roller properties
for ii = 1:N
    rollers(ii).ink = zeros(1,nBins(ii));
    rollers(ii).rotDirection = rotDir(ii);
    rollers(ii).connections = zeros(1,nBins(ii));
    rollers(ii).position = 1:nBins(ii);
end

for ii = 1:N
    for jj = 1:N
        if(ii~=jj)
            if(conn(ii,jj) == 1)
               connInd = getConnectionIndex(C,ii,jj,nBins(ii));
               rollers(ii).connections(connInd) = jj;
            end
        end
    end
end

% # Initialize averageAmountOfInk and calculate initial distribution
nTimeSteps = 1*nBins_max;
averageAmountOfInk = zeros(nTimeSteps,N);
inkPerSeg = zeros(nTimeSteps,N);
for ii = 1:N
    averageAmountOfInk(1,ii) = mean(rollers(ii).ink);
end

% # Iterate through timesteps
for tt = 1:nTimeSteps
        rollers(1).ink = ones(1,nBins(1));

        % # Rotate all rollers
    for ii = 1:N
                rollers(ii).ink(:) = ...
                circshift(rollers(ii).ink(:),rollers(ii).rotDirection);
    end

% # Update all roller-connections
for ii = 1:N
    for jj = 1:nBins(ii)
        if(rollers(ii).connections(jj) ~= 0)
            index1 = rollers(ii).connections(jj);
            index2 = find(ii == rollers(index1).connections);
            ink1 = rollers(ii).ink(jj);
            ink2 = rollers(index1).ink(index2);
            rollers(ii).ink(jj) = (ink1+ink2)/2;
            rollers(index1).ink(index2) = (ink1+ink2)/2;
        end
     end
end

% # Calculate average amount of ink on each roller
    for ii = 1:N
        averageAmountOfInk(tt,ii) = sum(rollers(ii).ink);
    end
end

    image(5:20) = (rollers(7).ink(5:20))./2;
    inkPerSeg1 = [rollers(1).ink]';
    inkPerSeg2 = [rollers(2).ink]';
    inkPerSeg3 = [rollers(3).ink]';
    inkPerSeg4 = [rollers(4).ink]';
    inkPerSeg5 = [rollers(5).ink]';
    inkPerSeg6 = [rollers(6).ink]';
    inkPerSeg7 = [rollers(7).ink]';
4

1 回答 1

1

这是一个扩展评论，而不是一个正确的答案，但评论框有点太小了......

你的代码让我不知所措，我只见树木不见森林。我建议你消除所有我们不需要看到的东西来帮助你解决你的直接问题（例如所有那些画线图）——我认为它会帮助你自己调试你的代码，把所有的东西放进去函数或脚本。

你的代码片段

For k = 1:nTimeSteps
    While nTimesSteps = mod(nTimeSteps, nBins_max) == 0 
    ‘Give me the ink on each segment at each time step in a matrix’
End

可能是（我不太明白你对while语句的使用，这个词While不是 Matlab 关键字，并且正如你所写的那样，语句返回的值不会随着迭代而改变）相当于

For k = 1:nBins_max:nTimeSteps
    ‘Give me the ink on each segment at each time step in a matrix’
End

您似乎错过了 Matlab 冒号运算符的一个基本功能......

1:8 = [1 2 3 4 5 6 7 8]

但

1:2:8 = [1 3 5 7]

也就是说，三元组中的第二个数字是连续元素之间的步幅。

您的矩阵在连接滚轮conn的位置为 1，在其他位置为 0。您可以像这样找到所有 s(row,col)的行和列索引：1

[ri,ci] = find(conn==1)

然后，您可以在没有循环嵌套和 if 语句开始的情况下获取 s的(row,col)位置1

for ii = 1:N
    for jj = 1:N
        if(ii~=jj)
            if(conn(ii,jj) == 1)

我可以继续，但不会，这足以发表一条评论。

于 2013-03-21T22:50:18.990 回答