基于MATLAB的各类混沌系统的计算机模拟(教学版)

一、概述Matlab是一种用于数学计算、数据分析和图形可视化的高级编程语言和工具，混沌理论是一种描述动态系统非线性行为的数学理论，而混沌分叉是混沌系统中特有的一种现象。

编写Matlab程序对混沌分叉进行研究即是将这两种领域结合起来，本文将介绍如何使用Matlab 编写混沌分叉程序。

二、混沌分叉理论简介混沌分叉理论是描述混沌系统行为的一个重要方面，其基本原理可简要描述如下：1. 系统的参数变化：在混沌系统中，当改变系统的某些参数时，系统的行为会发生变化，这种变化可能是渐变的，也可能是突然的。

2. 分岔现象：当系统的参数发生变化时，系统的稳定点可能会出现分叉现象，即稳定点从一个点分裂成多个点，这种分叉现象是混沌系统中一个显著的特征。

3. 分形结构：在混沌分叉中，分岔现象可能形成分形结构，这种结构在混沌系统中具有重要的理论和实际意义。

三、Matlab基础知识使用Matlab编写混沌分叉程序需要具备一定的Matlab基础知识，包括但不限于以下内容：1. Matlab基本语法：了解Matlab的基本语法规则，包括变量定义、数组操作、逻辑运算等。

2. Matlab图形绘制：掌握Matlab绘制图形的基本方法，包括绘制曲线、散点图等。

3. Matlab函数编写：了解如何在Matlab中编写自定义函数，并且能够熟练运用函数调用和参数传递。

四、混沌分叉程序的编写编写混沌分叉程序的基本步骤如下：1. 设定系统参数：首先需要定义混沌系统的参数，包括系统的微分方程、初始条件以及需要变化的参数。

2. 编写微分方程：根据所研究的具体混沌系统，编写系统的微分方程，通常为非线性微分方程。

3. 参数变化循环：对系统的某些参数进行变化，并且循环计算系统的轨迹，观察混沌分叉现象。

4. 图形绘制：绘制分叉图、分岔图等用于展示混沌分叉现象的图形。

五、程序示例下面为一个简单的混沌分叉程序示例，以具体的混沌系统为 Logistic映射为例，程序演示了 Logistic 映射参数 r 变化时轨迹的分叉现象。

二、系统整体结构框架本系统的基本结构，我们遵循你所提供的结构进行，整个框图如下所示：三、各个模块介绍和代码说明首先实现使用MATLAB 进行算法的仿真，我们得到的仿真结果如下所示： 运行MATLAB 程序：这个是混沌模型的基本仿真，说明公式以及初始值选取的正确性。

运行MATLAB程序：这个程序是混沌加密调制解调系统的MATLAB浮点仿真结果图，说明上面的结果是正确的。

运行MATLAB程序：从上面的仿真结果可知，以定点进行仿真，只要量化宽度满足一定要求之后，其完全不影响系统的精度。

根据上面的介绍，我们可以编写了如下的程序：自上而下，依次为：系统顶层文件——加密调制模块————加密子模块，lorenz混沌序列产生模块，组帧模块，并串模块。

——解密解调模块————解密子模块，Lorenz混沌序列产生模块，搜帧模块，串并模块。

下面对各个模块进行介绍。

接口含义分别为：系统时钟系统复位测试并行信号的产生使能信号测试并行信号加密模块的使能信号串行输出串行信号组帧输出加密输出解密输出解密输出信号的符号判决搜帧模块的相关峰输出搜帧模块的使能输出搜帧模块的数据输出最后串并转化的使能最后串并转化的数据输出顶层模块，主要功能就是产生测试数据，16位宽的并行数据：然后是分别调用加密调制模块和解密解调模块：具体的接口含义，我们在介绍这些子模块的时候进行说明。

加密模块的接口含义是时钟，复位，输入并行信号的使能信号，输入并行信号，串行输出，组帧后的输出，加密信号：其中：加密输出信号o_T_signal反馈给lorenz混沌模块进行迭代，迭代后的yn输出进行加密。

分别为系统时钟，复位，数据输入，混沌y输出，加密信号输出。

其基本原理，是通过信号和混沌信号的叠加来实现的。

其接口是系统时钟，复位，加密反馈信号，混沌xyz输出。

