资源描述:
《Matlab简单实例学习.doc》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、Matlab程序代码绘制的函数图象。fvclear;t=0:0.02:10;f1=10/sqrt(7.75).*exp(-1.5*t);f2=sin(sqrt(7.75).*t);y=f1.*f2;plot(t,y,'-k',t,y,'ok');xlabel('t');ylabel('y(t)');title('函数图像')axis([-210-0.52])拉氏变换clear;clc;symsstfs1fs2fs3ft1ft2ft3;L=1,C=0.1,R=[1.535];h1=1/(L*C*s^2+R(1)*C*s+1);h2=1/(L*C*s^2+R(2)*C*s+1);h3
2、=1/(L*C*s^2+R(3)*C*s+1);fs1=h1*(1/s);fs2=h2*(1/s);fs3=h3*(1/s);ft1=ilaplace(fs1,s,t);ft2=ilaplace(fs2,s,t);ft3=ilaplace(fs3,s,t);ezplot(t,ft1);holdon;ezplot(t,ft2);holdon;ezplot(t,ft3);信号编码对[11011101001]进行编码。clear;clc;c=[11011101001]fori=1:length(c)ifi==1d1(i)=0;d2(i)=0;elseifi==2d1(i)=c(i-1)
3、;d2(i)=c(i-1);elseifi==3d1(i)=mod(c(i-1)+c(i-2),2);d2(i)=c(i-1);elsed1(i)=mod(c(i-1)+c(i-2),2);d2(i)=mod(c(i-1)+c(i-3),2);endendd1d2迭代法使用一般迭代法求解方程的解。第一根:clear;clc;x0=10;err=1;while(err>10^-6)F1=x0;F2=log(3*x0^2);err=abs(F1-F2);x0=F2;endx0第二根:clear;clc;x=1;err=1;while(err>10^-6)F1=x;F2=sqrt(ex
4、p(x)/3);err=abs(F1-F2);x=F2;endx第三根:clear;clc;x=10;err=1;while(err>10^-6)F1=x;F2=-sqrt(exp(x)/3);err=abs(F1-F2);x=F2;endx牛顿迭代法使用牛顿迭代法求解方程的解。symsxfxfx1;fx=2*x^3-4*x^2+3*x-6;fx1=diff(fx)err=1;k=1;x0=1.5;while(err>10^-6&&k<=1000)x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0))x0=x1;err=abs(x0-x1)k=k+1;end
5、kx0牛顿迭代法求解使用牛顿迭代法求解方程的解。第一解:symsxfxfx1;fx=x.^3-sin(x)-12*x+1;fx1=diff(fx);err=1;k=1;x0=-4;while(err>10^-6&&k<=1000)x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0));err=abs(x0-x1);x0=x1;k=k+1;endkx0第二解:symsxfxfx1;fx=x.^3-sin(x)-12*x+1;fx1=diff(fx);err=1;k=1;x0=0;while(err>10^-6)&&(k<=1000)x1=x0-(subs(fx
6、,x,x0))/(subs(fx1,x,x0));err=abs(x0-x1);x0=x1;k=k+1;endkx0第三解:symsxfxfx1;fx=x.^3-sin(x)-12*x+1;fx1=diff(fx);err=1;k=1;x0=4;while(err>10^-6)&&(k<=1000)x1=x0-(subs(fx,x,x0))/(subs(fx1,x,x0));err=abs(x0-x1);x0=x1;k=k+1;endkx0普通迭代法求解使用普通迭代法求解方程的解。clear;clc;err1=1;err2=1;err3=1;x01=3;x02=0;x03=-3;
7、k=3;while(err1>10^-6)&&(err2>10^-6)&&(err3>10^-6)&&(k<1000)F11=x01;F12=(sin(x01)+12*x01-1)^(1/3);err1=abs(F11-F12);x01=F12;%第一解F21=x02;F22=(x02^3-sin(x02)+1)/12;err2=abs(F21-F22);x02=F22;%第二解F31=x03;F32=-sqrt((sin(x03)-1)/x03+12);err3=abs(F31