MATLAB习题.docx

第一次上机作业：计算y=2sin0,3e31+5lg6clc;cleary=(2*sin(0.3*pi)*exp(3)/(1+sqrt(5)*log10(6)1. 生成等差矩阵a=1,3,5,747,49;在经度122,128间生成70个等差数列。a=(1:2:50) %生成等差矩阵a=1,3,5,747,49b=linspace(122,128,70) %在经度122,128间生成70个等差数列2. 生成34均匀随机矩阵b；提取b的第3行第2列元素赋值为5，求b中最小值、平均值；提取b的第三列；b右边增添1列为 2 3 4,；删除b的第二列。b=rand(3,4) %生成34均匀随机矩阵bb(3,2)=5 %提取b的第3行第2列元素赋值为5m=min(b(:) %求b中的最小值n=mean(b(:) %求b中的平均值a=b(:,3) %提取b的第三列b(:,5)=2 3 4 %在b右边增添1列为 2 3 4b(:,2)= %删除b的第二列3. 在区间x=0,2内做sin(x)图；按x、sin(x)、cos(x)存成3列到data.txt文档；再读取文件data.txt，用读取的数据做cos(x)图。x=0:0.1*pi:2*pim=sin(x)n=cos(x)plot(x,m,r-*) %做sin(x)图像a=x;sin(x);cos(x)save -ascii data.txt a %按x、sin(x)、cos(x)存成3列到data.txt文档q=load(data.txt) %读取文件data.txtplot(q(:,1),q(:,3) %用读取的data.txt做cos(x)图第二次上机作业：1. 生成34，-1,1间的随机矩阵A；求A矩阵中的元素总个数，将A中小于0的元素等于1；再取出A中大于0.2并小于0.7的元素。clc;clearA= -1+(1+1)*rand(3,4) %生成34，-1,1间的随机矩阵Alength(A(:) %求A矩阵中的元素总个数A(A0.2&A0是x-y折线图。clc;clearx=0:0.1:2*pi % X在0:0.1: 2内y=cos(x.2) %求y=cos(x2 )值a=find(y0) % y0m=x(a);n=y(a);plot(m,n,r-*) %作x-y折线图3. 首先由计算机随机产生一个1,100之间的一个整数，然后由用户猜测所产生的这个数。根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won！”，同时退出这个游戏，如果用户7次还没有猜出，则显示“You lost！”，并显示出那个随机数。（用while编写）clc;clearx=round(1+(100-1)*rand(1,1); %随机产生一个1,100之间的一个整数n=1;while nx disp (High) end if y2)&(data(:,2)=124)&(lat=28);lat(b);lon(b);plot(lat(b),lon(b),.)hold onboundaxis(120 130 27.5 32.5)title(28N以北鲐鱼幼鱼地理分布图,fontsize,20)2. 文件cjk.nc为模型输出的长江口netcdf数据，绘制24个小时node=5019所有层平均的水温波动折线图；在地图上绘制cell=10661(第2层)和cell=1231(第7层)的24个小时的流场变动图；绘制time=9时刻整场的温度等值线(表层)、盐度假色彩图(所有层平均)、流场图(第3层)（标注尺度大小）。clc;clear;clfaddpath(genpath(netcdf_tools)nc=netcdf(cjk.nc,nowrite);wd=nctemp(:,:,5019);pwd=mean(wd);sj=1:24;plot(sj,pwd)title(水温波动折线图)xlabel(时间)ylabel(平均水温)lon=nclonc(10661);lon1=nclonc(1231);lat=nclatc(10661);lat1=nclatc(1231);uu=ncu(:,2,10661);uu1=ncu(:,7,1231);vv=ncv(:,2,10661);vv1=ncv(:,7,1231);lon=repmat(lon,1,24)lon1=repmat(lon1,1,24)lat=repmat(lat,1,24)lat1=repmat(lat1,1,24)myquiver=quiver(lon,lat,uu,vv)hold onbound;myquiver=quiver(lon1,lat1,uu1,vv1)axis(120 125 30 35)title(流场变动图)xlabel(经度)ylabel(纬度)t=nctemp(9,1,:)lon2=nclon(:)lat2=nclat(:)xlon=min(lon2):0.1:max(lon2)xlat=min(lat2):0.1:max(lat2)xt=griddata(lon2,lat2,t,xlon,xlat,cubic)c h=contour(xlon,xlat,xt)clabel(c,h)hold onbound;axis(120 125 26 35)title(温度等值线)xlabel(经度)ylabel(纬度)yd=ncsalinity(9,:,:)lon3=nclon(:)lat3=nclat(:)xlon=min(lon3):0.1:max(lon3)xlat=min(lat3):0.1:max(lat3)yd=mean(yd)yd=griddata(lon3,lat3,yd,xlon,xlat,cubic)pcolor(xlon,xlat,yd)shading interphold onboundaxis(120 125 26 35)title(盐度假色彩图)xlabel(经度)ylabel(纬度)v=ncv(9,3,:)u=ncu(9,3,:)lonc=nclonc(:)latc=nclatc(:)m_sca=0.2quiver(lonc,latc,u*m_sca,v*m_sca,autoscale,off)hold onboundquiver(120,29,1*m_sca,0,autoscale,off)text(120,3,29,1m/s)axis(120 125 26 35)title(流场图)xlabel(经度)ylabel(纬度) close(nc);第五次上机作业：1. 将文件cjk.nc中24个时刻、经度【121 123】、纬度【30 33】、第3层流场数据用fprintf写入文件，文件名为“vector001.dat”、“vector002.dat”、.、“vector024.dat”；将最后时刻上5层平均温度和下3层平均盐度数据写入Execel文件ts.xls里。（每个文件都包括经度、纬度以及数据）。clc;clear;addpath(genpath(netcdf_tools) ;nc=netcdf(cjk.nc,nowrite);lon=nclonc(:);lat=nclatc(:);m=find(lon=121&lon=30&lat=33);lon=lon(m);lat=lat(m);layer=3;for time=1:24 u=ncu(time,layer,:); v=ncv(time,layer,:); u=u(m); v=v(m); dat=lon lat u v; fid=fopen(tmp/vector,num2str(time,%3,3d),.dat,wt); fprintf(fid,%8.3f %8.3f %8.3f %8.3fn,dat); fclose(fid);endclose(nc);clc;clear;addpath(genpath(netcdf_tools) ;nc=netcdf(cjk.nc,nowrite);time=24;lon=nclon(:);lat=nclat(:);layer=1:5;temp=nctemp(time,layer,:);temp=mean(temp);layer=8:10;sal=ncsalinity(time,layer,:);sal=mean(sal);a=lon lat temp sal;xlswrite(ts.xls,a)close(nc);2. 将鲐鱼仔鱼输运轨迹文件larvae_tacking.dat，用textscan和fprintf改写成能够用load直接读入的单个文本文件，文件名格式为“larvae_0001.txt”、“larvae_0002.txt”、 “larvae_0018.txt”、；再在地图上作出仔鱼输运的动图clc;clear;fid1=fopen(larvae_tracking.dat,rt);dat0=textscan(fid1,%s %s %s %s %s %s,1);for i=2:1000 if feof(fid1) break; end dat1=textscan(fid1,%s %d %d %d %d %s %d,1); n=dat17; data=textscan(fid1,%f %f,n); dat1=data1 data2; fid2=fopen(tmp2/larvae_,num2str(i,%3.3d),wt); fprintf(fid2,%10.5f %10.5f n,dat1);