注意，这里我们采用的是基于浮点原理的定点计算方法，由于在FPGA中直接采用浮点计算算法比较复杂，而且比较耗资源，而传统的定点计算精度不够，因此，这里我们采用的是基于浮点的定点计算，即进行位宽扩展，相乘之后进行截位的方法。

常微分方程及常微分方程组地求解
( )
转载▼
分类：编程之
标签：
杂谈
最近参加了数学建模，对于老师说地算法地不同步长地精度不一样，编写了一个函数文件来实现这个精度地比较，把函数附上：文档收集自网络，仅用于个人学习
[] ()
这里使用可变输出输入函数地
{}为求解常微分方程地表达式
{}为求解常微分方程地定解条件
需要给出地变量有常微分方程地范围，({}{})
为对这个区间地分割({})文档收集自网络，仅用于个人学习
写于月日
取得算法需要地变量，并附上容易理解地含义变量
{};
{}自变量地范围
{}区间地分割次数
( )步长
{}; 常微分方程地表达式
{}; 常微分方程地定解条件表达式文档收集自网络，仅用于个人学习
首先求出所给常微分方程问题地精确解
();
();
; ;
(,'');
() ;
() (,{},{()});
:()
() () ;
() ((,{},{()}));
文档收集自网络，仅用于个人学习
利用方法求解近似数值微分解
();
();
;
() ;
() (,{},{})获得原方程地初解
:()
() () ;
() () .* (((),{},{()()}))特别记录中地字符串操作，提取子字符串即（：）...
文档收集自网络，仅用于个人学习
返回经过算法算出与地值
;
;
画图进行误差比较
(,'');
;
(,'');文档收集自网络，仅用于个人学习特此记录，以后写了新地算法再分享。

科技信息计算机与网络基于ｍａｔＩａｂ对混沌掩盖通信帕仿真硼究北京邮电大学信息与通信工程学院尉志青［摘要】本文利用ｍａｄ曲软件对混沌掩盖通信进行了仿真研究，具体研究了一维信号以及二维信号（数字图像）的混沌掩盖通信问题。

讨论了压缩系教对传输质量的影响。

并且对数字图像传输中的混沌加密进行了一些修正，使之加密效果更好。

［关键词】混沌加密（方法）混沌掩盖Ｌ０画如ｃ混沌序列、．混沌掩盖也称混沌遮掩或混沌隐藏。

混沌掩盖通信的基本原理是利用具有逼近丁｜高斯噪声统计特性的混沌信号对有用的信息进行混沌掩盖，接收端利用同步后的混沌信号进行解调，从而恢复出混沌掩盖之前的信号，接收有用信息。

实现通信的系统框图如图ｌ所示：图１混沌通信系统框图映射的混沌掩盖通信系统框图本文中所用的混沌掩盖方式为相加。

在接收端用同步后的混沌信号进行与之相廊的逆运算即可恢复出有用的信息。

下面针对基于Ｌ—ｇｉ８ｔ记映射的混沌掩盖通信系统进行仿真，研究混沌掩盖通信的一些特性。

ｋＩｇｉｓｔｉｃ混沌序列的迭代式为：ｘ（ｘ＋１）＝ＩＩ《ｎ）·（１一ｘ（ｎ））。

其中０≤《ｎ）≤ｌ，Ｏ≤ｕ≤４。

当ｕ取值在３．５８和４之间时是混沌的。

通信系统框图如图２图２基于Ｉ埘觚ｃ映射的混沌掩盖通信系统框图以上框图中ｓ（ｎ）为输入的信息序列，ｅ（ｎ）为已调序列，ｓ（ｎ）为解调输出序列。

令ｅ（ｎ）＝ｘ‘（ｎ＋１）：ｘ（ｎ＋１）＋ｋ×ｓ（ｎ）ｘ（ｎ＋１）：ｘ’（ｎ＋１）Ｍ０ｄ（１）ｋ为压缩系数，一般取ｋ＜１／５０。

我们知道，混沌掩盖的加法方式使消息信号浮于混沌载体之上，这就决定了只能采用很小的信号，而加性的小信号极易受到同样是加性的信道噪声的干扰，稍大的噪声会使接收端的误码率大大升高；如果用大的消息信号，该信号就会从掩盖的混沌码中浮出『ｆＩｉ在公开的信道上被窃取．所以ｋ不宜太大，也不宜太小。

