继续获取有关我的ODE函数的错误消息

6个视图（过去30天）

显示旧的评论

克里斯汀·辛西娅·汉塞尔 2019年11月27日

回答：克里斯汀·辛西娅·汉塞尔2019年12月2日

接受的答复： 1月

我目前正在一个间歇式反应器中对聚乳酸进行编码，不断收到错误信息。我的教授帮助我的小组达到了这一点，但每次我们运行它时，我们都会收到，

使用horzcat时出错

连接的数组的维度不一致。

我们已经想尽一切办法来解决这个问题，但我们不知道为什么会发生这种情况，也不知道如何解决。

sad_girl_项目中的错误（第46行）

u0=[P0零（N，1）0]；

                    作用悲伤女孩项目
                   
                    T=180+273.15；%温度[K]
                   
                    Keq=10；
                   
                    %%计算
                   
                    lk1=-（6810.3/T）+6.1928
                   
                    k1 = exp (lk1)聚合速率常数(L/mol*s)
                   
                    km1=k1/Keq%解聚速率常数（L/mol*s）
                   
                    P0=10；
                   
                    N=100；%极化程度
                   
                    %u=[P1，P2，….PN，W]
                   
                    作用dudt=deq（t，u）
                   
                    P=u（1:N）；
                   
                    W=u（N-1）；
                   
                    sum3=总和（P）；
                   
                    sum4=0；
                   
                    对于ii=1:N
                   
                    sum1=0；
                   
                    对于jj=1:ii-1
                   
                    sum1=sum1+P（ii-jj）*P（jj）；
                   
                    终止%对于jj
                   
                    %disp（[ii，P（1:2），sum1]）
                   
                    sum2 =总和(P (2 + 1: N));
                   
                    dudt（ii）=k1*sum1+2*km1*W*sum2-...
                   
                    2 * k1 * P (ii) * sum3——km1 * W * (ii-1) * P (ii);
                   
                    sum4=sum4+（ii-1）*P（ii）；
                   
                    %disp（[ii，sum1，sum2，sum3]）
                   
                    终止%第二次
                   
                    dudt（N+1）=k1*sum3^2-km1*W*sum4；
                   
                    dwdt（N+1）=0；
                   
                    dudt=[dudt；dwdt]；
                   
                    终止%德克
                   
                    m=[024*3600]；
                   
                    u0=[P0零（N，1）0]；
                   
                    [t u] = ode15 (@deq, m,情况);
                   
                    %%密谋
                   
                    图（1）
                   
                    图（t，u）
                   
                    头衔(“时间与分子量”)
                   
                    xlabel(“时间[]”)
                   
                    伊拉贝尔(“浓度[]”)
                   
                    终止