想请板上各位英雄一起玩玩看.m与.mdl联系问题的程序实现

武神裝攻

如果看不到图请滑鼠右键选显示图片(H) (像小弟按了几次就看到了.....)
不好意思 在各位百忙之中打扰一下  小弟纯是数值分析 不需用到专业
就是小弟近期用.m档写计算程式，有用到叠代方式求收敛值，最后不再改变收敛的值就是我的解答
当中有个数学式只局限单种机器，所以想到将值丢到simulink的.mdl档，但是目前有个问题就是
(1)我的.m档的值要如何丢给.mdl。
(2)这中间如何反覆执行?

小弟就想说先做一个简单的试验
就是每次叠代中都需要与.mdl做计算，那就假设一个矩阵里面的值就是每次.m档叠代后结果
clear all
clc
A = [10;12;14;16;18;20];
B = [11;13;15;17;19;21];
C = [10;13;14;17;18;21];
%(以上为例第一次ABC叠代结果是10 11 10，
                  第二次ABC叠代结果是12 13 13，
                  第三次ABC叠代结果是14 15 14，
                  第四次ABC叠代结果是16 17 17......依此类推
想用for回圈将每次ABC叠代结果分别传到如附件与图所示的test_mtomdl.mdl档内

将每次求得的输出out1  out2   out3
存到 I 矩阵
最后可以求到I=
I = [2  2.2  2;2.4  2.6  2.6;2.8  3  2.8;3.2  3.4  3.4;3.6  3.8  3.6;4  4.2  4.2];
%(以上为例第一次out123结果是2  2.2  2，
                   第二次out123结果是2.4  2.6  2.6，
                   第三次out123结果是2.8  3  2.8，......依此类推

想请问板上要如何利用.m与.mdl求 I ?    不好意思 在此先谢谢

武神裝攻

小弟想了一整天才想到.....
clear all
clc

A = [10;12;14;16;18;20];
B = [11;13;15;17;19;21];
C = [10;13;14;17;18;21];
ANSI = [2 2.2 2;2.4 2.6 2.6;2.8 3 2.8;3.2 3.4 3.4;3.6 3.8 3.6;4 4.2 4.2];

for k = 1:6
   
    set_param('test_mtomdl/DC Voltage Source','Amplitude','A(k,1)');
    set_param('test_mtomdl/DC Voltage Source1','Amplitude','B(k,1)');
    set_param('test_mtomdl/DC Voltage Source2','Amplitude','C(k,1)');
    sim('test_mtomdl')
    I(1,1) = out1(1,2);
    I(1,2) = out2(1,2);
    I(1,3) = out3(1,2);
    I2(k,:) = I;
end
小弟有在mdl上多加Clock原件
http://img.waffleimages.com/94eb9bbcec7477decc4b11c83d791812bb16b5be/01010102.JPG