ｘ（ｎ＋１）取ｘ’（ｎ＋１）对ｌ的模是为了使得调制输出的序列幅度不大于ｌ，可以保证ｋ画ｓｔｉｃ混沌序列不脱离混沌状态，可以抵御攻击者的破译。

基于Matlab的混沌特性分析1. 引言1.1 研究背景混沌理论起源于1960年代，是一种描述复杂系统行为的新理论，揭示了非线性系统中存在的一种无序、不可预测的动态行为。

混沌系统具有高度敏感性和非周期性，表现出随机性和确定性的结合，对于许多领域的研究具有重要的理论和实际意义。

在现代科学和工程领域，混沌系统的分析和控制已经成为一个热门的研究方向。

随着计算机技术的发展，基于Matlab的混沌特性分析方法成为研究混沌系统的有力工具。

Matlab提供了丰富的算法和库函数，可以方便地进行混沌系统建模、仿真和分析。

利用Matlab进行混沌特性分析，可以更深入地理解混沌系统的动力学行为，为系统的控制与优化提供理论支持。

1.2 研究目的研究目的的主要目标是通过基于Matlab的混沌特性分析，探讨混沌系统的特征和建模方法，并利用Matlab提供的分析工具对混沌系统进行详细分析。

通过深入研究混沌系统的特性和行为，可以更好地理解和预测混沌系统的运动规律和特点，为相关领域的研究和应用提供理论支持和参考依据。

本研究旨在探讨基于Matlab的混沌特性分析方法的有效性和可行性，为混沌系统的研究和应用提供一种新的分析途径和工具。

通过对混沌系统的特性进行深入分析和实验研究，可以揭示混沌系统背后的规律和内在机制，为相关领域的发展和应用提供新的思路和方法。

本研究的目的在于通过基于Matlab的混沌特性分析，深入探讨混沌系统的特性和行为，为相关领域的研究和应用提供新的视角和研究方法。

1.3 研究意义混沌系统在现代科学和工程中具有广泛的应用，例如在通信、控制、密码学等领域都有重要的作用。

对混沌系统进行特性分析，能够帮助我们更好地理解和掌握系统的行为规律，为系统的设计和优化提供重要的参考。

混沌系统的特性分析不仅可以帮助我们更好地理解系统的动态行为，还可以为混沌系统的控制和应用提供理论基础。

通过本文基于Matlab的混沌特性分析，我们可以更深入地探索混沌系统的特性和规律，为未来混沌系统的应用和发展提供重要参考。

混沌同步模型驱动系统和响应系统都是Lorenz System,只不过初值不同。

