原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
基础物理实验报告
第 PAGE 30页
非线性电路中的混沌现象
学号
姓名：蔡正阳
日期：2009年3月2
五：数据处理：
1．计算电感L
本实验采用相位测量。根据RLC谐振规律，当输入激励的频率时，RLC串联电路将达到谐振，L和C的电压反相，在示波器上显示的是一条过二四象限的45度斜线。
测量得：f=32.8kHz；实验仪器标示：C=1.095nF
由此可得：
估算不确定度：
估计u(C)=0.005nF，u(f)=0.1kHz
则：
即
最终结果：
2．用一元线性回归方法对有源非线性负阻元件的测量数据进行处理：
（1）原始数据：
R
V
R
V
R
V
71200
-12
2044.9
-8
1753.4
-4
21000
-11.8
2036.2
-7.8
1727.5
-3.8
12150
-11.6
2027.2
-7.6
1699.6
-3.6
8430
-11.4
2017.8
-7.4
1669.4
-3.4
6390
-11.2
2007.9
-7.2
1636.7
-3.2
5100
-11
1997.5
-7
1601.2
-3
4215
-10.8
1986.7
-6.8
1562.4
-2.8
3564
-10.6
1975.3
-6.6
1519.7
-2.6
3070
-10.4
1963.4
-6.4
1472.3
-2.4
2680
-10.2
1950.9
-6.2
1420
-2.2
2369
-10
1937.6
-6
1360.9
-2
2115
-9.8
1923.7
-5.8
1295.1
-1.8
2103.1
-9.6
1909
-5.6
1281.8
-1.6
2096.8
-9.4
1893.4
-5.4
1276.7
-1.4
2090.2
-9.2
1876.9
-5.2
1270.1
-1.2
2083.4
-9
1859.5
-5
1261.1
-1
2076.3
-8.8
1840.9
-4.8
1247.8
-0.8
2068.9
-8.6
1821.2
-4.6
1226
-0.6
2061.2
-8.4
1800.1
-4.4
1148.9
-0.4
2053.3
-8.2
1777.6
-4.2
1075
-0.2
VRR1
V
R
R1
根据可以得出流过电阻箱的电流，由回路KCL方程和KVL方程可知：
由此可得对应的值。
对非线性负阻R1，将实验测得的每个（I，U）实验点均标注在坐标平面上，可得：
图中可以发现，（0.0046336，-9.8）和（0.0013899，-1.8）两个实验点是折线的拐点。故我们在、、这三个区间分别使用线性回归的方法来求相应的I-U曲线。
使用Excel的Linest函数可以求出这三段的线性回归方程：
经计算可得，三段线性回归的相关系数均非常接近1（r=0.99997），证明在区间内I-V线性符合得较好。
应用相关作图软件可以得出非线性负阻在U0区间的I-U曲线。
将曲线关于原点对称可得到非线性负阻在U0区间的I-U曲线：
3．观察混沌现象：
（1）一倍周期：
一倍周期 Vc1-t
（2）两倍周期：
两倍周期 Vc1-t
（3）四倍周期：
四倍周期 Vc1-t
（4）单吸引子：
单吸引子 阵发混沌
三倍周期 Vc1-t
(5)双吸引子：
双吸引子 Vc1-t
4．使用计算机数值模拟混沌现象：
（1）源程序（Matlab代码）：
算法核心：四阶龙格库塔数值积分法
文件1：chua.m
function [xx]=chua(x,time_variable,aaa,symbol_no)
h=0.01;
a=h/2;
aa=h/6;
xx=[];
for j=1:symbol_no;
k0=chua_map(x,time_variable,aaa);
x1=x+kO*a;
k1=chua_map(xl,time_variable,aaa);
xl=x+k1*a;
k2=chua_map(x1,time_variable,aaa);
x1=x+k2*h;
k3=chua_map(x1,time-variable,aaa);
x=x+aa*(kO+2*(k1+k2)+k3);
xx=[xx x];
end
文件2：chua_initial.m:
function [x0]=chua_initial(x,aaa)
h=0.01;a=h/2;aa=h/6;
x=[-0.03 0.6 -0.01];
k0
文档评论（0）