完整Matlab软件应用与开发实验指导书(柏宏斌).docVIP

《Matlab软件应用与开发》实验指导书 实验一 Matlab 操作基础及矩阵运算 (一)实验目的： （1）熟悉MATLAB软件中关于矩阵运算的各种命令； （2）学会运用MATLAB软件自定义函数，并求出函数值； （3）学会在MATLAB环境下编写函数。 (二)预备知识： 计算机操作基础，Matlab软件操作基础。 (三)实验内容： 3.1 Matlab软件操作（具体内容见专题讲座PPT） 3.1.1 建立随机矩阵： (1) 在区间[20,50]内均匀分布的5阶随机矩阵。 (2) 均值为0.6、方差为0.1的5阶正态分布随机矩阵。 命令如下： x=20+(50-20)*rand(5) y=0.6+sqrt(0.1)*randn(5) 3.1.2 建立一个字符串向量，然后对该向量做如下处理：(1) 取第1～5个字符组成的子字符串。(2) 将字符串倒过来重新排列。(3) 将字符串中的小写字母变成相应的大写字母，其余字符不变。(4) 统计字符串中小写字母的个数。 在Editor里输入如下源代码，并保存为m文件char_1.m，再点击运行。 ch=ABc123d4e56Fg9 subch=ch(1:5) %取子字符串 revch=ch(end:-1:1) %将字符串倒排 k=find(ch=ach=z) %找小写字母的位置 ch(k)=ch(k)-(‘z-‘a) %将小写字母变成相应的大写字母 char(ch) length(k) %统计小写字母的个数 3.1.3 编写读取data.txt文件数据的程序。 首先生成data.txt文本文件保存在Matlab的工作目录，文件具体内容如下： 1 2 3 4 5 6 7 8 9 10 11 12 方法一：编写以字符串形式读数据的命令文件scan_data.m，源代码如下： fid=fopen(data.txt,r) [A,COUNT]= fscanf (fid,%s) fclose(fid) 方法二：编写按行读取文本文件数据的命令文件flread.m，源代码如下： fid=fopen(data.txt,rt); if fid0 sprintf(failure for reading file data.txt) return end row=0; while ~feof(fid) str=fgets(fid); row=row+1; ftell(fid) if row1 m(row,:)=str2num(str);%某些行少了数据会报错,该行可改写为。 else m=str2num(str); end end fclose(fid); m 3.2矩证运算与M文件使用 线性代数中的矩阵运算； 本实验所用MATLAB命令 矩阵输入格式：A=[a11 a12;a21 a22];b=初始值:步长:终值 求A的转置：A’ 求A加B: A+B 求A减B: A-B 求A乘B: A*B 求A的行列式: det(A) 求A的逆: inv(A) 求A的秩: rank(A) 建立函数文件：fun.m，建立主函数文件main_fun.m (四)实验操作 输入矩阵A,B,b; ，， C1=A’，C2=A+B,C3=A-B,C4=AB D1=|A|,D2=|B| E1=R(A),E2=R(B) 定义函数 f(x1,x2)=100(x2-x12)2+(1-x1)2并计算f(1,2) (五)、源代码 A=[21 24 85 4;2 35 15 34; 21 35 31 54;21 72 15 52] B=[12 45 1 24;18 72 53 35; 48 1 15 35;46 56 25 23] b=21:2:27 C1=A,C2=A+B,C3=A-B,C4=A*B D1=det(A),D2=det(B) E1=rank(A),E2=rank(B); 〉〉[F H]=eig(A)，G=inv(rand(3)) 第一步：建立函数文件：fun.m function f=fun(x) f=100*(x(2)-x(1)^2)^2+(1-x(1))^2; 第二步：建立主程序文件main_fun.m x=[1 2] fun(x) 第三步：在Command Window输入main_fun回车运行主程序。 (六)结果分析 对求解的不同方法与速度进行分析，检验结果地正确性。 (七)总结 （1）实验中出现的问题，自己的解决方案。 （2）对自己的实验成果评价和实验心得 实验