驱动系统： dx/dt=a*(y-x)dy/dt=r*x-y-xzdz/dt=x*y-b*z初值（0.1，0.1，0.1）输出信号令S(t)=x(t)响应系统：将S(t)代替x(t)作为激励信号dx/dt=a*(y-x)dy/dt=r*x-y-xzdz/dt=x*y-b*z初值（0.1，0.1，1）最后求响应系统的输出x(t),y(t),z(t)程序：function [Y1] = Lorenz_response(tspan);%%计算处于响应地位的Lorenz系统的数值解，并由此画出其相图yinit = [0.1,0.1,1];% 初始化输入y(1:3) = yinit;tstart = 0; % 时间初始值tstep = 1e-1; % 时间步长wholetimes = 1e2; % 总的循环次数steps = 1; % 每次演化的步数iteratetimes = wholetimes/steps; % 演化的次数S=output;for i=1:iteratetimes;tspan = tstart:tstep:(tstart + tstep*steps);[T,Y1] = ode45(@Lorenz_driven, tspan, y);y = Y1(size(Y1,1),:);y(1)=S(i,1);% 重新定义起始时刻tstart = tstart + tstep*steps;endfigure(1)plot3(Y1(:,1),Y1(:,2),Y1(:,3))function s=output;tstart = 0; % 时间初始值tstep = 1e-1; % 时间步长wholetimes = 1e2; % 总的循环次数% options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);tspan=tstart:tstep:wholetimes*tstep[T,Y] = ode45(@Lorenz_driven,tspan,[0.1 0.1 0.1]);s=Yfigure(3)plot3(Y(:,1),Y(:,2),Y(:,3))function dY=Lorenz_driven(t,Y);a=10;b=8/3;r=60;dY=zeros(3,1);dY=[a*(Y(2)-Y(1));-Y(1)*Y(3)+r*Y(1)-Y(2);Y(1)*Y(2)-b*Y(3)]MatLab常微分方程及常微分方程组的求解(2011-07-08 23:01:48)转载▼分类：编程之Matlab标签：杂谈最近参加了数学建模，对于老师说的Euler算法的不同步长的精度不一样，编写了一个M 函数文件来实现这个精度的比较，把函数附上：function [x,y]= Euler(varargin)%这里使用可变输出输入函数的%varargin{1}为求解常微分方程的表达式%varargin{2}为求解常微分方程的定解条件%需要给出的变量有常微分方程的范围a，b(varargin{3},varargin{4})%n为对这个区间的分割(varargin{5})%xlt写于7月8日%取得算法需要的变量，并附上容易理解的含义变量a = varargin{3};b = varargin{4};%自变量的范围n = varargin{5};%区间的分割次数h = (b - a)/n;%步长Dy = varargin{1}; %常微分方程的表达式y0 = varargin{2}; %常微分方程的定解条件表达式%首先求出所给常微分方程问题的精确解x1 = zeros(n+1,1);y1 = zeros(n+1,1);syms f1; syms x;f1 = dsolve(Dy,y0,'x');x1(1) = a;y1(1) = subs(f1,{x},{x1(1)});for i = 2:(n+1)x1(i) = x1(i-1) + h;y1(i) = double(subs(f1,{x},{x1(i)}));end%利用Euler方法求解近似数值微分解x2 = zeros(n+1,1);y2 = zeros(n+1,1);syms y;x2(1) = a;y2(1) = subs(f1,{x},{a});%获得原方程的初解for i = 2:(n+1)x2(i) = x2(i-1) + h;y2(i) = y2(i-1) + h .* double(subs(Dy(5:end),{x,y},{x2(i-1),y2(i-1)}));%特别记录Matlab中的字符串操作，提取子字符串即A（3：6）...end%返回经过Euler算法算出x与y的值x = x2;y = y2;%画图进行误差比较plot(x1,y1,'r');hold on;plot(x2,y2,'b');特此记录，以后写了新的算法再分享文 - 汉语汉字编辑词条文，wen，从玄从爻。

基于MATLAB的各类混沌系统的计算机模拟―――《混沌实验教学平台的设计与实现》初期报告物电05级1A班张丹伟20050003101摘要：本文利用数学软件MATLAB对Lorenz系统等六个重要的混沌模型进行数值计算，同时模拟出各类混沌系统的独特性质，如混沌吸引子，倍周期，初值敏感性，相图，分岔图等。

通过观察和分析上述特性，加深了我们对混沌现象的理解。

关键词：混沌；微分方程；MA TLAB；引言．混沌探秘混沌是非线性系统所独有且广泛存在的一种非周期运动形式, 其覆盖面涉及到自然科学和社会科学的几乎每一个分支。

1972年12月29日，美国麻省理工学院教授、混沌学开创人之一E.N.洛伦兹在美国科学发展学会第139次会议上发表了题为《蝴蝶效应》的论文，提出一个貌似荒谬的论断：在巴西一只蝴蝶翅膀的拍打能在美国得克萨斯州产生一个龙卷风，并由此提出了天气的不可准确预报性。

为什么会出现这种情况呢？这是混沌在作怪！“混沌”译自英语中“chaos”一词，原意是混乱、无序，在现代非线性理论中，混沌则是泛指在确定体系中出现的貌似无规则的、类随机的运动。

混沌现象是普遍的，就在我们身边，是与我们关系最密切的现象，我们就生活在混沌的海洋中。

一支燃着的香烟，在平稳的气流中缓缓升起一缕青烟，突然卷成一团团剧烈搅动的烟雾，向四方飘散；打开水龙头，先是平稳的层流，然后水花四溅，流动变的不规则，这就是湍流；一个风和日丽的夏天，突然风起云涌，来了一场暴风雨。

一面旗帜在风中飘扬，一片秋叶从树上落下，它们都在做混沌运动。

可见混沌始终围绕在我们的周围，一直与人类为伴。

