搜索内容 搜索用户

# 14 麦片

050

10主题 0最佳答案

 >> syms  d22 b22   %a是前浓缩比 for ce0=1:0.5:4   a=0.5;   c0=10;   v0=16;   n1=c0*v0-a*v0*ce0;   h=(c0*v0-n1)/(c0*v0);     for x=0.01:h/5:h      ce0x=ce0-x*c0/a;      b1=-log(ce0/c0);      d1=-log(ce0x/ce0);      b1x=-log(ce0x/c0);      f1=b1+d1;      c2=vpa(x*c0/d1);         for i=2:1:5;            c(i)=c2;            cei1=c(i)-(c(i)-c0)*exp(-d1);             b2=vpasolve(cei1*exp(-(b22-d1))-ce0==0,b22);             d2=vpasolve(cei1*exp(-(b2+d22-d1))-ce0x==0==0,d22);             cei2=cei1*exp(-(b2-d1));             cei2x=cei1*exp(-(b2+d2-d1));             ni2=c0*v0+x*c0*v0-a*v0*cei2;             ni2x=c0*v0+x*c0*v0-a*v0*cei2x;             fprintf('%d %4.2f %d %f %f %f %f %f %f %f %f %f\n',n1,x,i,c(i),d1,cei1,b2,d2,cei2,cei2x,ni2,ni2x);             c2=vpa((ni2x-ni2)/(a*v0*d1));          end    end end 152 0.01 2 0.448142 0.223144 8.089628 2.313726 0.223144 1.000000 0.800000 153.600000 155.200000 152 0.01 3 0.896284 0.223144 8.179257 2.324745 0.223144 1.000000 0.800000 153.600000 155.200000 152 0.01 4 0.896284 0.223144 8.179257 2.324745 0.223144 1.000000 0.800000 153.600000 155.200000 152 0.01 5 0.896284 0.223144 8.179257 2.324745 0.223144 1.000000 0.800000 153.600000 155.200000 152 0.02 2 0.391523 0.510826 6.156609 2.328352 0.510826 1.000000 0.600000 155.200000 158.400000 152 0.02 3 0.783046 0.510826 6.313218 2.353471 0.510826 1.000000 0.600000 155.200000 158.400000 152 0.02 4 0.783046 0.510826 6.313218 2.353471 0.510826 1.000000 0.600000 155.200000 158.400000 152 0.02 5 0.783046 0.510826 6.313218 2.353471 0.510826 1.000000 0.600000 155.200000 158.400000 152 0.03 2 0.327407 0.916291 4.196444 2.350528 0.916291 1.000000 0.400000 156.800000 161.600000 152 0.03 3 0.654814 0.916291 4.392888 2.396278 0.916291 1.000000 0.400000 156.800000 161.600000 152 0.03 4 0.654814 0.916291 4.392888 2.396278 0.916291 1.000000 0.400000 156.800000 161.600000 152 0.03 5 0.654814 0.916291 4.392888 2.396278 0.916291 1.000000 0.400000 156.800000 161.600000 152 0.04 2 0.248534 1.609438 2.198827 2.397362 1.609438 1.000000 0.200000 158.400000 164.800000 152 0.04 3 0.497068 1.609438 2.397654 2.483929 1.609438 1.000000 0.200000 158.400000 164.800000 152 0.04 4 0.497068 1.609438 2.397654 2.483929 1.609438 1.000000 0.200000 158.400000 164.800000 152 0.04 5 0.497068 1.609438 2.397654 2.483929 1.609438 1.000000 0.200000 158.400000 164.800000 152 0.05 2 0.000000 Inf 0.000000       在赋值 A(:) = B 中，A 和 B 中的元素数目必须相同。

# 论坛优秀回答者

4主题 189最佳答案
• 关注者： 92

 已经说得很清楚了啊 A(:)如果是n个的，那B也必须是n个

# 14 麦片

050

10主题 0最佳答案
楼主| 发表于 4 天前 | 显示全部楼层
 coolchen302 发表于 2020-9-15 16:39 已经说得很清楚了啊 A(如果是n个的，那B也必须是n个 我把c2改成c1就好了;P，虽然还是有点懵，但至少弄出来了

# 14 麦片

050

10主题 0最佳答案
楼主| 发表于 前天 09:45 | 显示全部楼层
 还是这个问题 syms  d22 b22  c1 %a是前浓缩比 for ce0=1:0.5:4   a=0.5;   c0=10;   v0=16;   n1=c0*v0-a*v0*ce0;   h=(c0*v0-n1)/(c0*v0);     for x=0.01:0.02:h      ce0x=ce0-x*c0/a;      b1=-log(ce0/c0);      d1=-log(ce0x/ce0);      b1x=-log(ce0x/c0);      f1=b1+d1; c1=x*c0/d1/a;       for i=2:1:5        c(i)=c1;            cei1=c(i)-(c(i)-c0)*exp(-d1);             b2=vpasolve(cei1*exp(-(b22-d1))-ce0==0,b22);             d2=vpasolve(cei1*exp(-(b2+d22-d1))-ce0x==0==0,d22);             cei2=cei1*exp(-(b2-d1));             cei2x=cei1*exp(-(b2+d2-d1));             ni2=c0*v0+x*c0*v0-a*v0*cei2;             ni2x=c0*v0+x*c0*v0-a*v0*cei2x; fprintf('%f %f    \n',x,c1);           c1=(ni2x-ni2)/(a*v0*i); end    end end 0.010000 0.896284    0.010000 0.100000    0.010000 0.066667    0.010000 0.050000    0.030000 0.654814    0.030000 0.300000    0.030000 0.200000    0.030000 0.150000    0.050000 0.000000    在赋值 A(:) = B 中，A 和 B 中的元素数目必须相同。 >>

# 14 麦片

050

10主题 0最佳答案
楼主| 发表于 前天 09:46 | 显示全部楼层
 agoodgirl 发表于 2020-9-17 09:45 还是这个问题 syms  d22 b22  c1 %a是前浓缩比 我先给c(i)一个初始值，在一个循环后，重新赋予一个值，为什么说两个数组的元素不一样呢？

 您需要登录后才可以回帖 登录 | 注册 本版积分规则 回帖后跳转到最后一页