MATLAB与仿真上机练习题目

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

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

3、的解 ： ax2 bx c 0要求（1）a、b 、c 由键盘输入（ 2）输出的时候区分相同的根、实数根和复数根 ，即先输出“该二元一次方程具有相同的实数根/ 复数根 / 两个不同 / 实数根 / 复数根”，再输出具体的值 。. 学习帮手 .专业整理 .实验 4：自定义一个函数 mymax ，实现 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)c1c2e tc3 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，即第一次落地时所需时间t1t11=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=roots(0.

7、5.*g,v2,0);%当地面距离为0，即第二次落地时所需时间t2t2=-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-t11).2+v2*(t-N*t2-t11);.

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

9、lear 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: 求得的最大值% msg=n

10、argchk(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) maxvaluemaxvalue=x(ii);maxposition=ii;endend同时可判断一维、二维的情况function maxvalue,maxposition=mymax(x)% 构造自定义的max 函数 %x:输入向

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

12、 size_xndim_x=ndims(x);%求出 x 的维度 ：ndim_xif 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) maxvaluemaxvalue=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;endendendendendend实验 5：t=0,0.2,0.4,1.3,5.6;y=0.8,0.7,0.6,0.5,0.