一．混沌的基本概念1. 混沌: 目前尚无通用的严格的定义, 一般认为,将不是由随机性外因引起的, 而是由确定性方程(内因)直接得到的具有随机性的运动状态称为混沌。

2. 相空间: 在连续动力系统中, 用一组一阶微分方程描述运动, 以状态变量(或状态向量)为坐标轴的空间构成系统的相空间。

毕业设计（论文）原创性声明本人郑重声明：所提交的毕业设计（论文），是本人在导师指导下，独立进行研究工作所取得的成果。

除文中已注明引用的内容外，本毕业设计（论文）不包含任何其他个人或集体已经发表或撰写过的作品成果。

对本研究做出过重要贡献的个人和集体，均已在文中以明确方式标明并表示了谢意。

论文作者签名：日期：年月日摘要混沌在现代科学与工程学领域的应用十分广泛，混沌现象存在于自然界各个领域，包括通讯领域、气象学领域、生物学领域、医学诊断疾病等方面。

学习混沌理论在未来的发展过程对我们是很有帮助的。

在非线性的世界里，通过混沌理论洞察所有的非线性运动,对其进行控制和掌握。

通过非线性电路对混沌系统进行分析和理解，进而构造出符合二阶混沌系统的非线性电路和函数模型。

Duffing 方程就是典型的二阶非线性方程。

运用MATLAB/Simulink对其混沌系统进行仿真实现，验证混沌系统的基本特性。

关键词：混沌；非线性；Duffing方程; MATLAB/SimulinkABSTRACTChaos widely used in modern science and engineering and chaos phenomenon exists in various fields of nature, including the communications field, the field of meteorology, biology, medical diagnosis of diseases. Learning Chaos Theory is very helpful to us in the development of this course in the future. In a nonlinear world, insight into the chaos theory, We can control and master non-linear movement. We analyze and understand the chaotic system via nonlinear circuit, and then construct a second-order chaotic systems of nonlinear circuits and function model. Duffing equation is a typical second-order nonlinear equation. Using MATLAB/Simulink, we complete the chaotic system simulation and test the basic characteristics of chaotic systems.Key words：Chaos；nonlinear；Duffing equation；MATLAB/Simulink目录第一章绪论 (1)1.1混沌理论 (1)1.2混沌的应用 (2)第二章二阶混沌系统的仿真实现 (5)2.1混沌系统 (5)2.1.1混沌产生的数学模型 (5)2.1.2 奇异吸引子与分形 (6)2.1.3 混沌系统的特征 (7)2.1.4 研究混沌的主要方法 (8)2.2 二阶混沌系统的实现 (9)第三章二阶非线性电路仿真实现 (15)3.1 Simulink仿真 (17)3.2 MATLAB语句命令演示模拟 (19)第四章结论 (22)致谢 (25)参考文献 (26)附录A (27)第一章绪论1.1混沌理论什么是混沌？现代科学意义上是很难得出确切的定义，之所以这样是因为：到目前为止，还没有足够和统一数学定理可以将混沌理论完全表达出来，在数学理论的基础上通过混沌系统所表现出的普遍现象总结归纳出混沌的本质。

毕业设计（论文）题目一个超混沌系统在MATLAB环境下的仿真实现系（院）物理与电子科学系专业物理学班级2005级1班学生姓名XXX学号2005080119指导教师XXX职称二〇一一年六月十八日独创声明本人郑重声明：所呈交的毕业设计(论文)，是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议。

尽我所知，除文中已经注明引用的内容外，本设计（论文）不含任何其他个人或集体已经发表或撰写过的作品成果。

对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。

本声明的法律后果由本人承担。

作者签名:二〇一〇年六月一十八日毕业设计（论文）使用授权声明本人完全了解滨州学院关于收集、保存、使用毕业设计（论文）的规定。

本人愿意按照学校要求提交学位论文的印刷本和电子版，同意学校保存学位论文的印刷本和电子版，或采用影印、数字化或其它复制手段保存设计（论文）；同意学校在不以营利为目的的前提下，建立目录检索与阅览服务系统，公布设计（论文）的部分或全部内容，允许他人依法合理使用。

（保密论文在解密后遵守此规定）作者签名:二〇一〇年六月一十八日一个超混沌系统在MATLAB环境下的仿真实现摘要在混沌、超混沌理论研究成果的基础上，利用外加驱动信号方法改进一个四阶超混沌系统，通过对外界驱动信号频率的控制，实现系统的动力学特性。

