MATLAB程序设计实验题答案

MATLAB程序设计实验

解方程和最优化问题求解

姓名：龙刚 学号：1404110427 班级：电信1104班 一．实验目的

1. 掌握线性方程组的数值求解方法。 2. 掌握常微分方程的数值求解方法。

3. 掌握非线性方程以及最优化问题的求解方法。

二．实验内容

1. 对下列线性病态方程组：

?1/21/31/4??x1??0.95??1/31/41/5??x???0.67? ???2?????1/41/51/6????x3????0.52??a) 求方程组的解。

b) 将方程组右边的向量元素b3（0.52）修改为0.53，再次求方程组的解。比较b3的

变化和解的相对变化。

c) 求系数矩阵的条件数，并分析条件数对线性方程组求解的影响。

%1.a

A=[1/2 1/3 1/4;1/3 1/4 1/5;1/4 1/5 1/6]; B=[0.95 0.67 0.52]'; disp('1.a方程组的解为：') xa=A\\B %1.b

C=[0.95 0.67 0.53]'; disp('1.b方程组的解为：') xb=A\\C chazhi=xa-xb %1.c t1=cond(

MATLAB程序设计实验

解方程和最优化问题求解

姓名：龙刚 学号：1404110427 班级：电信1104班 一．实验目的

1. 掌握线性方程组的数值求解方法。 2. 掌握常微分方程的数值求解方法。

3. 掌握非线性方程以及最优化问题的求解方法。

二．实验内容

1. 对下列线性病态方程组：

?1/21/31/4??x1??0.95??1/31/41/5??x???0.67? ???2?????1/41/51/6????x3????0.52??a) 求方程组的解。

b) 将方程组右边的向量元素b3（0.52）修改为0.53，再次求方程组的解。比较b3的

变化和解的相对变化。

c) 求系数矩阵的条件数，并分析条件数对线性方程组求解的影响。

%1.a

A=[1/2 1/3 1/4;1/3 1/4 1/5;1/4 1/5 1/6]; B=[0.95 0.67 0.52]'; disp('1.a方程组的解为：') xa=A\\B %1.b

C=[0.95 0.67 0.53]'; disp('1.b方程组的解为：') xb=A\\C chazhi=xa-xb %1.c t1=cond(

系立项建设 实验指导书

MATLAB程序设计与应用实验

主 编：王 春、段守付

物理与电子信息工程系

目 录

前 言 .................................................. 2 实验一 熟悉MATLAB 环境................................... 3 实验二 数值数组创建、应用及可视化 ........................ 7 实验三 字符串数组的使用、简单脚本文件和函数的编写 ....... 11 实验四 数据可视化方法 ................................... 18 实验五 简单数据分析 ..................................... 24 实验六 MATLAB 仿真初步 .................................. 28 实验七 电路模型仿真应用 ................................. 31 实验八 简单图形用户界面设计 ............................. 33 附录：MATLAB 主要函数指令表（按字母

中南大学MATLAB实验

报告五

2014年6月2日

MATLAB程序设计实验5

数据分析和多项式计算

一．实验目的

1. 掌握数据统计和分析的方法。

2. 掌握数值插值与曲线拟合的方法和应用。 3. 掌握快速傅立叶变换的应用方法。 4. 掌握多项式的常用运算

二．实验内容

1. 利用randn函数生成符合正态分布的10x5随机矩阵，进行以下操作：

a) 矩阵中各列元素的均值和标准方差。 解：程序及结果如下： A=randn(10,5); pjz=mean (A) fc=std(A,0,1) pjz =

-0.1199 -0.1975 -0.2014 -0.6891 -0.2486 fc =

0.5407 0.9103 1.0901 0.8234 0.8529 b) 矩阵的最大元素和最小元素。 解：程序及结果如下： A=randn(10,5); zdz=max(max(A)) zxz=min(min(A)) zdz =

1.5374 zxz =

-2.5200

c) 求矩阵每行元素的和以及全部元素之和。 解：程序及结

电子信息工程系实验报告

课程名称： MATLAB语言 成 绩： 实验项目名称： 实验3 MATLAB程序设计 实验时间：2013-4-6 指导教师（签名）： 1. 实验目的 学习函数的编制，掌握MATLAB的编程应用。通过不同的程序结构和不同的实际编程问题，掌握MATLAB的编程方法。 2.实验环境

硬件：PC机，基本配置CPU PII以上，内存4G以上； 软件：Matlab 版本7.1

3.实验原理

MATLAB是一个专门的数学软件，特别是对矩阵进行操作，可以方便快捷的进行矩阵运算，从而简化了人们的脑力，同时能更加准确地计算结果。

4.实验内容及结果分析

4.4编写一个程序，计算出坐标系中用户指定两点(X1,Y1)和(X2,Y2)之间的距

离。要求有输入、输出及其相关提示。 编写程序如下： clc; clear;

X1=input('please input X1=');

Y1=input('Y1='); X2=input('X2=');

Y2=input('Y2=');

s=sqrt((X1-X2)^2+(Y1-Y2)^2);

