標(biāo)題: matlab太陽(yáng)能光伏并網(wǎng)逆變仿真 [打印本頁(yè)]
作者: wee770    時(shí)間: 2019-10-28 11:46
標(biāo)題: matlab太陽(yáng)能光伏并網(wǎng)逆變仿真
matlab(2016a版本)-太陽(yáng)能光伏并網(wǎng)逆變仿真

單片機(jī)源程序如下:
  1. function [sys,x0,str,ts] = sfun_PV_array_MPPT(t,x,u,flag)
  2. switch flag,
  3.   case 0,
  4.     [sys,x0,str,ts] = mdlInitializeSizes;
  5.   case 2,                                                
  6.     sys = mdlUpdate(t,x,u);
  7.   case 3,                                                
  8.     sys = mdlOutputs(t,x,u);
  9.   case {1,4,9},                                                
  10.     sys =[];
  11.   otherwise
  12.     error(['unhandled flag = ',num2str(flag)]);
  13. end

  15. function [sys,x0,str,ts] = mdlInitializeSizes
  16. sizes = simsizes;
  17. sizes.NumContStates  = 0;
  18. sizes.NumDiscStates  = 3;
  19. sizes.NumOutputs     = 1;
  20. sizes.NumInputs      = 2;
  21. sizes.DirFeedthrough = 1;
  22. sizes.NumSampleTimes = 1;
  23. sys = simsizes(sizes);
  24. x0  =[0,0,0];
  25. str =[];
  26. ts  =[1e-6 0];

  28. function sys = mdlUpdate(t,x,u)
  29. k1=0.01;
  30. dV=u(1)-x(1);
  31. dI=u(2)-x(2);
  32. if dV==0
  33.     if dI==0
  34.     Vref = x(3);
  35.     elseif dI > 0
  36.     Vref = x(3)+k1;            
  37.     elseif  dI < 0
  38.     Vref = x(3)-k1;  
  39.     end
  40. else
  41.     if  (dI/dV) ==(-u(2)/u(1))
  42.         Vref = x(3);   
  43.     elseif (dI/dV) >(-u(2)/u(1))
  44.         Vref = x(3)+k1;   
  45.     elseif  (dI/dV) <(-u(2)/u(1))
  46.         Vref = x(3)-k1;
  47.     end
  48. end
  49. x(1)=u(1);
  50. x(2)=u(2);
  51. x(3)=Vref;
  52. sys=[x(1),x(2),x(3)];

  54. function sys = mdlOutputs(t,x,u)

  56. sys=x(3);
復(fù)制代碼

所有資料51hei提供下載:
太陽(yáng)能光伏并網(wǎng)逆變仿真.rar (14.22 KB, 下載次數(shù): 62)


作者: 時(shí)光之誠(chéng)    時(shí)間: 2020-4-9 23:09
謝謝樓主分享



歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1