matlab河北大学,河北大学matlab实验一报告.doc

河北大学matlab实验一报告河北大学matlab实验一报告

实验一 MATLAB 基本操作

一、 实验目的

1． 熟悉MATLAB 的使用方法及特点；学会建立MATLAB 搜索路径；熟悉MATLAB 工作空间、MATLAB 集成环境、命令窗口；掌握MATLAB 的通用命令、管理命令和函数、管理变量和工作空间的使用方法；

2． 掌握MATLAB 基本操作及矩阵基础知识，包括：输入矩阵、矩阵的转置、矩阵元素求和、矩阵下际、矩阵连接、矩阵行列删除、矩阵产生和操作、逻辑和关系运算、操作符和特殊字符、基本矩阵和矩阵操作、基本矩阵和阵列；掌握特殊变量和常数；掌握基本数学函数。

二、实验内容

1.1 建立工作目录MYBIN和MYDATA。并在MYBIN和BIN(MATLAB系统自动生成)中分别建立test.m文件，然后在MATLAB中键入test，这时检查执行的是哪一个test.m文件？如改变搜索路径，结果又会如何？如果在MATLAB环境中建立一个名为test的变量，再次执行结果会如何？text

text=3

text

text

Warning: Function D:\Program Files\MATLAB71\mybin\text.m has the same name as a MATLAB builtin. We suggest you rename the function to avoid a potential name conflict.

b =

1

>> text

Warning: Function D:\Program Files\MATLAB71\bin\text.m has the same name as a MATLAB builtin. We suggest you rename the function to avoid a potential name conflict.

a =

2

>> text=3

text =

3

>> text

text =

3

1． 利用基本矩阵产生 3x3 和15x8 的单位阵，全1 阵，全0 阵，均匀分布的随机阵([-1，1]之间)，正态分布随机阵(方差4，均值1)

a=eye(3)

a=eye(15,8)

a=zeros(3)

a=zeros(15,8)

a=ones(3)

a=ones(15,8)

a=1-2*rand(3)

a=1-2*rand(15,8)

a=2*randn(3)+1

a=2*randn(15,8)+1

2．利用diag()函数和rot90()产生下列矩阵：

然后求解a 阵的逆矩阵aa 及b 阵的特征值和对应特征向量，并利用reshape 将aa 阵变换成行向量。

v1=[8 -7 2];

v2=[5 3];

a =diag(v1)+diag(v2,1)

a=rot90(a)

a=rot90(a)

a=rot90(a)

v1=[2 5 8];

v2=[4 0 7];

v3=diag(v2)

v3=rot90(v3)

b=diag(v1)+v3

3． 产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。

a=5-round(10*10*rand(50,2))/10

4．编程实现当α∈[-π，π]，间隔为1o 时，求解正弦和余弦的值，并利用plot()函数绘制正弦，余弦曲线。

>> x=-pi:pi/180:pi;

>> y=sin(x);

>> plot(x,y)

>> y=cos(x);

>> plot(x,y)

5．利用rand 函数产生(0,1)间均匀分布的10x10 随机矩阵a,然后统计a 中大于等于0.6 的元素个数。

a=rand(10,10)

b=a>=0.6

c=b(:)

d=sum(c)

6．利用randn 函数产生均值为0，方差为1 的10x10 正态分布随机阵，然后统计其

中大于-0.5，小于0.5 的元素个数。

a=randn(10,10)

b=a>-0.5&a<0.5

c=b(:)

d=sum(c)

7．编程实现下表功能：

5

a b b<=0.5 b>0.5

a<1 显示‘case 1’ 显示‘case 2’

a>=1 显示‘case 3’ 显示‘error’

8．有一矩阵a,找出矩阵中其值大于1 的元素，并将他们重新排列成列向量b。

a=2*rand(3,3);

[i,j]=find(a>1)

b=a(find(a>1))

9．在一保定市区9 月份平均气温变化测量矩阵temp_Baoding_sep 中(48x30)，存在

有奇异值(大于42o C，小于0o C),编程实现删除奇异值所在的行。

a=-1+44*rand(48,30)

[i,j]=find(a>42|a<