电科132王健的matlab课程资料.docVIP

序号: 学号: 课 程 报 告 课程报告名称： Matlab应用软件 学 生 姓 名： 王健 专 业 班 级： 电科 132 学 院（系）： 信息科学与工程学院 数理学院 校内指导教师： 姜忠义 专业技术职务： 讲师 2015 年8月 1日 —— 2015 年 12月 28日 目 录 目 录 2 1． 第一部分 3 2． 第二部分 4 3． 第三部分 6 4． 第四部分 7 5． 第五部分 7 6． 第六部分 8 7． 第七部分 9 8． 第八部分 10 9． 第九部分 11 10． 心得体会 15 第一部分 用matlab画彗星运动 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%使当前轴及图形不在具备被刷新的性质，关闭在此基础上的再画图 运行结果： 第二部分 图形处理 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);%显示网格线在二行二列的第二行第二列的坐标图上 运行结果： 第三部分 Matlab 函数指令 f1=sin(x^2) f2=3*a^2+5*a+3 a1=2; a2=100 diff(f1,x)%对f1函数关于x一次求导 isstr(a1)%判断a1是不是字符串类型的 subs(f1,x,b)%将f1函数中的x用b赋值 int(f2,a)%将f2关于a的积分 double(a1)%将a1转为双精度浮点类型 运行结果： 第四部分 function a=b)%定义一个功能函数判定一个数是不是素数 for i=2:(b-1) a=mod(b,i); if a==0 a=不是素数; break; else %continue; a=是素数; end end end 运行结果： 第五部分 function x=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 end end end 运行结果： 第六部分 function c=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 end end 运行结果： 第七部分 function y=n)%定义一个非不拉列函数 if n==1||n==2 y=1; else y= n-1)+n-2); end 运行结果： 第八部分 文件的输写a=0:1:10; f=[a;3*a+2]; fid=fopen(b.txt,w);%将指定文件写入b.txt fprintf(fid,%d %d\n,f);%将f