2011年数学应用软件操作测验参考答案.docVIP

2011年数学应用软件操作测验 A 卷 编写一个函数，（1）计算，，（2）画出在上图像。 源码： %% 计算函数 function y=f(x) if x=0 y=0; elseif x2 y=0.5+1/x; else y=0.5; end %% 画图 x=-4:0.05:4; n=length(x); y=zeros(1,n); for i=1:n y(i)=f(x(i)); end plot(x,y,*) 结果： y=f(1) y = 0.5000 y2=f(f(4)) y2 = 0.5000 任给20个数，找出其中最小的数，并说明第几个数最小？要求两种方法（调用Matlab函数或写一段程序）。 法I：[xmin, nmin]=min(x) 法II：源码： function [xmin, nmin]=mymin(x) xmin=x(1); nmin=1; n=length(x); for i=2:n if xminx(i) xmin=x(i); nmin=i; end end 将图像窗口分成2幅，分别给出下列函数图形：（1），（2）。 源码： clc clear all clf subplot(1,2,1); x=-pi/2:0.05:pi/2; y=sqrt(cos(x)); plot(x,y) subplot(1,2,2); x=-2:0.1:2; y=-3:0.1:3; %x=linspace(-2,2,20); %y=linspace(-3,3,20); [X,Y]=meshgrid(x,y); Z=X.^2./4+Y.^2./9; %plot3(X,Y,Z); %surf(X,Y,Z); meshc(X,Y,Z); 已知在区间上11个采样点的函数值： 1 11 21 31 41 51 0 1.0414 1.3222 1.4914 1.6128 1.7076 61 71 81 91 101 1.7853 1.8513 1.9085 1.9590 2.0043 试求的5次拟合多项式，并绘出以为坐标的点和在区间上的函数曲线。要求加上标题：实验数据点及拟合曲线图例：*数据点 ---拟合曲线坐标轴标识：横坐标，纵坐标 法I、II，源码： clc clear all clf x=[1 11 21 31 41 51 61 71 81 91 101 ]; y=[0 1.0414 1.3222 1.4914 1.6128 1.7076... 1.7853 1.8513 1.9085 1.9590 2.0043 ]; Method I p=polyfit(x,y,5); % Method II: p(x)=p1*x^5+p2*x^4+...+p5*x+p6 %x=x; y=y; %A=[x.^5, x.^4,x.^3,x.^2,x,ones(length(x),1)]; %p=A\y; x1=1:101; y1=polyval(p,x1); plot(x,y,*,x1,y1,:); title(实验数据点(x_i,y_i) 及拟合曲线p(x) ); legend(数据点(x_i,y_i) ,拟合曲线p(x)) xlabel(x); ylabel(y); 法III：源码： clc clear all clf x=[1 11 21 31 41 51 61 71 81 91 101 ]; y=[0 1.0414 1.3222 1.4914 1.6128 1.7076... 1.7853 1.8513 1.9085 1.9590 2.0043 ]; plot(x,y,*); 然后，点击figure窗体Tools菜单中Basic Fitting，选择5th degree polynomial,OK. 用for循环计算 源码： clear all clc sum=0; fact=1; for n=1:10 fact=fact*n; sum=sum+fact; end disp([结果为,num2str(sum)]); 结果为4037913 列出100以内所有素数，即只能被1或本身整除的自然数。 法I： clear all clc n=input(Input a integer: ); disp(All prime numbers in [1,n] are ); for i=2:n yes=1; for j=2:fix(sqrt(i)) if rem(i,j)==0 yes=0; break;