《最新Modeling and simulation of systems using MATLAB and Simulink Chapter 2教學(xué)講義ppt課件.ppt》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、ModelingandsimulationofsystemsusingMATLABandSimulinkChapter2ProblemSituationProblemDefinitionDirectExperimentationModelBuildingModelSimulationAnalysisDecisionMakingModelBuildingStep–1ObservetherealsystemanditsinteractionsbetweencomponentsCollectdataonitsbehavior.Consult

2、totheexpertsfordetailedknowledge.Asmodeldevelopmentproceeds,newquestionsmayraiseandmodeldeveloperwillreturntothisstep.ContinuousSystemSimulation(WaterReservoirSystem)ReservoirRiverInflowDirectRainfallSeepagelossEvaporationlossSupplytomeetdemandProblemDefinitionDetermin

3、ethevolumeofwaterstoredinthereservoir,sothatitwillnotspilloverornoshortageinsummer.Seepagelossisfunctionofvolume.Seepageloss=f(Volume)Evaporationlossisfunctionofareaofexposedsurfaceandcoefficientofevaporation.-Evaporationloss=f(SurfaceArea,Cevap)Thestatevariableisthevol

4、umewhichchangeswithtime.ModelSimulationtime=100years(1200months)Stepsize=1monthCalculatetotalinputVin=Rain+RiverFlowGrossVolume=Vin+EarlierVolumeTotalloss=Seepage+EvaporationNetVolume=GrossVolume–TotalLossIfDemand>=NetVolumeThenReservoirisdryandshortageofwater(=Dem-Vnet

5、)OtherwiseDiff=Demand–VnetIfDiff>CapThenSpillover.AChemicalReactorInacertainchemicalreactionsubstancesAandBproducethirdchemicalsubstanceC.A(1gm)+B(1gm)C(2gm)RateofformationofCproportionaltoamountsofAandB(Forwardreaction).RateofCDecompositionofCisproportionaltoamountofCp

6、resentinthemixer(Backwardreaction).ModelDifferentialEquationsda/dt=K2c–K1abdb/dt=K2c–K1abdc/dt=2K1ab–2K2cInitialconditionsK1,K2,a,baregiven.c=0,t=0,choosesuitablevalueofdt.Usesomeintegrationtechniquetogetthetimeprofileofa,bandc.SelectionofIntegrationMethodEuler’smethodS

7、implestmethodofintegrationVnew=Vold+dt*dv/dt.Thismethodisrarelyusedbecauseofrapidaccumulationoferrorandalsonotverystable.RungeKuttamethodVnew=Vold+K2whereK2=dt*f(tn+0.5dt,V+0.5k1)K1=dt*f(tn,V)FourthOrderRungaKuttaMethodVnew=Vold+K1/6+K2/3+K3/3+K4/6whereK1=dt*f(tn,V)K2=d

8、t*f(tn+0.5dt,V+0.5K1)K3=dt*f(tn+0.5dt,V+0.5K2)K4=dt*f(tn+dt,V+K3)Simpson’sRuleOne-ThirdRuleV=dt*(K1+4K2+K3)/3w