20801008 发表于 2009-9-9 13:13

关于分岔的困惑

最近一直在研究分岔的问题，程序也编了很多，感觉还是问题多多，出现的结果和理想的还是有很大的差距，不知道能否有高人给予指点，谢谢了
下面是我的程序
vi = 10;
L = 0.001;
C = 12e-6;
T = 100e-6;
R = 20;
k = 1/(2*R*C);
w = sqrt(1/(L*C)-k^2);
for iref = 0.5:0.01:5.5
    vn = 5;
    in = 0.25;
   
      for n = 1:500
            tn = abs(L*(iref-in)/vn)
            tnn = T-tn;
            if tn<T
            vn = vi-exp(-k*tnn)*((k*vn*exp(-2*k*tn)-k*vi-(iref-vi/R)/C)*sin(w*tnn)/w+(vi-vn*exp(-2*k*tn))*cos(w*tnn));
            in = exp(-k*tnn)*((k*L*(iref-vi/R)+vi-vn*exp(-2*k*tn))/(w*L)*sin(w*tnn)+(iref-vi/R)*cos(w*tnn))+vi/R;
            else
                in = in +vi*T/L;
                vn = vn*exp(-T/R*C);
            end
   
            if (in>0&&in<6)
               plot(iref,in,'k.','markersize',2)
               hold on
            end
      
      end
end


下面是我得到的图形

20801008 发表于 2009-9-9 13:16

下面是理想的图像

讨论求知 发表于 2009-9-9 16:37

我也顶一下

无水1324 发表于 2009-9-11 16:38

vn = vi-exp(-k*tnn)*((k*vn*exp(-2*k*tn)-k*vi-(iref-vi/R)/C)*sin(w*tnn)/w+(vi-vn*exp(-2*k*tn))*cos(w*tnn));
            in = exp(-k*tnn)*((k*L*(iref-vi/R)+vi-vn*exp(-2*k*tn))/(w*L)*sin(w*tnn)+(iref-vi/R)*cos(w*tnn))+vi/R;

这里应该有问题吧，另外你的方程式一个离散系统？

20801008 发表于 2009-9-13 15:15

回复 地板 无水1324 的帖子

系统并非离散系统，只是将系统离散化了
因为看到资料上是将系统做了这样的处理
不知道这个方程的问题应该在什么地方

20801008 发表于 2009-9-13 15:21

还有另一个分岔图，不知道是不是分岔了，是关于vi的分岔

无水1324 发表于 2009-9-14 08:42

回复 6楼 20801008 的帖子

那你不要离散化，直接求解方程然后做分岔图，详细地看看论坛里面的帖子！
页: [1]
查看完整版本: 关于分岔的困惑