如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
83ye4、编程实现当t∈[-π,π],间隔为1°时求解正弦和余弦值。解:t=(-1*pi:pi/180:pi);y1=sin(t)y2=cos(t)5、利用rand函数产生(0,1)间的均匀分布的10*10随机矩阵A,然后统计A中大于等于0.6的元素的个数。解:A=rand(10);B=A>=0.6;C=sum(B);count=sum(C)运行结果(每次运行结果是不同的,仅作参考):count=326、利用randn函数产生均值为0,方差为1的10*10随机矩阵A,然后统计A中大于-0.5且小于0.5的元素的个数。解:A=randn(10);B=(A<0.5)&(A>-0.5);C=sum(sum(B))运行结果(每次运行结果是不同的,仅作参考):C=488、有一矩阵A,找出矩阵中值大于1的元素,并将它们重新排列成列向量B。解:A=2*rand(4);k=find(A<=1);A(k)=[];%删除下标为k的元素B=A'9、在一测量矩阵A(100*3)中,存在有奇异值(假设大于100的值认为是奇异值),编程实现删去奇异值所在的行。解:A=120*randn(100,3);[i,j]=find(A>100);A(i,:)=[]%删去存在奇异值的行运行结果(每次运行结果是不同的,仅作参考):A=49.5355-23.7550-73.0269-118.435439.3214-88.447291.1482-28.5962-209.985521.2336-74.0239-9.5871-15.818472.132222.044471.442911.077034.89484、在给定的100*100矩阵中,删去整行为0的行,删去整列为0的列。解:A=diag([1234],1)B=any(A)[i,j]=find(B==0)A(:,i)=[]%删除全为0的列B=any(A')[i,j]=find(B==0)A(j,:)=[]%删除全为0的行运行结果:初始值:A=0100000200000300000400000操作后:A=1000020000300004146ye2、将图形窗口分成两个,分别绘制正割和余割曲线,并加上标注。程序为:x=0:pi/10:2*pi;figure(1)subplot(2,1,1)plot(x,sec(x),'k-'),gridonlegend('\ity=sec(x)')title('y=sec(x)')xlabel('x'),ylabel('y')subplot(2,1,2)plot(x,csc(x),'k-'),gridonlegend('\ity=csc(x)')title('y=csc(x)')xlabel('x'),ylabel('y')3、设有函数y=exp(x+5)+x.^3,在半对数坐标系中绘制曲线。程序为:x=1:0.01:10;y=exp(x+5)+x.^3;figure(1)subplot(2,1,1)plot(x,y,'r-'),gridonlegend('\ity=exp(x+5)+x.^3')title('平面坐标')xlabel('x'),ylabel('y')subplot(2,1,2)semilogx(x,y,'k-'),gridon%半对数坐标轴legend('\ity=exp(x+5)+x.^3')title('半对数坐标')xlabel('x'),ylabel('y')5、将图形窗口分成两个,分别绘制以下函数在[-3,3]区间上的曲线,并利用axis调整轴刻度,使他们具有相同缩放尺度。y1=2x+5;y2=x2-3x+1。程序为:x=-3:0.1:3;y1=2*x+5;y2=x.^2-3*x+1;figure(1)subplot(2,2,1)plot(x,y1,'r-'),gridonlegend('\ity1=2*x+5')title('y1=2x+5')xlabel('x'),ylabel('y1')subplot(2,2,2)plot(x,y2,'k-'),gridonlegend('\ity2=x.^2-3*x+1')title('y2=x^2-3x+1')xlabel('x'),ylabel('y2')subplot(2,2,3)plot(x,y1,'r-'),gridonlegend('\ity1=2*x+5')title('调整后的y1=2x+5')axis([-33-1010])xlabel('