第三章 MATLAB期末作业.docx

MATLAB期末作业1．从键盘输入一个4位整数，按如下规则加密后输出。加密规则：每位数字都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 a=input(请输入一个四位数a=?);while (a1000|a9999); disp(输入错误，请重新输入！); a=input(请输入一个四位数a=?); end%要求输入的是4位整数，输入错误时要求重新输入 a1=fix(a/1000);%求a的千位数字 a2=rem(fix(a/100),10);%求a的百位数字 a3=rem(fix(a/10),10); %求a的十位数字 a4=rem(a,10);%求a的各位数字 b1=rem(a3+7,10);%以下用加密规则进行加密 b2=rem(a4+7,10); b3=rem(a1+7,10); b4=rem(a2+7,10);disp([b1,b2,b3,b4]);运行程序为请输入一个四位数a=?1326（此为随机输入的数字） 9 3 8 0实验四第4题.已知求中：最大值、最小值、各数之和。正数、零、负数的个数。解：f(1)=1;f(2)=0;f(3)=1;for n=4:100f(n)=f(n-1)-2.*f(n-2)+f(n-3);enddisp([max=,num2str(max(f))])disp([min=,num2str(min(f))])disp([sum=,num2str(sum(f))])k1=0;k2=0;k3=0;for n=1:100 if f(n)0 k1=k1+1; end if f(n)==0 k2=k2+1; end if f(n)0 k3=k3+1; end enddisp([正数个数为,num2str(k1)])disp([零的个数,num2str(k2)])disp([负数个数为,num2str(k3)])实验六第4题：绘制极坐标曲线，并分析参数对曲线形状的影响。解：a=input(a=?);b=input(b=?);n=input(n=?);theta=0:0.01:2*pi;rho=a.*sin(b+n*theta);polar(theta,rho,k)结果为：当输入a=1,b=1,c=1时，图像为当输入a=1,b=1,c=2时，图像为当输入a=1,b=1,c=3时，图像为当输入a=1,b=1,c=4时，图像为当输入a=1,b=2,c=1时，图像为当输入a=1,b=3,c=1时，图像为当输入a=1,b=4,c=1时，图像为当输入a=2,b=1,c=1时，图像为当输入a=3,b=1,c=1时，图像为综上：采用控制变量法绘制图像，参数a是控制图像的极经；参数b对图像影响不大，只是对图像取值位置有所影响；参数c是控制函数周期，数值越大，画出的图像越密集。1.利用MATLAB提供的rand函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：（1）均值和标准方差。（2）最大元素和最小元素。（3）大于0.5的随机各数占总数的百分比。解：A=rand(300,100);a=max(A);a1=max(a);b=min(A);b1=min(b);c=mean(A);c1=mean(c);d=var(A(:));t=sum(sum(A0.5));disp([大于0.5的个数为,num2str(t)])disp([最大元素为,num2str(a1)])disp([最小元素为,num2str(b1)])disp([平均值为,num2str(c1)])disp([标准方差为,num2str(d)])结果为：大于0.5的个数为15015最大元素为0.99997最小元素为4.8345e-005平均值为0.49866标准方差为0.0834773.某气象观测站测得某日6:00~18:00之间每隔2h的室内外温度(℃)如实验表1所示时间h681012141618室内温度t118.020.022.025.030.028.024.0室外温度t215.019.024.028.034.032.030.0试用三次样条插值分别求出该日室内外6:00~17:30之间每间隔2h各点近似温度(℃)。解：H=6:2:18;X=6.5:2:17.5;F1=[18.0,20.0,22.0,25.0,30.0,28.0,24.0];F2=[15.0,19.0,24.0,28.0,34.0,32.0,30.0];F3=interp1(H,F1,X,spline)F4=interp1(H,F2,X,spline)结果为