华科电气matlab第一次作业

第一次作业（2015.10.27）目的：熟悉MATLAB的基本应用要求：用editor对程序进行编辑和调试，形成相应的脚本文件（*.m）。1.1若已知某电路参数，输入电压，绘制其输出电压的波形如图1，已知：R1=3Ω，R3=6Ω，C=1uF；us=18V，is=3A。在零初始条件下，开关s开路，当t=0时，开关S闭合，求uc，并画出波形；图1下列程序仅供参考（注：本例中的RC电路非图1所示的电路）：已知一阶RC电路的端电压的表达式为10s~0,62teut，试绘制电压波形u。t=0:0.1:10;u=6*exp(-2*t);plot(t,u);title('RCresponse');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('RCresponse');%以RCresponse为标题xlabel('t/s');%以t/s为横坐标ylabel('V/V')%以V/V为纵坐标grid2.已有两组测试数据不等的曲线，对其进行拟合和插值处理。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));%随机矩阵生成yp=polyfit(x,y,4)%拟合成4阶多项式y1=polyval(p,x);%求多项式y1的值plot(x,y,'+',x,y1,'-r')%绘制y与x的函数关系曲线title('fitting')%曲线命名为fitting2.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，是0到2pi，步长是0.1y1=interp1(xs,ys,xi,'linear');%使用linear法插值plot(xs,ys,'+k')%绘制ys与xs的曲线，黑色+号holdon%图形保持不变plot(xi,y1,':r')%同一坐标系中画出插值点的曲线，红色虚线legend('sampledpoint','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和脚本文件的应用。输入functiondy=vdp1000(t,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=1000*(1-y(1)^2)*y(2)-y(1);clearall[T,Y]=ode15s(@vdp1000,[03000],[20]);plot(T,Y(:,1),'-o')上交的作业文件包括：1.例程文件（需要有注释部分），2.对例程进行改动后的文件（注释说明那部分有修改即可）。