matlab编程考试题及答案

matlab编程考试题及答案一、选择题（每题2分，共20分）1.以下哪个变量名符合MATLAB命名规则？A.2nd_varB.@funcC.my-varD.result_32.执行`A=[12;34];B=A(2,:);`后，B的维度是？A.1×2B.2×1C.2×2D.1×13.对于逻辑数组`L=[truefalsetrue]`，`find(L)`的输出是？A.[13]B.[02]C.[24]D.[123]4.要绘制散点图，应使用以下哪个函数？A.plotB.scatterC.barD.histogram5.执行`x=linspace(0,pi,5);`后，x的第3个元素是？A.π/4B.π/2C.3π/4D.π6.关于MATLAB函数句柄，以下说法错误的是？A.用`@`符号创建B.可作为函数参数传递C.只能指向内置函数D.可用于匿名函数7.读取CSV文件`data.csv`的正确函数是？A.loadB.csvreadC.fopenD.readmatrix8.向量化操作相比循环的主要优势是？A.代码更易读B.计算速度更快C.支持更复杂的逻辑D.内存占用更少9.结构体数组`S(1).name='Alice';S(2).name='Bob';`中，`[S.name]`的结果是？A.'AliceBob'B.{'Alice','Bob'}C.['Alice','Bob']D.错误10.捕获异常时，应使用以下哪组关键字？A.if-elseB.for-endC.try-catchD.switch-case二、填空题（每题3分，共30分）1.提供一个3×3的全1矩阵的命令是__________。2.计算矩阵`A=[12;34]`的转置矩阵的命令是__________。3.执行`x=[1357];y=x(2:end-1);`后，y的值为__________。4.编写匿名函数计算`f(x,y)=x²+y³`，应使用__________。5.绘制正弦曲线`y=sin(x)`（x从0到2π，步长0.1）的代码段为__________。6.读取文本文件`log.txt`并存储为字符串变量的命令是__________。7.计算积分`∫₀¹x²dx`的数值解，应使用__________函数。8.求解微分方程`dy/dt=-2y`（初始条件y(0)=1）的数值解，应调用__________函数。9.单元数组`C={1,'apple',[246]}`中，提取字符串'apple'的索引是__________。10.对数组`v=[5281]`进行升序排序后，排序结果为__________，原索引为__________。三、编程题（共50分）1.（10分）编写函数`stats_calculate`，输入一个一维数值数组，输出其均值、方差、最大值和最小值。要求：输入参数为数组`x`（长度≥2）；输出参数为结构体`result`，包含字段`mean`、`var`、`max`、`min`；不使用内置函数`mean`、`var`、`max`、`min`，需自行计算。2.（12分）某传感器采集了一组温度数据`T=[23.1,24.5,22.8,25.3,21.7,26.0,23.9,24.1,22.5,25.8]`（单位：℃）。要求：绘制温度随时间变化的折线图，x轴标注“时间（小时）”，y轴标注“温度（℃）”，标题为“温度变化趋势”；计算温度超过25℃的时刻（索引），并在图中用红色五角星标记这些点；计算所有时刻的平均温度（保留2位小数），并在图中添加文本标注显示平均值。3.（12分）已知二维矩阵`A=[413;257;869]`，要求：找出所有大于5的元素的位置（行、列索引）；将这些位置的元素替换为对应行的均值；输出修改后的矩阵`A_new`。4.（16分）某信号处理场景中，采集到离散信号`y(n)=2sin(0.2πn)+0.5cos(0.8πn)+w(n)`，其中`n=0:100`，`w(n)`为高斯白噪声（均值0，方差0.1）。要求：提供信号`y`，并绘制其时域波形图；对信号进行快速傅里叶变换（FFT），绘制幅频特性图（频率范围0到0.5倍采样频率）；设计一个低通滤波器（截止频率0.3π），对信号进行滤波，绘制滤波前后的时域波形对比图（叠加显示）。答案一、选择题1.D（变量名需以字母开头，可包含字母、数字、下划线）2.A（`A(2,:)`提取第2行所有列，结果为1×2行向量）3.A（`find`返回逻辑为真的索引，从1开始计数）4.B（`scatter`用于散点图，`plot`为折线图）5.B（`linspace(0,pi,5)`提供5个点，间隔为π/4，第3个元素为0+2×(π/4)=π/2）6.C（函数句柄可指向自定义函数或匿名函数）7.D（`readmatrix`是新版本推荐的读取表格数据的函数，`csvread`已逐渐弃用）8.B（向量化利用MATLAB底层优化，通常比循环快）9.A（`[S.name]`将结构体数组的`name`字段按顺序连接为字符串）10.C（`try-catch`用于异常处理）二、填空题1.`ones(3)`或`ones(3,3)`2.`A'`或`transpose(A)`3.`[35]`（`2:end-1`表示第2个到倒数第2个元素）4.`f=@(x,y)x.^2+y.^3`5.`x=0:0.1:2pi;y=sin(x);plot(x,y);`6.`str=fileread('log.txt');`7.`integral(@(x)x.^2,0,1)`8.`ode45(@(t,y)-2y,[010],1)`（时间范围可根据需求调整）9.`C{2}`10.`[1258]`；`[4213]`（`[sorted,idx]=sort(v)`）三、编程题1.函数实现：```matlabfunctionresult=stats_calculate(x)n=length(x);%计算均值sum_x=0;fori=1:nsum_x=sum_x+x(i);endmean_val=sum_x/n;%计算方差（无偏估计）sum_var=0;fori=1:nsum_var=sum_var+(x(i)mean_val)^2;endvar_val=sum_var/(n1);%计算最大值和最小值max_val=x(1);min_val=x(1);fori=2:nifx(i)>max_valmax_val=x(i);endifx(i)<min_valmin_val=x(i);endend%构建结构体result.mean=mean_val;result.var=var_val;result.max=max_val;result.min=min_val;end```2.温度数据分析代码：```matlabT=[23.1,24.5,22.8,25.3,21.7,26.0,23.9,24.1,22.5,25.8];t=1:length(T);%时间索引（小时）%绘制折线图figure;plot(t,T,'b-o','LineWidth',1.5);xlabel('时间（小时）');ylabel('温度（℃）');title('温度变化趋势');gridon;%标记超过25℃的点idx=find(T>25);holdon;scatter(t(idx),T(idx),80,'r','filled','Marker','');holdoff;%计算平均温度mean_T=round(mean(T),2);%添加文本标注（位置可调整）text(0.8max(t),max(T)0.95,sprintf('平均温度：%.2f℃',mean_T),'FontSize',10);```3.矩阵操作代码：```matlabA=[413;257;869];%找出大于5的元素位置[row,col]=find(A>5);%计算每行均值row_means=mean(A,2);%3×1向量%替换对应位置元素fori=1:length(row)r=row(i);c=col(i);A(r,c)=row_means(r);endA_new=A;%输出结果（示例输出）%A_new=%4.00001.00003.0000%2.00005.00004.3333（原第2行均值为(2+5+7)/3=14/3≈4.3333）%7.66677.66679.0000（原第3行均值为(8+6+9)/3=23/3≈7.6667）```4.信号处理代码：```matlab%提供信号n=0:100;L=length(n);w=0.1randn(1,L);%高斯白噪声y=2sin(0.2pin)+0.5cos(0.8pin)+w;%时域波形图figure;subplot(3,1,1);plot(n,y);xlabel('n');ylabel('y(n)');title('原始信号时域波形');%FFT分析（假设采样频率Fs=1）Fs=1;Y=fft(y);P2=abs(Y/L);P1=P2(1:L/2+1);P1(2:end-1)=2P1(2:end-1);f=Fs(0:(L/2))/L;subplot(3,1,2);plot(f,P1);xlabel('频率（Hz）');ylabel('幅度');title('幅频特性');xlim([00.5]);%低通滤波器设计（截止频率0.3π，对应归一化频率0.3）fc=0.