MATLAB与仿真上机练习题目

1、.专业整理 .MATLAB与仿真练习实验 1：橡皮球从一个高度为 h0的高处，以初速度 v0 下落到地面，其高度 h、 速度 v 与时间 t 的关系由下式给出：12h(t) gt 2 v0t h02 0 0(1)v(t) gt v0(2)其中 g 为重力加速度( -9.81m/s 2)。用 matlab 程序编写， t 从 0 到500，以 1 为间隔，并绘出图像。 v0 、 h0应由键盘输入。. 学习帮手 .专业整理 .实验 2：一个 RLC振荡电路的负载公式是VR2 V0R2其中，2 f , f 为频率，单位是 hz。假设 L 0.1mH ， C 0.25nF ，R 50 ， V0 10m

2、V 。1. 绘制VR与 f 相关的图形，并根据图形给出 f 为何值时，VR取最大值，并计算最大值是多少？此时的 f 成为谐振频率 f 02. 若接收频率 f 大于 f0的 10%， VR为多少？3. 若VR为 VR最大值的一半，则此时接收频率 f 为多少？. 学习帮手 .专业整理 .实验 3：1. 请实现以下分段函数xyx0,y0x2 yx0,y0f(x,y) 2xyx0,y02x2 yx0,y0要求（ 1）.用 if 实现（2）. 用逻辑关系式实现绘制 f（x,y）与x,y的图像， x,y从-4 到 4, 步长为 0.1。2. 求以下一元二次方程的解： ax2 bx c 0要求（ 1）a、

3、b、c 由键盘输入（2）输出的时候区分相同的根、 实数根和复数根， 即先输出“该 二元一次方程具有相同的实数根 / 复数根 / 两个不同 / 实数根 / 复数 根”，再输出具体的值。. 学习帮手 .专业整理 .实验 4：自定义一个函数 myma，x 实现 matlab 函数 max 的功能，要求输出一 个参数时， 则输出最大值； 若要求输出两个参数， 输出最大值和最大 值的位置；若要求输出其他参数，则报错。1. 若输入向量为一个常数，则输出该数本身；2. 若输入向量为一维(行向量或者列向量) ，则输出最大值；3. 若输入向量为二维，则输出每一列的最大值4. 若输入向量为三维以上，则输出每一按照

4、列取出的维度的最大值(课后思考)实验 5：对 一 组 测 量 数 据 t=0,0.2,0.4,1.3,5.6 ， y=0.8,0.7,0.6,0.5,0.2， 用 以 下 函 数 来 拟 合 这 组 数 据 ：y(t) c1 c2e c3 cost ，画出拟合结果，并给出 c1 c2 c3 的值。实验 6：对正弦波 sin(2 * * t) ，完成如下任务： 1对该正弦波一个周期采样 1024 点，画出波形； 2对每个采样点量化成 10 位，并按照二进制补码进行编码；. 学习帮手 .专业整理 .3将编码后的 1024 个样点保存在 sample1024.txt 文件中； 提示： 正弦波的处理过

5、程如下：正弦波采样 量化 编码采样：将正弦波一个周期 0,2*pi ，采样 1024 点，采样后的数据是 浮点数，数据范围是 -1 ，1 ； 量化：将采样后的数据从浮点型转换成定点型， 即转换成 10bit 整型 数，数据范围将变成 -512,511 ；编码：将量化后的 10bit 整型数按照二进制补码规则进行编码；. 学习帮手 .专业整理 . 学习帮手 .专业整理 .实验结果实验 1：workxiangpiqiu.mclcclear all;%v0=input( 请输入初始速度 v0);%h0=input( 请输入初始高度 h0); %以正上方为正速度，以下方为负速度 v0=10;h0=50

