MATLAB).doc

MATLAB一function c=f2(a,b)k=input(enter 0 or 1 please:);a=input(enter matrix a:);b=input(enter matrix b:);c=zeros(size(a,1),size(b,2);if k=0 if size(a,2)=size(b,1) c=a*b; else disp(a与b无法相乘) endelse if size(a)=size(b) for i=1:size(a,1) for j=1:size(a,2) c(i,j)=a(i,j)*b(i,j); end end else disp(只有同型矩阵才能点乘) endend二t=0:0.01:2*pi;%在同一坐标中，可以绘制3个同心圆，并加坐标控制。x=exp(i*t);y=x;2*x;3*x;plot(y)grid on; %加网格线box on; %加坐标边框 axis equal %坐标轴采用等刻度三x,y=meshgrid(-8:0.5:8);%在xy平面内选择区域-8,8,绘制4种三围曲面图。z=sin(sqrt(x,2+y,2),/sqrt(x,2+y,2+eps);subplot(2,2,1);mesh(x,y,z);title(mesh(x,y,z)subplot(2,2,2);meshc(x,y,z);title(meshc(x,y,z)subplot(2,2,3);meshz(x,y,z)title(meshz(x,y,z)subplot(2,2,4);surf(x,y,z);title(surf(x,y,z)四subplot(2,2,1);%隐函数绘图。ezplot(x2+y2-9);axis equalsubplot(2,2,2);ezplot(x3+y3-5*x*y+1/5)subplot(2,2,3);ezplot(cos(tan(pi*x),0,1)subplot(2,2,4);ezplot(8*cos(t),4*sqrt(2)*sin(t),0,2*pi)五subplot(2,2,1);%隐函数绘图。ezplot(x2+y2-9);axis equalsubplot(2,2,2);ezplot(x3+y3-5*x*y+1/5)subplot(2,2,3);ezplot(cos(tan(pi*x),0,1)subplot(2,2,4);ezplot(8*cos(t),4*sqrt(2)*sin(t),0,2*p六%3种图形者色方式的效果展示。x,y,z=sphere(20);colormap(copper);subplot(1,3,1);surf(x,y,z);axis equalsubplot(1,3,2);surf(x,y,z);shading flat;axis equalsubplot(1,3,3);surf(x,y,z);shading interp;axis equal七%绘制peaks函数曲面并且将它绕z轴旋转。x,y,z=peaks(30);surf(x,y,z)axis(-3,3,-3,3-10,10)axis off;shading interp;colormap(hot);m=movieim(20); %建立一个20列大矩阵for i=1:20 view(-37,5+24*(i-1),30) %改变视点 m(:,i)=getframe; %将图形保存到m矩阵endmovie(m,2); %播放画面2次八%光照处理后的球面。x,y,z=phere(20);subplot(1,2,1);surf(x,y,z);axis equal;light(posi,0,1,1);shading interp;hold on;plot3(0,1,1,p);text(0,1,1,light);subplot(1,2,2);surf(x,y,z);axis equal;light(posi,1,0,1);shading interp;hold on;plot3(1,0,1,p);text(1,0,1light);九%取第15个字符组成的子字符串。ch1=input(Enter your charaters:,s);one=ch1(1:5) b=length(ch1);ch2=;c=b; for k=1:b ch2=ch2,ch1(c); c=c-1; end十%将字符串倒过来重新排列。ch3=; for i=1:b if (a=ch1(i)&(ch1(i)=z) th=upper(ch1(i); ch3=ch3,th; else ch3=ch3,ch1(i); end end十一%将字符串中的小写字母变成相应的大写字母，其余字符不变。d=0; for j=1:b if (a=ch1(j)&(ch1(j)=z) d=d+1; end enddisp(小写字母个数为,num2str(d); end十二%从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。a=input(输入一个数);shuzu=;while(a=0) shuzu=shuzu,a; a=input(输入一个数);endle=length(shuzu);sum=0;for x=1:le sum=sum+shuzu(x);endsumavr=sum/leend十三%求100，200之间第一个能被21整除的整数。for k=100:1:200 a=rem(k,21); %取余是否为0 if a=0 break; endend disp(100到200之间第一个能被21整除的数是,num2str(k);end 十四%设 y=cos0.5+(3sinx)/(1+x2) 把x=02间分为101点，画出以x为横坐标，y为纵坐标的曲线；x=0:(2*pi)/100:2*pi;y=cos(0.5+(3*sin(x)./(1+x.2);plot(x,y)xlabel(x);ylabel(y);title(y=cos(0.5+(3*sin(x)./(1+x.2); 十五%在-10，10；-10，10范围内画出函数 的三维图形。x=-10:1/20:10;y=x; x,y=meshgrid(x,y); z=sin(x.2+y.2).0.5)./(x.2+y.2).0.5); surf(x,y,z) xlabel(x); ylabel(y); zlabel(z); title(z=sin(x.2+y.2).0.5)./(x.2+y.2).0.5)十六%利用Matlab求函数f(x)=-x2+ex+lnx的导数、不定积分和1到10区间的定积分及趋于1的极限。clc; clear all; close all;syms xfx = -x2 + exp(x) + log(x);dfx = diff(fx, x)ifx = int(fx, x)iabfx = int(fx, x, 1, 10)limfx = limit(fx, x, 1十七%用matlab函数解方程：x5+3=0；y5(t)+2*y=5t 。clear all; syms x y teq1 = x5 + 3;eq2 = y5 + 2*y - 5*t;solve(eq1, x)solve(eq2, y)十八%求方阵A所对应的行列式的值的函数是det(A)。b=1 2 3 4;2 3 4 1;3 4 1 2;4 1 2 3;det(b)det(b)c=b;c(:,1)=2*b(:,1);det(c)det(b(:,3 2 1 4)d=b;d(2,:);det(d)十九function y=fl(n)n=input(enter the order of a random matrix,please:);a=fix(100*rand(n)b=zeros(1,n*n),c=zeros(1,n*n),d=zeros(1,n*n);k=1;for i=1:n for j=1:n if rem(a(i,j),15)=0 & a(i,j)=0 b(k)=a(i,j);c(k)=i;d(k)=j; k=