1、 -5.8845 5.3549 -0.2028 6.3554 -5.6596 4.3293 -6.6325 6.2078 9.0241B/A 30.2855 19.2643 -8.5614 33.4394 21.1547 -9.3974 -0.7443 -0.3938 0.2830(4)A,B 3 54 2 1 -2 67 34 -45 7 2 8 74 87 90 15 9 3 0 A(1,3,:);B2 3 54 2 87 90 15 600 183 -81 684 282 726 15 6 8255.a=1+2i;b=3+4i;c=exp(pi*i/6)c = 0.8660 + 0.50
2、00ic+a*b/(a+b) 1.6353 + 1.8462i第2章3.s=0;for n=0:63 s=s+2n;enddisp(s)n=0:63;s=sum(2.n)1.x=input(输入一个四位整数:y=fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)z=mod(y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)2. gh=input(输入工号h=input(输入工时dj=84;if h120 gz=dj*120+1.15*dj*(h-120);elseif h60 gz=dj*h-70
3、0;else gz=dj*h;format bank;display(gh,gz)3. 循环结构n=input(input n:s=0;for k=1:n s=s+1/k2;display(sqrt(s*6)向量运算k=1:n;display(sqrt(sum(1./k.2)*6)4.y=0;k=0;while y=1e-5 & k500 x0=x; x=a/(b+x0);display(k,x);display(-b+sqrt(b2+4*a)/2,(-b-sqrt(b2+4*a)/2);6.y=fun(40)/(fun(30)+fun(20)(1)函数文件fun.mfunction f=fu
4、n(n)f=n+log(n2+5);(2)函数文件fun.ma=1:f=sum(a.*(a+1);第3章4. t=0:0.01:1.5;y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);x=-10:10;y=linspace(-6,6,size(x,2)z=x.3+3*x.*y.2;plot3(x,y,z)6. x=100:100:400;y=100:z=636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334;X,Y=meshgrid(x,y);mesh(X,Y,z)1. (1)x=
5、-10:plot(x,x-x.3/3/2)(2)plot(x,exp(-x.*x/2)/2/pi)(3)x=-8:8;plot(x,sqrt(64-x.*x)/2)(4)t=0:pi/100:6*pi;plot(t.*sin(t),t.*cos(t)2. (1)x1=linspace(0,1,100);y1=2*x1-0.5;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);plot(x1,y1,r-,x,y,b:text(0.4,1,y=2x-0.5text(-0.4,0.5,x=sin(3t)cos(t)text(-0.4
6、,0.3,y=sin(3t)sin(t)(2) subplot(1,2,1);scatter(x1,y1,10) title(subplot(1,2,2);scatter(x,y,10) subplot(1,21); x=1:1:100; y=sin(1./x);plot(x,y)fplot(sin(1/x),1,100)subplot(2,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3);stem(t,y)subplot(2,2,4);semilogy(t,y);theta=linspace(-pi,pi,100);ro=5.*c
7、os(theta)+4;polar(theta,ro);fi=linspace(0,2*pi,100);a=1r=a.*(1+cos(fi);polar(fi,r);pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);1;x=t;y=t.2;z=t.3;7. x=-30:0;y=0:30;x,y=meshgrid(x,y);z=10.*sin(sqrt(x.2+y.2)./sqrt(1+x.2+y.2);meshc(x,y,z);8. x=linspace(-3,3,100);y=linspace(-3
8、,3,100);x y=meshgrid(x,y);fxy=-5./(1+x.2+y.2);i=find(abs(x)=0.8 & abs(y)0.5)/size(A)h=466,715,950,1422,1635;w=7.04,4.28,3.40,2.52,2.13;hh=500,900,1500;ww=interp1(h,w,hh,splinex=linspace(1,10,50);y=log(x);f=polyfit(x,y,5);yy=polyval(f,x);plot(x,y,x,yy,g.N=64; % 采样点数T=5; % 采样时间终点t=linspace(0,T,N); % 给
9、出N个采样时间ti(I=1:N)x=3*exp(-t); % 求各采样点样本值xdt=t(2)-t(1); % 采样周期f=1/dt; % 采样频率(Hz)X=fft(x); % 计算x的快速傅立叶变换XF=X(1:N/2+1); % F(k)=X(k)(k=1:N/2+1)f=f*(0:N/2)/N; % 使频率轴f从零开始plot(f,abs(F),-*) % 绘制振幅-频率图xlabel(Frequencyylabel(|F(k)|(1)p1=1 2 0 0 7;p2=1 -2;p3=1 0 0 5 1;p12=conv(p1,p2);p=p12+zeros(1,size(p12,2)-
10、size(p3,2),p3;roots(p)(2)A=-1,4,3;2,1,5;0,5,6;Px=polyval(p,A)Pxm=polyvalm(p,A)(1) z=fzero(3*x-sin(x)+1,0)(2) 建立函数文件myfun.mfunction F=myfun(X)x=X(1);y=X(2);F(1)=x*x+y*y-9;F(2)=x+y-1;在命令窗口中输入以下命令:x=fsolve(myfun,3,0,optimset(Display,off)第6章fx=inline(1./(1+x.2)I,n=quad(fx,-100000,100000,1e-3);I,n=quadl(
11、fx,-100000,100000,1e-3);x=-100000:100000;y=1./(1+x.*x);trapz(x,y);(1) fx=inline(-2*y+2*x*x+2*xt,y=ode23(fx,0,0.5,1)(2) fx=inline(y-exp(x)*cos(x)t,y=ode23(fx,0,3,1)for x=1: fx=x,x2,x3;1,2*x,3*x;0,2,6*x; diff(fx)(1)x=0:y=x.10+10.x+1./log10(x);dy=diff(y)/0.01;(2)x=0:y=log(1+x);dy=diff(y,2)/0.01;plot(x(
12、1:99),dy)(1)fx=inline(x.2.*sqrt(2*x.*x+3) quad(fx,1,5)(2)fx=inline(x./sin(x).2quad(fx,pi/4,pi/3)(3)fx=inline(abs(cos(x+y)dblquad(fx,0,pi,0,pi)(4)syms x y;fx=x*y;int(int(fx,y2,y+2),-1,2)x的积分区间为【0,2】时x.*ydblquad(fx,0,2,-1,2)x=0.3:0.2:y=0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325;trapz(x,y)(1)yp=i
13、nline(-(1.2+sin(10*x)*y t,y=ode23(yp,0,5,1);(2)令,则可写出原方程的状态方程形式:-建立函数文件ztfun.mfunction y=ztfun(t,x)b=0;cos(t);y=0,1,0;0,0,1;-1/(3+sin(t),-1,5*cos(2*t)/(t+1)2*x+b;解微分方程组t,y=ode23(ztfun,0,5,1;2);function yy=ztfun(t,y)yy=y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2);t,y=ode23(ztfun,0,5,0;1)第7章(1)数值积分exp(x).*(1+
14、exp(x).2quad(fx,0,log(2)符号积分f=sym(exp(x)*(1+exp(x)2v=int(f,0,log(2);eval(v)(2)略A=sym(1,2,3;x,y,z;3,2,1rank(A)inv(A)det(A)(1) y=sym(sqrt(x+sqrt(x+sqrt(x)y1=diff(y)y2=diff(y,x,2)(2) syms x y;fxy=sin(x2*y)*exp(-x2-y);diff(diff(fxy,x),y)(1) syms xint(1/(1+x4)(2) syms x tint(-2*x*x+1)/(2*x*x-3*x+1)2,x,co
15、s(t),exp(2*t)syms n xsymsum(1/(2*n+1)/(2*x+1)(2*n+1),n,0,inf)symsum(1/(2*n+1)/(2*x+1)(2*n+1),n,0,5)taylor(exp(x)+exp(-x)/2,5,0)(2) syms a xtaylor(exp(-5*x)*sin(3*x+pi/3),5,a)(1)x=solve(sym(x3+a*x+1=0(2) x y=solve(sym(sqrt(x2+y2)-100=0,3*x+5*y-8=07. 方程转化为:符号解y1,y11=dsolve(Dy=y1,Dy1+y=1-t2/piy(-2)=5,y
16、1(-2)=5t数值解编写函数文件ztfun.myy=y(2);1-t2/pi-y(1);在命令窗口输入以下命令t,y=ode45(ztfun,-2,7,-5;5);t=linspace(-2,7,49)y2=y x,y=dsolve(Dx=3*x+4*y,Dy=-4*x+3*yx(0)=0,y(0)=1第9章(1)新建一个Blank GUI。(2)在窗口上添加一个Static Text、一个Edit Text,一个Push Botton,两个Check Box,一组Radio Button。各控件的属性设置如表9.1所示。表9.1 界面控件属性控 件属 性 名属 性 值静态文本框String
17、实验一:编辑框按钮设置字体大小单选按钮1red单选按钮2yellow单选按钮3blue单选按钮4blackValue1.0复选框1倾斜复选框2加粗(3)在按钮的Callback函数中添加以下代码:s=eval(get(handles.edit1,String);set(handles.text1,FontSize,s);(4)在复选框1的Callback函数中添加以下代码:if (get(hObject,Value) = get(hObject,Max set(handles.text1,FontAngleitalicnormal在复选框2的Callback函数中添加以下代码:FontWeig
18、htbold(5)在单选按钮组的SelectionChangeFcn函数中添加以下代码:ForegroundColor,get(eventdata.NewValue,string(2)用菜单编辑器编辑菜单,各菜单项属性如表9.2.表9.2 菜单项属性菜单项Label 属 性 值Tag 属 性 值MenuLineLabelm_lineHelpm_helpLine Menu ItemLine Stylem_linestyleLine Widthm_linewidthLine Colorm linecolorHelp Menu ItemAboutm_about(3)绘图。在窗口的OpeningFcn
19、函数添加以下代码:global h;6;h=plot(t.*sin(t),t.*cos(t)(4)在Line Style菜单项的Callback函数中添加以下代码:set(h,LineStyle-.在Line Width菜单项的Callback函数中添加以下代码:LineWidth,2);在Line Color菜单项的Callback函数中添加以下代码:colorred在About菜单项的Callback函数中添加以下代码:msgbox(Help Itema=2;b=3;r=a+b*t;h_draw=polar(t,r,h_menu=uicontextmenu; %uimenu(h_menu,LabelCallset(h_draw,Col
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1