Well, I confess I don’t follow the intentions of your code, but I think you may
have some misunderstandings of what we mean by a “scenario” in MOST. Think of
it as a possible “state” of the system at a particular time (not a set of
trajectories through time). So to model wind, you will want to attempt to model
the joint distribution of all of your wind locations as a set of probability
weighted “scenarios” or states for that moment in time.
Ray
> On Sep 14, 2018, at 8:52 AM, Aamir Nawaz <[email protected]> wrote:
>
> Dear all,
> I want to find expected dispatch for IEEE 188 Bus system with 30% Wind farms
> where wind scenrios is nj=400 and for time period nt=12 hours. Please help
> where this code have problem.
>
> define_constants
> mpopt = mpoption('verbose', 0);
> mpc = loadcase(case118);
> nWF=fix(sum(mpc.gen(:,9))*0.3/100);%30% wind farms
> PVbus=find(mpc.bus(:,BUS_TYPE)~=1);
> xgd = loadxgendata(ex_xgd118, mpc);
> for i=1:nWF
> [iwind(i), mpc, xgd] = addwind(ex_wind_dynmc118, mpc,xgd);
> end
> nt =12;nj=400;
> weibullpdf=wblpdf(0.01:0.01:nj*0.01,1:nj,1);
> transmat = {repmat(weibullpdf,nt,1)};
> profiles = getprofiles(uniformwindprofile(nj, nt), iwind);
> profiles = getprofiles('ex_load_profile118', profiles);
> mdi = loadmd(mpc, transmat, xgd, [], [], profiles);
> mpopt = mpoption('verbose', 0, 'most.dc_model', 1);
> mdo = most(mdi, mpopt);
> EPg = mdo.results.ExpectedDispatch
>
> But its giving expected dispatch for 1 hour only, instead of 12 hours. Can
> you point out my mistake?
>
>
>
> Regards,
> Aamir Nawaz
