非線性電路中的混沌現(xiàn)象實驗報告-資料下載頁
2025-04-03 01:20本頁面
【正文】 [xx x]。 end 文件 2： : function [x0]=chua_initial(x,aaa) h=。a=h/2。aa=h/6。 x=[ ]39。 k0=chua_map(x,1,aaa)。 x1=x+k0*a。 k1=chua_map(xl,1,aaa)。 x1=x+k1*a。 k2=chua_map(x1,1,aaa)。 x1=x+k2*h。 k3=chua_map(x1,1,aaa)。 x=x+aa*(k0+2*(kl+k2)+k3)。 for k=2:400 kO=chua_map(x,k,aaa)。 x1=x+k0*a。 k1=chua_map(x1,k,aaa)。 x1=x+k1*a。 基礎(chǔ)物理實驗報告 第 28 頁 k2=chua_map(x1,k,aaa)。 x1=x+k2*h。 k3=chua_map(xl,k,aaa)。 x=x+aa*(kO+2*(k1+k2)+k3)。 end x0=x。 文件 3： : function[x]=chua_map(xx,time_variable,aaa) m0=1/。 m1=2/。 if xx(1)=1 hx=m1*xx(1)+m0m1。 elseif abs(xx(1))=1 hx=m0*xx(1)。 else hx=m1*xx(1)m0+m1。 end A=[0 0 O aaa 0]。 x=A*xx。 x=x+[9*hx 0 O]39。 文件 4： x0=*randn(3,1)。 [x0]=chua_initial(x0,100/7)。 [xx]=chua(x0,1,100/7,20210)。 plot(UVI(1,1:end),UVI(2,1:end))。 xlabel(39。Uc1 (V)39。)。ylabel(39。Uc2 (V)39。)。 figure。 plot3(UVI(3,1:end),UVI(2,1:end),UVI(1,1:end)) xlabel(39。I (V)39。)。ylabel(39。Uc1 (V)39。)。zlabel(39。Uc2 (V)39。)。 (2) 基礎(chǔ)物理實驗報告 第 29 頁 對于本實驗，其微分方程組的求解還可以采用離散化的處理。具體代碼如下：（ Matlab 代碼） function discrete_chai dt=。 c1=1/9。 c2=1。 L=1/7。 G=。 N=10000。 a0=。a1=。 MT=[1dt*G/c1,dt*G/c1,0。dt*G/c2,(1dt*G/c2),dt/c2。0,dt/L,1]。 UVI=zeros(3,N)。 UVI(:,1)=[。]。 for k=1:N1。 Bd=[dt/c1*a0*UVI(1,k)*(a1^2*UVI(1,k)^2/31)。0。0]。 UVI(:,k+1)=MT*UVI(:,k)+Bd。 end plot(UVI(1,1:end),UVI(2,1:end))。 xlabel(39。Uc1 (V)39。)。ylabel(39。Uc2 (V)39。)。 figure。 plot3(UVI(3,1:end),UVI(2,1:end),UVI(1,1:end)) xlabel(39。I (V)39。)。ylabel(39。Uc1 (V)39。)。zlabel(39。Uc2 (V)39。)。 經(jīng)驗證：該代碼的執(zhí)行效率比 四階龍格庫塔數(shù)值積分法要高，但初始精度稍差。 基礎(chǔ)物理實驗報告 第 30 頁 （ 2）數(shù)值仿真結(jié)果： 改變 G 的值，當 G= 時，數(shù)值仿真出現(xiàn)雙吸引子： Uc1Uc2 圖 使用 matlab 的 Plot3 可以做出 IUc1Uc2 的三維圖： IUc1Uc2 圖 基礎(chǔ)物理實驗報告 第 31 頁 同時可以使用 Plot 做出 I、 Uc1 和 Uc2 對時間的曲線： 基礎(chǔ)物理實驗報告 第 32 頁 改變 G 值，使 G=，數(shù)值仿真出現(xiàn)單吸引子： Uc1Uc2 圖 使用 matlab 的 Plot3 可以做出 IUc1Uc2 的三維圖： 基礎(chǔ)物理實驗報告 第 33 頁 同時可以使用 Plot 做出 I、 Uc1 和 Uc2 對時間的曲線： 基礎(chǔ)物理實驗報告 第 34 頁 基礎(chǔ)物理實驗報告 第 35 頁 基礎(chǔ)物理實驗報告 第 36 頁
高考資料相關(guān)推薦
鄂ICP備17016276號-1