Problem 04 Rev 1 - Arx

5

‫ﺷﮕﺮﺩﻫﺎﻱ ﻋﺪﺩﻱ ﺩﺭ ﭘﺎﺳﺨﻴﺎﺑﻲ ﭘﺮﺳﺸﻬﺎﻱ ﻣﻬﻨﺪﺳﻲ‬

‫ ﻛﻪ ﺩﺭ ﺁﻥ ﺩﺭﻭﻧﺪﺍﺩ ﺩﺭ ﺩﻭ ﺍﻧﺪﺍﺯﻩﻱ ﺑﺨﺘﻲ ﺗﻐﻴﻴﺮ‬،‫ ﻧﻘﻄﻪﻱ ﺩﺭﻭﻧﺪﺍﺩ ـ ﺑﺮﻭﻧﺪﺍﺩ ﺍﺯ ﻳﻚ ﻓﺮﺍﻳﻨﺪ ﺭﺍ‬70 ‫ ـ ﮔﻤﺎﻥ ﻛﻨﻴﺪ‬7.New.1 ‫ﻧﻤﻮﻧﻪ‬ y ،‫ ﺩﺭﻭﻧﺪﺍﺩ ﺩﺳﺘﮕﺎﻩ‬u) ‫ ﺍﻳﻦ ﺩﺍﺩﻩﻫﺎ‬.‫ ﺍﺯ ﺁﺯﻣﺎﻳﺶ ﺑﺮ ﺭﻭﻱ ﻳﻚ ﺩﺳﺘﮕﺎﻩ ﻣﻮﺧﺸﻚﻛﻦ ﺩﺳﺘﻲ ﺑﻪ ﺩﺳﺖ ﺁﻭﺭﺩﻩ ﺑﺎﺷﻴﻢ‬،‫ﻣﻲﻛﻨﺪ‬

‫ ﺑﺮﺍﻱ ﺷﻨﺎﺳﺎﻳﻲ ﻭ ﭘﻴﺶﺑﻴﻨﻲ ﺭﻓﺘﺎﺭ‬.‫ ﺫﺧﻴﺮﻩ ﺷﺪﻩﺍﻧﺪ‬d20301 ‫ﻓﺎﻳﻞ‬mat ‫ ﻧﻤﻮ ﺯﻣﺎﻥ ﻧﻤﻮﻧﻪﺑﺮﺩﺍﺭﻱ ﺍﺳﺖ( ﺩﺭ‬T ‫ﺑﺮﻭﻧﺪﺍﺩ ﺩﺳﺘﮕﺎﻩ ﻭ‬ :‫ ﺭﺍ ﻧﻮﺷﺘﻪﺍﻳﻢ‬prdct ‫ ﺗﺎﺑﻊ‬،‫ﺩﺳﺘﮕﺎﻩ‬ function [tsim,usim,ysim]=prdct(t,u,y,T,ni,np) %PRDCT Special Sample of ARX Model for System Identification. % [TSIM,USIM,YSIM]=PRDCT(T,U,Y,SAMPLT,NI,NP) makes predicted output % for sampling data of a test. % % Output: % TSIM: Time Samples % USIM: Random Input % YSIM: Estimated Output of Simulation % Input: % T : Test Time Samples % U : Test Input Samples % Y : Test Output Samples % SAMPLT: Sampling Time Increment % NI : Number of Identification Samples % NP : Number of Predicted Samples % % Numerical Recipes for Solving Engineering Problems using MATLAB % By: Mehdi Mosafer, 2nd ed., 2008. t=((1:length(u))-1)*T; A=[-y(4:ni-1),-y(3:ni-2),u(2:ni-3),u(1:ni-4)]; b=y(5:ni); x=A\b usim=[u;(max(u)-min(u))*round(rand(np,1))+min(u)]; tsim=((1:length(usim))-1)*T; ysim=zeros(size(usim)); ysim(1:4)=y(1:4); n=length(usim); for k=5:n ysim(k)=-x(1)*ysim(k-1)-x(2)*ysim(k-2)+x(3)*usim(k-3)+x(4)*usim(k-4); end

‫ ﺑﺨﺸﻲ ﺍﺯ ﺩﺍﺩﻩﻫﺎ ﺭﺍ ﺑﺮﺍﻱ ﺷﻨﺎﺳﺎﻳﻲ ﻭ ﺑﺨﺸﻲ ﺭﺍ‬.‫ﺑﺎ ﺍﺳﺘﻔﺎﺩﻩ ﺍﺯ ﺍﻳﻦ ﺗﺎﺑﻊ ﻧﻤﻮﺩﺍﺭ ﺩﺭﻭﻧﺪﺍﺩ ﻭ ﺑﺮﻭﻧﺪﺍﺩ ﺭﺍ ﻧﺴﺒﺖ ﺑﻪ ﺯﻣﺎﻥ ﺑﻜﺸﻴﺪ‬ ‫ ﺍﻓﺰﻭﻥ ﺑﺮ ﺁﻥ ﺑﻪ ﻳﺎﺭﻱ ﻋﺪﺩﻫﺎﻱ ﺑﺨﺘﻲ ﺍﻧﺪﻛﻲ ﭘﺲ ﺍﺯ ﭘﺎﻳﺎﻥ ﺯﻣﺎﻥ ﺁﺯﻣﺎﻳﺶ‬.‫ﺑﺮﺍﻱ ﺁﺯﻣﻮﺩﻥ ﺩﺭﺳﺘﻲ ﻣﺪﻝ ﺷﻨﺎﺧﺘﻪﺷﺪﻩ ﺑﻪ ﻛﺎﺭ ﺑﺮﻳﺪ‬ .‫ﺭﺍ ﭘﻴﺶﺑﻴﻨﻲ ﻛﻨﻴﺪ‬