华科电气matlab第一次作业.docVIP

第一次作业（2015.10.27） 目的：熟悉MATLAB的基本应用 要求：用editor对程序进行编辑和调试，形成相应的脚本文件（*.m）。 1.1若已知某电路参数，输入电压，绘制其输出电压的波形 如图1，已知：R1=3Ω，R3=6Ω，C=1F；us=18V，is=3A。 在零初始条件下，开关s开路，当t=0时，开关S闭合，求uc，并画出波形； 图1 下列程序仅供参考（注：本例中的RC电路非图1所示的电路）： 已知一阶RC电路的端电压的表达式为 ，试绘制电压波形u。 t=0:0.1:10; u=6*exp(-2*t); plot(t,u); title(RC response); xlabel(t/s); ylabel(V/V) grid 解得u = -12+12/exp(t/(2*(e - 6))) t=0:1e-8:1e-5; %以10-8为间隔,在0到10-5区间内取点 u=12*exp(-2e6*t)-12; %算出t时刻对应的电压值 plot(t,u); %以t为横坐标,u为纵坐标画图 title(RC response); %以RC response为标题 xlabel(t/s); %以t/s为横坐标 ylabel(V/V) %以V/V为纵坐标 grid 2. 已有两组测试数据不等的曲线，对其进行拟合和插值处理。 2.1给例程2.1加入注释，说明程序的用途和调试结果，尝试修改x，y，或polyfit的不同设置。熟悉fitting工具 例程2.1： clear;clc;close; %删除工作空间的所有变量，删除命令窗口的所有变量，关掉图形 x=0:0.1:2*pi; %在0到2pi范围内，以0.1为步长，随机生成样本点 y = sin(x)+0.5*rand(size(x)); %随机矩阵生成y p = polyfit(x,y,4) %拟合成4阶多项式 y1= polyval(p,x); %求多项式y1的值 plot(x,y, +,x,y1, -r) %绘制y与x的函数关系曲线 title(fitting) %曲线命名为fitting 2.2 给例程2.2加入注释，说明程序的用途和调试结果，尝试修改xs，ys，或xi，或interp1的不同设置，熟悉插值的用法。 例程2.2： clear;clc;close; %删除工作空间的所有变量，删除命令窗口的所有变量，关掉图形 xs=0:0.2:2*pi; %在0到2pi之间，以0.2为步长，随机生成样本点 ys= sin(xs); % 生成样本点ys=sin（xs） xi=0: 0.1: 2*pi; %生成插值点xi， y1=interp1(xs,ys,xi,linear); %使用linear法插值 plot(xs,ys,+k) %绘制ys与xs的曲线，黑色+号 hold on %图形保持不变 plot(xi,y1,:r) % 同一坐标系中画出插值点的曲线， legend(sampled point,linear,spline); %添加图例 title(INTERP) %添加题目 修改interp1: 即y1=interp1(xs,ys,xi,nearest) 3. 试用solve求解常系数微分方程。 3.1 熟悉dsolve的用法，并尝试例程3中的微分方程（组），再求解。 例程3： dsolve(Dy=1+y^2)y=dsolve(Dy=1+y^2,y(0)=1,x)x=dsolve(D2x+2*D1x+2*x=exp(t),x(0)=1,Dx(0)=0)S=dsolve(Df=3*f+4*g,Dg=-4*f+3*g) 在matlab中输入如上代码得出结果： 3.2 对比加深学习，阅读下面例程example2，学习数值微分方程的求解方法，体会与dsolve的不同，学会ode45应用，熟悉函数function和脚本文件的应用。 输入 function dy=vd