使用线性连接图、蛛网图或分枝与混沌图对参数 a 进行讨论与观察【数学实验】【matlab】

选择下列三个函数中的至少一个函数进行迭代，分析迭代产生的序列的收敛性。

使用线性连接图、蛛网图或分枝与混沌图对参数 a 进行讨论与观察，会得到什么结论？

本次选用研究 其他的如出一辙。

• 算法与编程

1.原始函数fun.m的编写。

function [output] = fun(a,x)

%主函数f(x)=a-(x-a'(1/2))^2

output=a-(x-sqrt(a)).^2;

end

2.线性连接图line_connect_pic.m的编写。

clf;

clc;

clear;

%绘出线性连接图。

a=1:1:4;

x=[];

%通过内外两次循环，在同一界面显示4次不同参数对应的线性连接图

for i=1:4  

    x(i,1)=0.5;

    for j=2:25

        x(i,j)=a(i)-(x(i,j-1)-sqrt(a(i)))^2;

    end

n=1:25;

subplot(2,2,i),plot(n,x(i,:)),title(['a=',num2str(a(i)),' x0=0.5时的线性连接图。'])

end

3.蛛网图图spider_pic.m的编写。

clf;

clc;

clear;

%循环做出a和x0不同值时对应的蛛网图。

a=linspace(1.5,4,4);

x=-1:0.25:5;

xx=[];

%通过内外两次循环，在同一界面显示4次不同参数或初始值对应的蛛网图

for i=1:4

    y=fun(a(i),x);

    %前两次a值时，x0取0.1；后两次a值时x0取0.3.

    if a(i)<=3

        xx(i,1)=0.1;

    else

        xx(i,1)=0.3;

    end

    %绘出函数图像，设置标题和横轴标签

    subplot(2,2,i),plot(x,y),xlabel('x'),title(['a=',num2str(a(i)),' x0=',num2str(xx(i,1))]);

    hold on;

    %绘出y=x的图像和y=0的图像

    subplot(2,2,i),plot(x,x)

    subplot(2,2,i),plot(x,zeros(size(x)));

    %迭代出折线

    for j=2:20

        xx(i,j)=a(i)-(xx(i,j-1)-sqrt(a(i)))^2;

        subplot(2,2,i),plot([xx(i,j-1),xx(i,j-1)],[xx(i,j-1),xx(i,j)]);

        subplot(2,2,i),plot([xx(i,j-1),xx(i,j)],[xx(i,j),xx(i,j)]);

    end

    hold off

end

4.含参数变量函数进行迭代的可调用程序iter.m的编写

function root = iter(x,a)

%root函数 用来产生含参数变量函数f(x,a)进行迭代的可调用程序

%   输入自变量x和参数a.

for i=2:100

     x(i)=a-(x(i-1)-sqrt(a))^2;

end

%返回100个迭代序列组成的数组

root=x;

5.分枝与混沌图的生成文件fgbm.m的编写。

clf;

%绘出费根鲍姆图

x=[];

x(1)=0.5;

hold on

for a=0:0.01:4

    %调用root函数生成迭代序列数组

    root=iter(x,a);

    %用.循环做出不同的参数a对应的序列图像

    plot(a.*ones(size(root(51:100))),root(51:100),'.');

end

xlabel('参数a');

ylabel('迭代序列（51-100）');