没有合适的资源?快使用搜索试试~ 我知道了~
非线性电路中的混沌现象实验报告.docx
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 1 下载量 114 浏览量
2022-05-29
15:05:22
上传
评论
收藏 1.16MB DOCX 举报
温馨提示
试读
21页
。。。
资源推荐
资源详情
资源评论
如有你有帮助,请购买下载,谢谢!
非线性电路中的混沌现象
学号:
姓名:蔡正阳
日期:2009 年 3 月 24 日
五:数据处理:
1.计算电感 L
本实验采用相位测量。根据 RLC 谐振规律,当输入激励的频率
1
f
时,RLC 串联电路将达到谐振,L 和 C 的电压反相,在示
2 LC
波器上显示的是一条过二四象限的 45 度斜线。
测量得:f=32.8kHz;实验仪器标示:C=1.095nF
由此可得:
估算不确定度:
估计 u(C)=0.005nF,u(f)=0.1kHz
则:
u(L) 0.16mH
即
L u(L) (21.5 0.2)mH
最终结果:
2.用一元线性回归方法对有源非线性负阻元件的测量数据进行处理:
(1)原始数据:
V
71200
21000
12150
8430
6390
5100
4215
3564
3070
2680
2369
2044.9
2036.2
2027.2
2017.8
2007.9
1997.5
1986.7
1975.3
1963.4
1950.9
1937.6
-8
1753.4
1727.5
1699.6
1669.4
1636.7
1601.2
1562.4
1519.7
1472.3
1420
-4
-3.8
-3.6
-3.4
-3.2
-3
-10.8
-10.6
-10.4
-10.2
-10
-6.8
-6.6
-6.4
-6.2
-6
-2.8
-2.6
-2.4
-2.2
-2
1360.9
1 页
如有你有帮助,请购买下载,谢谢!
1295.1
1281.8
1276.7
1270.1
1261.1
1247.8
1226
-1.8
-1.6
-1.4
-1.2
-1
2103.1
2096.8
2090.2
2083.4
2076.3
2068.9
2061.2
2053.3
1893.4
1876.9
1859.5
1840.9
1821.2
1800.1
1777.6
-8.8
-8.6
-8.4
-8.2
-4.8
-4.6
-4.4
-4.2
-0.8
-0.6
-0.4
-0.2
1148.9
1075
(2)数据处理:
R1
R
U
I
根据
R
可以得出流过电阻箱的
R
R
电流,由回路 KCL 方程和 KVL 方程可知:
由此可得对应的 I 值。
R1
对非线性负阻R1,将实验测得的每个(I,
U)实验点均标注在坐标平面上,可得:
图中可以发现,(0.0046336,-9.8)和
V
( 0.0013899 , -1.8 ) 两 个 实 验 点 是 折 线 的 拐 点 。 故 我 们 在
12 U 9.8V
、 9.8 U 1.8V 1.8 U 0V
、
这三个区间分
别使用线性回归的方法来求相应的 I-U 曲线。
使用 Excel 的 Linest 函数可以求出这三段的线性回归方程:
经计算可得,三段线性回归的相关系数均非常接近1(r=0.99997), 证
明在区间内 I-V 线性符合得较好。
应用相关作图软件可以得出非线性负阻在 U<0 区间的 I-U 曲线。
将曲线关于原点对称可得到非线性负阻在 U>0 区间的 I-U 曲线:
2 页
如有你有帮助,请购买下载,谢谢!
3.观察混沌现象:
(1)一倍周期:
一倍周期
Vc -t
1
(2)两倍周期:
两倍周期
Vc -t
1
(3)四倍周期:
四倍周期
(4)单吸引子:
单吸引子
Vc -t
1
阵发混沌
三倍周期
Vc -t
1
(5)双吸引子:
双吸引子
Vc -t
1
3 页
如有你有帮助,请购买下载,谢谢!
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=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;
4 页
如有你有帮助,请购买下载,谢谢!
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:chua_map.m:
function[x]=chua_map(xx,time_variable,aaa)
m0=-1/7.0;
m1=2/7.0;
if xx(1)>=1
hx=m1*xx(1)+m0-m1;
elseif abs(xx(1))<=1
hx=m0*xx(1);
else
hx=m1*xx(1)-m0+m1;
end
A=[0 9.0 0
1.0 -1.0 1.0
O aaa 0];
x=A*xx;
x=x+[-9*hx 0 O]';
文件 4:chua_demo.m
x0=0.05*randn(3,1);
[x0]=chua_initial(x0,-100/7);
[xx]=chua(x0,1,-100/7,20000);
plot(UVI(1,1:end),UVI(2,1:end));
xlabel('Uc1 (V)');ylabel('Uc2 (V)');
figure;
plot3(UVI(3,1:end),UVI(2,1:end),UVI(1,1:end))
xlabel('I (V)');ylabel('Uc1 (V)');zlabel('Uc2 (V)'); (2)
对于本实验,其微分方程组的求解还可以采用离散化的处理。具体代
码如下:(Matlab 代码)
function discrete_chai
dt=0.04;
c1=1/9;
c2=1;
5 页
剩余20页未读,继续阅读
资源评论
- lxnsuwod2024-01-02感谢大佬分享的资源,对我启发很大,给了我新的灵感。
苦茶子12138
- 粉丝: 1w+
- 资源: 6万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- Picasso_v3.1 2.ipa
- chromedriver-mac-arm64.zip
- 蓝zapro.apk
- chromedriver-linux64.zip
- UCAS研一深度学习实验-MNIST手写数字识别python源码+详细注释(高分项目)
- 基于Python和PyTorch框架完成的一个手写数字识别实验源码(带MINIST手写数字数据集)+详细注释(高分项目)
- 基于Matlab在MNIST数据集上利用CNN完成手写体数字识别任务,并实现单层CNN反向传播算法+源代码+文档说明(高分项目)
- NVIDIA驱动、CUDA和Pytorch及其依赖
- 基于SVM多特征融合的微表情识别python源码+项目说明+详细注释(高分课程设计)
- html动态爱心代码一(附源码)
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功