对新构建的超混沌系统的特性进行了详细分析，包括验证其超混沌性质，相空间轨迹分析，Lyapunov指数谱分析等，仿真结果关键词：超混沌；lyapunov指数；EWB；超混沌电路；MATLABIA hyperchaos circuit was simulated in MATLABsimulationAbstractBased on chaotic and the hyperchaotic theory research, by using plussing a drive signal to the fourth-order hyperchaos system to improve the fourth-order hyperchaos system, through to control the external drive signal frequency, realizeing the system’ dynamic characteristics. The construction of the new characteristics of hyperchaos system are analyzed in detail, including its hyperchaos nature, path analysis, phase space Lyapunov index and bifurcation diagram analysis and simulation results show that the system characteristics. Is abundant. Using the signal frequency control and drive can completely accurate control of the entire system dynamics characteristic. Finally,design a simulated circuit, and simulat in EWB environment, through the comparison of simulation results between MATLAB and EWB, Further verify the consistency between experiment results and numerical simulation.Keyword：hyperchaos；Lyapunov exponents；bifurcation；hyperchaotic CircuitII目录引言 (1)第一章动力系统形态及其分析 (2)1.1动力系统 (2)1.1.1动力学系统的基本概念 (2)1.1.2几种常见的平衡态 (4)1.1.3吸引子和结构稳定性 (6)1.2 分岔 (7)1.2.1分岔的基本概念 (7)1.2.2非线性映射及其分岔 (8)第二章混沌系统判别方法讨论 (10)2.1混沌的特性及其判别方法 (10)2.1.1混沌的定义 (10)2.1.2混沌运动的基本特征 (11)2.2超混沌特性及其判别方法 (11)2.3混沌电路研究方法 (12)第三章一个新的超混沌系统设计及其性能分析 (12)3.1超混沌系统 (12)3.1.1一个四阶的超混沌系统 (12)3.1.2平衡点及稳定性分析 (13)3.1.3系统相空间轨迹分析 (13)3.2一个新的超混沌系统 (15)3.2.1一个新的超混沌系统的设计 (15)3.2.2系统相空间轨迹分析 (16)3.2.3李雅普诺夫指数分析 (17)3.2.4系统的电路设计和实验结果 (17)结论 (23)参考文献 (24)i谢辞 (25)ii引言混沌科学是一门新兴的学科，混沌（Chaos）是一种貌似无规则的运动，指在确定性非线性系统中，不需要附加任何随机因素亦可出现的行为（内在随机性）。

MATLAB中的计算机模拟第一节 导言计算机科学技术的迅猛发展，给许多学科带来了巨大的影响．计算机不但使问题的求解变得更加方便、快捷和精确，而且使得解决实际问题的领域更加广泛．计算机适合于解决那些规模大、难以解析化以及不确定的数学模型．例如对于一些带随机因素的复杂系统，用分析方法建模常常需要作许多简化假设，与面临的实际问题可能相差甚远，以致解答根本无法应用，这时模拟几乎成为人们的唯一选择．在历届的美国和中国大学生的数学建模竞赛（MCM）中，学生们经常用到计算机模拟方法去求解、检验等．计算机模拟(computer simulation)是建模过程中较为重要的一类方法．本章将讨论如何利用计算机技术对连续系统和离散系统进行模拟．由于计算机以重复性运算见长，故它为研究模拟方法提供了极为合适的手段．计算机模拟是一种广义的数值计算方法．通过本章的学习，你将会了解蒙特卡洛方法的思想；初步掌握对连续系统或离散系统进行模拟的方法；掌握由实际问题怎样去建立计算机模拟模型以及应用MATLAB编程语言进行计算．第二节引例第三节 随机变量的抽样第四节 连续系统的模拟第五节 离散系统的模拟第六节 范例第七节 实验第二节 引例：葡丰投针问题在用传统方法难以解决的问题中，有很大一部分可以用概率模型进行描述．由于这类模型含有不确定的随机因素，分析起来通常比确定性的模型困难．有的模型难以作定量分析，得不到解析的结果，或者是虽有解析结果，但计算代价太大以至不能使用．在这种情况下，可以考虑采用Monte Carlo 方法。