6、;g=-9.81;t=0:0.01:50; % 以 0.01 为间隔，共 500s% h=0.5.*g.*t.2+v0*t+h0;%离地面距离 ht1=roots(0.5.*g,v0,h0);%当地面距离为 0 ，即第一次落地时所需时间 t1即第t11=min(t1(t10);% 对求解的根 t 做如下处理： 小于等于 0的t 舍弃， 大于 0的t 取最小值， 次落地时真正的时间 t11v1=g.*t11+v0;% 即第一次落地时的速度 v1%假%设%无%能量损耗，完全刚性，即以同样速度v1 反弹回空中v2=-v1;% h2=0.5.*g.*t.2+v2*t;% 离地面距离 h% t2=roo

7、ts(0.5.*g,v2,0);% 当地面距离为 0，即第二次落地时所需时间 t2 t2=-v2./(0.5.*g);% 当地面距离为 0，即第二次落地时所需时间 t2 。与上式等价%求%解%h和 v% %实际为一分段函数，当 tt11 时候， v=g.*t+v0;h=0.5.*g.*t.2+v0*t+h0;%当 t11tt11+t2 时候，v=g.*(t-t11)+v2;h=0.5.*g.*(t-t11).2+v2*(t-t11);. 学习帮手 .专业整理 .% 当 t11+N*t2tt11+(N+1)*t2 时候，v=g.*(t-N*t2-t11)+v2;h=0.5.*g.*(t-N*t2

8、-t11).2+v2*(t-N*t2-t11);将当 tt11 时候与其他时候的 t 分开% x=(tt11).*mod(t-t11,t2);x1=(tt11).*mod(t-t11,t2);v22=g.*x2+v2;h22=0.5.*g.*x2.2+v2*x2;%第%一%次%落地之后的 h和 v v=v11.*(tt11); h=h11.*(tt11);subplot(2,1,1)plot(t,v)xlabel( 时间 );ylabel( 皮球速度，上为正，下为负 ); subplot(2,1,2)plot(t,h)xlabel( 时间 );ylabel( 皮球离地面的高度 );. 学习帮手

9、 .专业整理 .实验 3：1(2)clear all;x,y=meshgrid(-4:0.1:4,-4:0.1:4);f1=(x+y).*(x=0&y=0)+(x.2+y.2).*(x0&y=0&y0)+(x.2+y).*(x =0);% mesh(x,y,f1);surf(x,y,f1);% f2=x.*(x=0)+(x.2).*(x=0)+(y.2).*(y0);% 两者效果相同实验 4：workmymax.m一维情况function maxvalue,maxposition=mymax(x) %构%造%自%定义的 max函数 % %x:输入向量，可以为任意维度%maxvalue: 求得的

10、最大值% msg=nargchk(1,2,nargout);% err(msg);. 学习帮手 .专业整理 .% %maxposition: 最大值所在的位置% if nargout2% error( 输出参数应该是 1 或者 2);% end% size_x=size(x);maxvalue=x(1); maxposition=1;for ii=1:size(x,2)if x(ii) maxvalue maxvalue=x(ii); maxposition=ii;endend同时可判断一维、二维的情况function maxvalue,maxposition=mymax(x) %构%造%自%

11、定义的 max函数 % %x:输入向量，可以为任意维度%maxvalue: 求得的最大值 %maxposition: 最大值所在位置% %maxposition: 最大值所在的位置% 如果 x 是一个向量，那么 maxvalue 便是这个向量元素中的最大值；%1*n 的矩% 如果 x 是一个二维矩阵，那么 maxvalue 便是选出每一列中的最大值，返回一个 阵，或者说是一个 n 维行向量。% msg=nargchk(1,2,nargout);% err(msg);% if nargout2% error( 输出参数应该是 1 或者 2);% end size_x=size(x);% 求出 x

12、 的 size ： size_x ndim_x=ndims(x);% 求出 x 的维度： ndim_x if ndim_x=1%为%一%个%常x数的情况 %. 学习帮手 .专业整理 .maxvalue=x; maxposition=1;elseif ndim_x=2%为%一%维%向%量x 的情况 % if size_x(1)=1|size_x(2)=1maxvalue=x(1); maxposition=1;for ii=1:length(x) if x(ii) maxvalue maxvalue=x(ii); maxposition=ii;endendelse%为%x二维矩阵的情况 %maxvalue=zeros(1,size_x(2);maxposition=zeros(1,size_x(2);for jj=1:size_x(2)y=x(:,jj);maxvalue(jj)=y(1);maxposition(jj)=1;for ii=1:size_x(1)if y(ii) maxvalue(jj)maxvalue(jj)=y(ii);maxposition(jj)=ii;endend end end end end实验 5：t=0,0.2,0.4,1.3,5.6;y=0.8,0.7,0.6,0.5,