MATLAB实验报告第八次.docVIP

MATLAB学实验报告 姓名： 李犁 学号： 1304120724 指导老师：易昆南 专业班级：统计1203 目录 1循环矩阵的性质 …………………………3 2航线连接问题 ……………………………6 3 生物群体模型…………………………….8 4学生成绩管理系统……………………….11 5调和级数问题……………………………14 6 3N+1问题 ……………………….….......17 7 航程与起点……………………………...20 8 ε-δ图的动画绘制……………….……27 9求π………………………………………32 10求函数f的一阶泰勒展式………………33 11蛛网模型的动画显示…………………..36 12音乐制作……..…………………………40 13用matlab做插值计算…………………49 14以g(x)为例作三种插值的比较………...53 15 函数的拟合………………...60 16 函数的拟合…...62 17 2011年数模竞赛A题的作图…………64 18混沌现象……….……………………….77 19分形图案的绘制………………………..79 20古典概型与随机问题的模拟问题……..84 21抽取小球问题……...…………………...86 22 追赶火车问题………………………….89 23 Galton钉板模型………………………..91 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 循环矩阵的性质 评 分 设计（实习）目的： 探究可逆循环矩阵的性质，运用matlab，初步了解程序的编写，再根据结果自行进行证明。 2、实验内容： 编写一个可逆循环矩阵程序，实现矩阵的随机性，在用随机产生的矩阵进行研究，从而研究任意一个可逆循环矩阵的性质，包括矩阵的和、积、逆的性质。 3.详细设计： 此程序只有一个主函数，储存在一个M文件中，实验时调用函数完成。 主函数（matrixrand） function output=matrixrand(n,a,b) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here r=fix(a+(b-a)*rand(1,n)); %r向量是一组随机数 N=zeros(n); %把矩阵N初始为一个n阶零矩阵 while rank(N)~=n %判定是否可逆 N(1,:)=r; %矩阵第一行是随机向量r for i=2:n %以下为循环矩阵的实现，从第二行到第n行 r=[r(1,n) r(1,1:n-1)]; N(i,:)=r; end end output=N; %输出可逆循环矩阵 end 4：实验结果： 在命令窗口输入 A=matrixrand(5,1,5) 输出 A = 1 2 3 4 4 4 1 2 3 4 4 4 1 2 3 3 4 4 1 2 2 3 4 4 1 1.在命令窗口输入 B=inv(A) 输出 B = -0.2412 -0.0776 -0.0048 0.2497 0.1406 0.1406 -0.2412 -0.0776 -0.0048 0.2497 0.2497 0.1406 -0.2412 -0.0776 -0.0048 -0.0048 0.2497 0.1406 -0.2412 -0.0776 -0.0776 -0.0048 0.2497 0.1406 -0.2412 2.在命令窗口输入 C=matrixrand(5,1,5); D=matrixrand(5,1,5); F=C+D 输出 F = 7 4 4 4 2 2 7 4 4 4 4 2 7 4 4 4 4 2 7 4 4 4 4 2 7 3.在命令窗口输入 G=C*D 输出 G = 24 2