fprintf('两点距

实验五 MATLAB程序设计

一、实验目的

1、熟练掌握MATLAB的程序流程控制结构；

2、熟练掌握M文件的结构；

3、熟练掌握函数调用和参数传递；

二、实验内容

Fibonacci数列的个元素为：1、1、2、3、5、8、……，满足以下关系F1=1，F2=1，Fn=Fn-1+Fn-2

用M函数文件实现，数列的元素个数为输入变量。

（1） 按M函数文件格式创建文件开头

function f=shiyan4(n) %SHIYAN4 Fibonacci % Fibonacci数列 %n 元素个数

%f 构成 Fibonacci数列向量 %copyright 2009-04-09

（2） 用while循环实现程序功能

f(1)=1;f(2)=1; i=2;

while i<=n

f(i+1)=f(i-1)+f(i); i=i+1;

end

在命令窗口输入调用命令，调用函数结果如下： >> f=shiyan4(10) f =

1 1 2 3 5 8 13 21 34

（3） 用for循环实现

f(1)=1;f(2)=1

2014实验课

第二周

1. 编程求两整数的和（要求 用函数调用方式实现）。

2. 编写程序，打印图案，图案可以自行设计，例如：画一条线，直角三角形等。（要求：输入数字，得出相应大小的图案。并用函数调用方式实现）。 3. 修改第二题中函数，使得输出的图案整体向右偏移若干列。

4. 输出200以内的所有素数，要求定义一个yes函数，用来判别一个数是不是素数。

5. 计算1+（1+2）+（1+2+3）+…+(1+2+3+…+n)。要求（1)n由键盘输入 （2）每项的计算由函数完成。

第三周

1. 定义如下分段函数sign(x)并在main函数中调用它。 1 (x>0) Sign(x)= 0 (x=0) -1 (x<0)

2. 输入一批正整数（以零或负数为结束标志），求其中的奇数和。要求定义和调用函数

even(n)判断数的奇偶性，当n为偶数时返回1，否则返回0。

3. 求s= 13+23+33+…+n3。要求：n的值在main函数中输入，结果也在main函数中输出。

其它功能用函数完成（提示：函数嵌套调用）。

4. 输入精度e，输入x，用下列公式求cos x的近似值，精确到最后一项

2014实验课

第二周

1. 编程求两整数的和（要求 用函数调用方式实现）。

2. 编写程序，打印图案，图案可以自行设计，例如：画一条线，直角三角形等。（要求：输入数字，得出相应大小的图案。并用函数调用方式实现）。 3. 修改第二题中函数，使得输出的图案整体向右偏移若干列。

4. 输出200以内的所有素数，要求定义一个yes函数，用来判别一个数是不是素数。

5. 计算1+（1+2）+（1+2+3）+…+(1+2+3+…+n)。要求（1)n由键盘输入 （2）每项的计算由函数完成。

第三周

1. 定义如下分段函数sign(x)并在main函数中调用它。 1 (x>0) Sign(x)= 0 (x=0) -1 (x<0)

2. 输入一批正整数（以零或负数为结束标志），求其中的奇数和。要求定义和调用函数

even(n)判断数的奇偶性，当n为偶数时返回1，否则返回0。

3. 求s= 13+23+33+…+n3。要求：n的值在main函数中输入，结果也在main函数中输出。

其它功能用函数完成（提示：函数嵌套调用）。

4. 输入精度e，输入x，用下列公式求cos x的近似值，精确到最后一项

MATLAB是一个运算性质的软件

实验三（2）、MATLAB程序设计

一、 实验目的

1、 掌握建立和执行M文件的方法

2、 掌握利用for语句和while语句实现选择结构的方法

二、 实验内容及步骤

1、 分别使用for循环语句和while循环语句找出最小的n值，

使得n!>10100,并求出n!。要求显示的结果为

n=70

n!=1.197857e+100>1e100

或

70!=1.197857e+100>1e100

用IF语句执行如下：

s=1;

for n=1:100

s=s*n;

if(s>10^100)

break;

end

end

str1=['n=',num2str(n)];

str2=['n!=',num2str(s)];

disp(str1)

disp(str2)

用WHILE语句执行如下：

n=1;

s=1;

while(s<=10^100)

s=s*n;

n=n+1;

end

str1=['n=',num2str(n-1)];

str2=['n!=',num2str(s)];

disp(str1)

disp(str2)

MATLAB是一个运算性质的软件

2、 已知 f1=1,n=1

f2=0,n=2

f3=1,n=3

fn=fn-1-2fn-2+fn-3,n>3

实验四 循环结构程序设计 1.根据??621123?21?21???21n，2x =

3.1406 求?的近似值。当n分别取100、1000、10000时，记过是多少？

n=100

程序设计：

n=100;

i=1:n;

f=1./i.^2;

x=sum(f);

y=sqrt(6*x)

运行结果：

x =

3.1321 n=1000

程序设计：

n=1000;

i=1:n;

f=1./i.^2;

x=sum(f);

y=sqrt(6*x)

运行结果：

n=10000

程序设计：

n=1000;

i=1:n;

f=1./i.^2;

x=sum(f);

y=sqrt(6*x)

运行结果：

x =

3.1415

2.根据y=1+1113?5???2n-1,求： （1）y<3时的最大n值。

（2）与（1）的n值对应的y值。程序设计：

y=1;n=1;

while(y<3)

n=n+2;

y=y+1/n;

end

1

y=y-1/n n=(n+1)/2 运行结果： y = 2.9944 n =

57 3.

考虑以下迭代公式：

10e-5&n<=500) x=a/(b+x); n=n+1; end x n

运行结果： x = 0.5616 n =

6

（2）如果迭代过程收敛于r，那

xn?1?a b?