MATLAB-09电子、机器人(A卷).docVIP

一、（20分） 在正态随机数发生器的状态为1的情况下，产生（10*10）的正态分布随机矩阵。给出该矩阵中在(0.5,0.6)范围内的元素所在的位置，分别写出“全下标”及“单下标”，并统计个数。（保存程序m文件，结果显示在指令窗口） clear; randn(state,1); %初始状态设定2分 A=randn(10,10); %产生矩阵3分 L=(A0.5)(A0.6);%关系运算3分 si=find(L) %单下标3分 [ri,ci]=ind2sub(size(A),si) %全下标3分 count=length(si) %个数统计2分 结果： si =21 29 r =1 9 c =3 3 文件名命名1分，数据结果显示在指令窗口 3分 二、（15分） 绘制函数 的曲线，求在[0,10]范围内的所有实数解，要求函数表达式用内联函数和匿名函数两种形式分别求解。（保存程序m文件，结果显示在指令窗口） clear; t=0:0.01:10; plot(t,-sin(t)+cos(t).^2-0.5*t+1); %绘图观察零点 2分 f=inline(-sin(t)+cos(t).^2-0.5*t+1); %内联对象定义3分 result1=fsolve(f,[1 3 4]) %求解4分 f=@(t)(-sin(t)+cos(t).^2-0.5*t+1); %匿名函数定义3分 result2=fsolve(f,[1 3 4]) 文件名命名1分，数据结果显示在指令窗口 2分 结果： result1 = 0.9749 2.6664 4.2138 result2 = 0.9749 2.6664 4.2138 三、（20分） 现有一个函数，求该函数的积分式，要求(1)采用数值计算方法求积分，并用双坐标轴绘出被积分函数曲线及积分结果的曲线，时间范围为[0,15]。给出横纵坐标名称标注，其中纵坐标名称与函数表达式一致。(2)采用数值计算方法求，要求保证结果的绝对精度为10-9（提示利用quad函数）。（保存程序m文件，数据结果显示在指令窗口） clear; x1=0:0.1:15; f=exp(-x1); y1=cumtrapz(f)*0.1;%累积积分式3分 h=plotyy(x1,f,x1,y1);%绘图指令 3分 set(get(h(1),Ylabel),string,f(x)=e^{-x}); %坐标轴设置3分 set(get(h(2),Ylabel),string,y(x)=\int_0^x f(t)dt) %坐标轴设置3分 xlabel(x); %坐标轴设置1分 format long %显示格式设置2分 y3=quad(@(x)exp(-x),0,0.5,1e-9) %积分指令3分 y3 = 0.39346934028827 文件名命名1分，数据结果显示在指令窗口 1分 四、（20分） 利用数值计算方法求解下面的微分方程： ，初始条件为 设定时间范围为[0,5]，要求：(1) 产生两个子图，第一个子图是函数的曲线，添加标题“微分方程解”；第二个子图是关系的相平面图，添加标题“相平面图”。要求两子图中给出横纵坐标名称。 (2)求的值，要求显示小数点后15位。（保存程序m文件，数据结果显示在指令窗口） clear; t=[0,5]; [t,y]=ode45(@fun,t,[2;0]);%指令3分 subplot(1,2,1);plot(t,y(:,1));%子图1分，绘图2分 xlabel(t);ylabel(x);%坐标名称1分 title(微分方程解);%标题1分 subplot(1,2,2);plot(y(:,1),y(:,2)); %绘图2分 xlabel(x);ylabel(dx/dt); %坐标名称1分 title(相平面图); %标题1分 [t1,y1]=ode45(@fun,[0,1.5],[2;0]); y1(end,1) %结果2分 fun.m函数文件 %函数文件格式2分 function dxdt= fun(t,x) dxdt=[x(2);x(2)-t*x(1)+1]; %微分方程表达2分 end ans = 2.135976645059067 文件名命名1分，数据结果显示在指令窗口 1分 五、（15分） 在matlab_final文件夹下有一个data.mat的数据文件，要求使用指令导入该文件中数据，并根据这组数据拟合5阶多项式。绘出数据点及拟合多项式曲线，添加图例对应拟合曲线或数据点，添加横坐标名称。（保存程序m文件） clear; data=importdata(data.mat); %导入数据指令2分