下面通过例子简单介绍Monte Carlo 方法的基本思想．Monte Carlo 方法是计算机模拟的基础，它的名字来源于世界著名的赌城——摩纳哥的蒙特卡洛，其历史起源于1777年法国科学家蒲丰提出的一种计算圆周 π的方法——随机投针法，即著名的蒲丰投针问题。

1)Monte Carlo 方法的基本思想是首先建立一个概率模型，使所求问题的解正好是该模型的参数或其他有关的特征量．然后通过模拟一统计试验，即多次随机抽样试验（确定m 和n ），统计出某事件发生的百分比．只要试验次数很大，该百分比便近似于事件发生的概率．这实际上就是概率的统计定义．利用建立的概率模型，求出要估计的参数．蒙特卡洛方法属于试验数学的一个分支．MATLAB 语言编程实现 l=1; n=1000;d=2； m=0； for k ＝l ：nx ＝unifrnd （0，d ／2）； p ＝unifrnd （0，pi ）；if )sin(15.0y x ××< m=m ＋1 elsc end end p=m/npi_m=1/p 运行，即得结果．蒙特卡洛方法适用范围很广泛，它既能求解确定性的问题，也能求解随机性的问题以及科学研究中的理论问题．例如利用蒙特卡洛方法可以近似地计算定积分，即产生数值积分问题．任意曲边梯形面积的近似计算一个古老的问题：用一堆石头测量一个水塘的面积．应该怎样做呢？测量方法如下：假定水塘位于一块面积已知的矩形农田之中．如图8．2所示．随机地向这块农田扔石头使得它们都落在农田内．被扔到农田中的石头可能溅上了水，也可能没有溅上水，估计被“溅上水的”石头量占总的石头量的百分比．试想如何利用这估计的百分比去近似计算该水塘面积？结合图8．2中的图形(1)分析，只要已知各种参数及函数（a ，b ，H ，f(x)），有以下两种方法可近似计算水塘面积．1．随机投点法1）赋初值：试验次数n=0，成功次数m=0；规定投点试验的总次数N ；2）随机选择m 个数对,1,,m i y x i i <<，其中H y b x a i i <<<<0,，置 n ＝n ＋l ； 3）判断，若是，转4，否则停止计算；N n ≤ 4）判断条件(表示一块溅水的石头)是否成立，若成立则置m=m+1，转2，否则转2；)(i i x f y < 5）计算水塘面积的近似值N m a b H S /)(×−×=．2．平均值估计法1）产生[a,b]区间的均匀随机数;,,2,1,N i x i "= 2) 计算;,,2,1),(N i x f i "= 3）计算∑=−=Ni i x f N a b S 1)()(。

基于Matlab的超混沌系统同步方法仿真
毕洪波;张玉波;王秀芳
【期刊名称】《实验室研究与探索》
【年(卷),期】2009(028)011
【摘要】为探讨超混沌同步机制,研究其在保密通信中的应用,以超混沌chen系统为研究对象,给出了超混沌系统的Matlab仿真步骤.利用复制、耦合等方法实现了耦合同步法仿真,借助于等间隔脉冲对收端系统进行控制,实现了脉冲同步法仿真.结果表明:超混沌耦合同步中,耦合系数大,则同步时间短;脉冲同步中,脉冲周期小,则同步时间短.
【总页数】4页(P39-42)
【作者】毕洪波;张玉波;王秀芳
【作者单位】大庆石油学院电气信息工程学院,黑龙江,大庆,163318;大庆石油学院电气信息工程学院,黑龙江,大庆,163318;大庆石油学院电气信息工程学院,黑龙江,大庆,163318
【正文语种】中文
【中图分类】TP391
【相关文献】
1.基于非线性反馈控制的超混沌系统同步方法 [J], 赵辽英;赵光宙;厉小润
2.基于状态观测器的超混沌系统高精度同步方法 [J], 张学义;李殿璞;陈实如;王文武
3.时变参数超混沌系统的同步方法与实现 [J], 王志坚
4.超混沌系统的多变量驱动误差反馈控制同步方法 [J], 孙克辉;张泰山
5.Rossler超混沌系统的非线性同步方法研究 [J], 朱夜明;乔宗敏
因版权原因，仅展示原文概要，查看原文内容请购买。