电科132王健的matlab课程报告.doc

序号: 学号: 13461220 课 程 报 告课程报告名称： Matlab应用软件 学 生 姓 名： 王健 专 业 班 级： 电科 132 学 院（系）： 信息科学与工程学院 数理学院 校内指导教师： 姜忠义 专业技术职务： 讲师 2015 年8月 1日 2015 年 12月 28日 目 录目 录21第一部分32第二部分43第三部分64第四部分75第五部分76第六部分87第七部分98第八部分109第九部分1110心得体会151 第一部分用matlab画彗星运动t = (1/2:2/5:9/2)*pi;%t函数x = cos(t); x = x(1:2:end-1) x(1);%t的横坐标函数y = sin(t); y = y(1:2:end-1) y(1);%t的纵坐标函数plot(x,y,r);%用红色线显示图像axis equal;%坐标轴的长度单位设成相等axis off%清除坐标刻度N=10;T=N*pi*(0:0.0005:1);%T函数X=sin(T);Y=cos(T);%编译X、Y关于T的函数plot(X,Y,g);%用蓝色线在原图象上显示动态图像hold on;%使当前轴及图形保持而不被刷新，再次基础上再画图comet(X,Y,0.02);%质点运动轨迹动画的绘制hold off%使当前轴及图形不在具备被刷新的性质，关闭在此基础上的再画图运行结果：2 第二部分图形处理a=-3:0.2:3;%定义横坐标b=1/sqrt(2*pi)*exp(-1/2*a.2);%定义纵坐标subplot(2,2,1),plot(a,b,r-*)%用点画线符号*显示在二行二列的第一行第一列的坐标图上subplot(2,2,2),plot(a,b,r*)%用细实线符号*显示在二行二列的第一行第二列的坐标图上X=-1:0.1:1;Y=X;X1=X.2;Y1=Y.2;x=ones(3,1);x=ones(length(X),1);%生成X长度行1列且所有元素均为1的矩阵y=ones(1,length(Y);%生成1行Y长度列且所有元素均为1的矩阵X1=x*X1;Y1=Y1*y;Z=X1+Y1;subplot(2,2,3),mesh(X,Y,Z);%显示网格线在二行二列的第二行第一列的坐标图上subplot(2,2,4),surf(X,Y,Z);%显示网格线在二行二列的第二行第二列的坐标图上运行结果：3 第三部分Matlab 函数指令syms x a%定义符号函数f1=sin(x2)f2=3*a2+5*a+3a1=2;a2=100diff(f1,x)%对f1函数关于x一次求导isstr(a1)%判断a1是不是字符串类型的subs(f1,x,b)%将f1函数中的x用b赋值int(f2,a)%将f2关于a的积分double(a1)%将a1转为双精度浮点类型运行结果：4 第四部分判定一个数是不是素数function a=b13461220(b)%定义一个功能函数判定一个数是不是素数for i=2:(b-1) a=mod(b,i); if a=0 a=不是素数; break; else %continue; a=是素数; endendend运行结果：5 第五部分将行列式化成最简行列function x=b13461220(x)%定义一个功能函数将行列式化成最简行列式a b=size(x);for k=1:b-1 for i=a:-1:k+1 for j=b:-1:k x(i,j)=x(i,j)-x(k,j)*(x(i,k)/x(k,k); end endendend运行结果：6 第六部分二分法计算函数function c=b13461222(f,a,b) %二分法计算函数while (b-a)1e-7 c=(a+b)/2; if (f(c)*f(b)0 b=c; else if (f(c)*f(b)0 a=c; else c; break; end endend运行结果：7 第七部分非不拉列函数function y=b13461220(n)%定义一个非不拉列函数if n=1|n=2 y=1;else y= b13461220(n-1)+b13461220(n-2);end运行结果：8 第八部分文件的输写a=0:1:10;f=a;3*a+2;fid=fopen(b.txt,w);%将指定文件写入b.txtfprintf(fid,%d %dn,f);%将f值输出到指定文件c,count=fscanf(fid,%d %d ,5 inf);%c,count读取5 inf的数值fprintf(1,%d %dn,c);%将c值输出fseek(fid,0,eof);%定文件第零个元素位置d=ftell(fid);%得到文件数据fprintf(1,File Size=%dn,d);%将d值输出frewind(fid);%位置指针移到文件首部d=ftell(fid);%得到文件数据fprintf(1,File Position=%dn,d);%将d值输出fclose(fid);%关闭文件运行结果：9 第九部分课程训练（数据类型函数库）clear student%清除student变数=马云;%加入name字段student.id=阿里巴巴;%加入id字段student.scores=500,900;%加入scores字段student%显示结果（student代表一个结构数组的第一个元素）运行结果：clear student%清除student变数=马云;%加入name字段student.id=阿里巴巴;%加入id字段student.scores=500,900;%加入scores字段student(2).name=马化腾;%加入name字段student(2).id=腾讯;%加入id字段student(2).scores=300,800;%加入scores字段student%显示结果（student代表一个结构数组的第一个元素）%此结构太过复杂，matlab并不将所有字段印出运行结果：clear student%清除student变数student=struct(name,叶良成,叶鹏飞,scores,10 5,90 99);student(1)%显示student(1)student(2)%显示student(2)运行结果：clear student%清除student变数student(1)=struct(name,叶良成,scores,10 5);%输入student1数据student(2)=struct(name,叶鹏飞,scores,90 99);%输入student2数据values=struct2cell(student)%传回的异质数组维度扩大了运行结果：clear student%清除student变数=马云;%加入name字段student.id=阿里巴巴;%加入id字段student.scores=500,900;%加入scores字段student(2).name=马化腾;%加入name字段student(2).id=腾讯;%加入id字段student(2).scores=300,800;%加入scores字段student(2).name=王石;student(1)%显示student(1)student(2)%显示student(2)运行结果：clear student%清除student变数student(1)=struct(name,叶良成,scores,10 5);%输入student1数据student(2)=struct(name,叶鹏飞,scores,90 99);%输入student2数据student(2).age=21;%加入新数据student(1)%显示student(1)student(2)%显示student(2)运行结果：clear student%清除student变数student=struct(name,叶良成,scores,10 5);student=rmfield(student,scores);%删除scores字段student%显示student运行结果：clear student%清除student变数student(1)=struct(name,叶良成,scores,10 5);%输入student1数据student(2)=struct(name,叶鹏飞,scores,90 99);%输入student2数据A=cat(2,student.scores)%A=cat(dim,structureField);dim=1代表上下并排以改变横列的维度；dim=2代表左右并排以改变直行的维度average=mean(cat(2,student.scores)%计算平均值运行结果：clear student%清除student变数student(1)=struct(name,叶良成,scores,10);%输入student1数据student(2)=struct(name,叶鹏飞,scores,90);%输入student2数据student(3)=struct(name,龙傲天,scores,11);%输入student3数据for i=1:length(student)%打印每个学生的名字 fprintf(student %g:%s scores %gn,i,student(i).name,student(i).scores);end运行结果：10 心得体会一学期的MATLAB学习已经结束，虽然时间很短，但我还是从中学到了很多，对MATLAB有了一定的了