《中南大学matlab课后习题2.pdf》由会员分享,可在线阅读,更多相关《中南大学matlab课后习题2.pdf(42页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精 选 试 题 中 南 大 学 matlab课 后 习 题 编 制:_审 核:_出 版:_第 二 章1 求 下 列 表 达 式 的 值。(1)w=sqrt(2)*(l+*10A-6)w=(2)a=;b=5;c=;x=(2*pi*a+(c+b)/(pi+a*b*c)-exp(2)/(tan(b+c)+a);x=(3)a=;b=;y=2*pi*aA2*(l-pi/4)*b-4)*a;y=(4)t=2,l-3i;5J;z=l/2*exp(2*t)*log(t+sqrt(l+tA2);z=+004*2,已 知 a,b,求 下 列 表 达 式 的 值。a=-l,5,-4;0,7,8;3,61,7;b=8
2、3rl;2,5,3;-3A 0;(1)a+6*bans=47 23-1012 37 26-15 73 7aA2-b+eye(3)ans=-18-217 1722 533 10921 867 526(2)a*bans=14 14 16-10 51 21125 328 180a.*bans=-8 15 40 35 24-9 122 0b*aans=-11 0-157 228 533-1 28(3)a/bans=(4)a,bans=-1 5-4 8 3-10 7 8 2 5 33 61 7-3 2 0a(l,3,:);bA2ans=-1 5-43 61 773 37 117 37 13-20 1 9
3、3.已 知 a,完 成 下 列 操 作。a=23,10”0;41 45,65,5;32,5,0,32;6,54J;(1)输 出 a 在 10,25范 围 内 的 全 部 元 素。k=find(a10&a25)a(k)k=1ans=23 取 出 a 前 3行 构 成 矩 阵 b,前 两 列 构 成 矩 阵 c,右 下 角 3*2子 矩 阵 构 成 矩 阵 d,b与 c 的 乘 积 构 成 矩 阵 e。b=a(l:3,:)b=00c=a(:,l:2)d=a(2:4,3:4)d=0e=b*ce=+003*(3)分 别 求 表 达 式 额 ed,e&d,e|d和 0k(1的 值。edans=00010
4、1e&dans=0e|dans=1 11 11 1-e|-dans=0 01 00 04,产 生 5 阶 希 尔 伯 特 矩 阵 H 和 5阶 帕 斯 卡 矩 阵 P,且 求 其 行 列 式 的 值 H h和 H p以 及 他 们 的 条 件 数 T h和 T p,判 断 哪 个 矩 阵 性 能 更 好,为 什 么 format ratH=hilb H=1 1/2 1/3 1/4 1/51/2 1/3 1/4 1/5 1/61/3 1/4 1/5 1/6 1/71/4 1/5 1/6 1/7 1/81/5 1/6 1/7 1/8 1/9P=pascal(5)P=I 1 1 1 11 2 3 4
5、 51 3 6 10 151 4 10 20 351 5 15 35 70Hh=det(H)Hh=1/0000Hp=det(P)Hp=Th=cond(H)Th=476607Tp=cond(P)Tp=178868/21P矩 阵 性 能 更 好,因 为 T p更 接 近 1.5.已 知 A,求 A 的 特 征 值 及 特 征 向 量,并 分 析 其 数 学 意 义。A=-29,6,18;20,5,12;-8,8,5A=-29 6 1820 5 12-8 8 5V,D=eig(A,nobalance,)V=1 361/1013 493/1574-1227/1438-1 1290/593 551/78
6、8 331/713D=-3595/142 0 00-3755/357 00 0 4697/279第 二 早 1.从 键 盘 输 入 一 个 3位 整 数,将 它 反 向 输 出。如 输 入 6 3 9,输 出 为 936.clearm=inputC请 输 入 一 个 三 位 数:,;ml=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);m=nil+m2*10+m3*100;disp(m);请 输 入 一 个 三 位 数:4533542.输 入 一 个 百 分 制 成 绩,要 求 输 出 成 绩 等 级 A、B、C、D、E o其 中 90 100分 为
7、A,80 89分 为 B,70 7 9分 为 C,60 6 9分 为 D,6 0分 以 下 为 E。要 求:(1)分 别 用 i f语 句 和 switch语 句 实 现。(2)输 入 百 分 制 成 绩 后 要 判 断 该 成 绩 的 合 理 性,对 不 合 理 的 成 绩 应 输 出 出 错 信 息。c=inputC请 输 入 一 个 百 分 制 分 数:,;if c=90&c=80&c=70&c=60&c=0&c60disp(E);elsedispC输 入 出 错 力 end请 输 入 一 个 百 分 制 分 数:88Bc=inputC请 输 入 一 个 百 分 制 分 数:,);swi
8、tch fix(c/10)%用 switch 语 句 实 现 case9,10m=A;case8m=B1;case7m=C;case6m=D;casenum2cell(0:5)m=E;otherwise输 入 出 错,;enddisp(m);请 输 入 一 个 百 分 制 分 数:98A3.输 入 2 0个 数,求 其 中 最 大 数 和 最 小 数。要 求 分 别 用 循 环 结 构 和 调 用 MATLAB的 m ax函 数、m in函 数 来 实 现。a=input。请 输 入 二 十 个 数 上 用 循 环 结 构 实 现 max=a(l);min=a(l);for i=l:20;if
9、 maxa(i);min=a(i);endendmaxmin请 输 入 二 十 个 数 3 4 6 5 3 44 5 5 5 6 6 5 5 4 3 32 2 2 2 6max=44min=24.丫=-$1111+0,3)+11#答,当 2取、一、时,求 各 点 的 函 数 值。要 求 分 别 用 顺 序 结 构 和 循 环 结 构 实 现。clearfor a=:%用 循 环 结 构 实 现 y=(exp*a)-exp*a)/2)*sin(a+log(+a)/2)endy=+y=+y=+y=+y=+y=+y=+y=+y=+y=+y=4-y=+y+y=+y=+y=+y=+y=+y=+y=+y=
10、+y=y+y=+y=+y=+y=+y=y=y=yy=yyy=y=y=yy=y=y=y=y=y=y=y=y=y=y=y=yyyyy=y=yy=y=y=y=5.当 n 分 别 取 100、1000、10000时,求 下 列 各 式 的 值:clearn=inputt输 入 一 个 数:,);用 循 环 结 构 实 现 a=0;for i=l:na=a+l/(iA2);enddisp(a);j=l:n;%用 sum函 数 实 现 b=sum(l./(j.A2);输 入 一 个 数:34b=-(2舒 2x(40*4(6二*67)?(,记)?今 clearn=input(,输 入 一 个 数:,);用
11、循 环 结 构 实 现 a=l;for i=l:na=(2*i*i*2)/(2*i-l)*(2*i+l)*a;enddisp(a);输 入 一 个 数:86.建 立 6 5矩 阵,要 求 输 出 矩 阵 第 n 行 元 素。当 n 值 超 过 矩 阵 的 行 数 时,自 动 转 为 输 出 矩 阵 最 后 一 行 元 素,并 给 出 出 错 信 息。cleara=rand(5,6);n=inputC请 输 入 一 个 数:,;if n5n=5a(n,:)dispC输 入 超 出 范 围 了 elsea(n,:)end请 输 入 一 个 数:4ans=请 输 入 一 个 数:75ans=输 入
12、超 出 范 围7.已 天 口 y-式 刈+侬)当 f(n)=n+101n(nA2+5)B,y 的 值 是 多 少。functiona=n2a(n)a=n+10*log(n*n+5)y=n2a(40)/(n2a(30)+n2a(20)a=y=(2)当 f(n)=l*2+2*3+3*4+?+n*(n+l)时,y 的 值 是 多 少。functionb=nlb(n)b=0;fori=l:nb=b+i*(i+l);endy=nlb(40)/(nlb(30)+nlb(20)y=第 四 章 i.绘 制 下 列 曲 线。y=x-x=-10:l();y=x-x.A3/6;plot(x,y)(2*+2y2=64
13、x=-10:10;ezplot(,xA2+2*yA2-64-8,8);grid on;2,设 iy=7T7T,-n tW n,在 同 一 图 形 窗 口 采 用 子 图 的 形 式 绘 制 条 形 图,阶 梯 图,杆 图 和 对 数 坐 标 图 等 不 同 图 形,并 对 不 同 图 形 加 标 注 说 明。t=-pi:pi/10:pi;y=l./(l+exp(-t);subplot(2,2,l);bar(t,y);titled 条 形 图(t,y),);axis(-pi,pi,O,lD;subplot(2,2,2);stairs(t,y,b);titled 阶 梯 图(t,y),;axis(
14、-pi,pi,O,l);subplot(2,2,3);stem(t,y,k);titled 杆 图(t,y),;axis(-pi,pi,O,l);subplot(2,2,4);loglog(t,y,y);title。对 数 坐 标 图(t,y);Warning:Negative data ignored In title at 42In title at 233.绘 制 下 列 极 坐 标 图。(1)p=5cos。+4t=0:pi/50:2*pi;r=5.*cos(t)+4;polar(t,r);title(,rho=5*costheta+4,);p=5*cos0+4/与、5 sin2(?兀
15、兀(2)yi=COSQ-3r T-3t=-pi/3:pi/50:pi/3;r=5.*(sin(t).A2)./cos(t);polar(t,r);4.绘 制 下 列 三 维 图 形。(1)x=eA-(t/20)costy=eA-(t/20)sint?0t2jtt=0:pi/50:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t;plot3(x,y,z);grid on;(2)z=5,|x|w5,|ylw5.要 求 应 用 插 值 着 色 处 理。x,y=meshgrid(-5:5);z=zeros(ll)+5;mesh(x,y,z);sh
16、ading interp;6axis off.shading 总 加;=moviein(20);i=l:20axis(U h i,i Dm(M=g“frame.end第 五 章 1,利 用 matlab提 供 的 randn函 数 生 成 符 合 正 态 分 布 的 10*5随 机 矩 阵 A,进 行 以 下 操 作:(1)A 各 列 元 素 的 均 值 和 标 准 方 差。(2)A 的 最 大 元 素 和 最 小 元 素。(3)求 A 每 行 元 素 的 和 以 及 全 部 元 素 的 和。(4)分 别 对 A 的 每 列 元 素 按 升 序,每 行 元 素 按 降 序 排 列。A=rand
17、n(10,5)x=mean(A)y=std(A)Max=max(max(A)Min=min(min(A)Sumhang=sum(A,2)SumA=sum(Sumhang)B=sort(A);C=sort(B,2,descend*);CAX=y=Max=Min=Sumhang=SumA=C=2,按 要 求 对 指 定 函 数 进 行 插 值 和 拟 合。(1)按 表 一 所 示 用 三 次 样 条 方 法 插 值 计 算 0 9 0 内 整 数 点 的 正 弦 值 和 0 7 5 内 整 数 点 的 正 切 值,然 后 用 5次 多 项 式 拟 合 方 法 计 算 相 同 的 函 数 值,并 将
18、 两 种 计 算 结 果 进 行 比 较。(2)按 表 二 所 示 用 3次 多 项 式 方 法 插 值 计 算 1100之 间 整 数 的 平 方 根。表 一 特 殊 角 的 正 弦 与 正 切 值 表 a(度)0 15 30 45 60 75 90sin a 0 07071tan a 0表 二 1100内 特 殊 值 的 平 方 根 表 N I 4 9 16 25 36 49 64 81 100VN 1 2 3 4 5 6 7 8 9 10(1)a=0:15:90;b=a./180.*pi;s=sin(b)c=0:15:75;d=c./180.*pi;t=tan(d)e=inputC请 输
19、 入 想 计 算 的 值:,)S=sin(e/18()*pi)T=tan(e/180*pi)Sl=interpl(a,s,e,Spline1)Tl=interpl(c,t)e/spline,)Pl=polyfit(a,s,5);P2=polyfit(c,t,5);S2=polyval(Pl,e)T2=polyval(P2,e)00请 输 入 想 计 算 的 值:s=T=SI=T1=S2=T2=(2)n=l,9,16,25,36,49,64,81,100;N=sqrt(n);x=inputC计 算 值:*);interpl(iijNyX/cubic*)计 算 值:ans=53,己 知 h=e%t
20、,t2 0,取 N=6 4,对 t 从。5 秒 采 样,用 FFT作 快 速 傅 里 叶 变 换,并 绘 制 相 应 的 振 幅-频 率 图。N=64;T=5;t=linspace(O,T,N);h=exp(-t);dt=U2)t;f=l/dt;X=fft(t);F=X(l:N/2+l);f=f*(0:N/2)/N;plot(f,abs(F),-*t)4,有 两 个 多 项 式 P(x)=2x八 4 3x八 3+5x+13,Q(x)=x八 2+5x+8,试 求 P(x),P(x)Q(x)以 及 P(x)/Q(x)的 导 数。P=2,.3,0,5,13;Q=U,5,8;p=polyder(P)q
21、=polyder(P,Q)a,b=polyder(P,Q)P=8-9 0 5q=12 35 4-57 76 105a=4 27 34-77-26-25b=1 10 41 80 645,有 3 个 多 项 式 PI(x)=xA4+2xA3+4xA2+5,P2(x)=x+2,P3(x)=x2+2x+3,试 进 行 下 歹 操 作:(1)求 P(x)=Pl(x)+P2(x)P3(x)o(2)求 P(x)的 根。(3)当 x 取 矩 阵 A 的 每 一 元 素 时,求 P(x)的 值。其 中 A=-l20 5(4)当 以 矩 阵 A 为 自 变 量 时,求 P(x)的 值。其 中 A 的 值 与(3)
22、相 同。Pl=l,2,4,0,5;P2=0,l,2;P3=l,2,3;P=Pl+conv(P2,P3)X=roots(P)A=-l;,2,;0,5,;p=polyval(P,A)Undefined function or variable*x*.第 六 章 1,下 面 是 一 个 线 性 病 态 方 程 组:1 12 31 13 41 14 5(1)求 方 程 的 解。(2)将 方 程 右 边 向 量 元 素 b3改 为,再 求 解,并 比 较 b3的 变 化 和 解 的 相 对 变 化。(3)计 算 系 数 矩 阵 A 的 条 件 数 并 分 析 结 论。A=l/2,l/3,l/4;l/3,
23、l/4,l/5;l/4,l/5,l/6;P=,J;x=ApA=l/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6J;P=J;x=Apcond(A)x=;xl 0.95-*x2=0.67i x3 0.52x=ans=+0032,求 下 列 方 程 的 解。(1)x4 1+x3+l=0,x 0=-lfunction fx=fun(x)fx=x.A41+x.A3+l;xO=fzero(,fun*,-l)x0=function fx=fun(x)fx=x-sin(x)/x;xOMzeroCfun*,-1)sinx+y2+ln z-7=0(3)3x+2y-z3+1=0,初 值 xO=
24、l,yO=l,zO=l.,x+y+z-5=0function q=myfun(d)x=d(l)y二 d z=d(3)q(1)=sin(x)+yA2+log(z)-7;q(2)=3*x+2Ay-zA3+l;q(3)=x+y+z-5options=optimset(Display,ofr);x=fsolve(myfun,1,1,1,options)4,求 函 数 在 指 定 区 间 的 最 大 值。1+f(x)=7 7 G,x(0,2)function x y=qiujizhi()h=;%解 的 精 度 x=;y=(l+xA2)/(l+xA4);while(l+(x+h)A2)/(l+(x+h)A
25、4)yx=x+h;y=(l+xA2)/(l+xA4);endfprintfC极 大 值 时 x取 值,极 大 值 时 y=,%x,y);极 大 值 时 x 取 值,极 大 值 时 y=,+000ans=5,设 有 4 00万 元 资 金,要 求 4 年 内 使 用 完,若 在 一 年 内 使 用 资 金 x 万 元,则 可 得 效 益 J x 万 元(效 益 不 能 再 使 用),当 年 不 用 的 资 金 可 存 入 银 行,年 利 率 为 10%。试 制 定 出 资 金 的 使 用 计 划,以 使 4 年 效 益 之 和 最 大。(1)先 建 立 M 文 件,定 义 目 标 函 数:fun
26、ctionf=fun4(x)f=-(sqrt(x(1)+sqrt(x(2)+sqrt(x(3)+sqrt(x(4);(2)再 建 立 M 文 件 定 义 非 线 性 约 束:functionfg,ceq=mycon(x)g(l)=x(1)-400;g(2)=*x(1)+x(2)-440;g(3)=*x(1)+*x(2)+x(3)-484;g(4)=*x(1)+*x(2)+*x(3)+x(4);ceq=0(3)主 程 序 为:x 0=l;l;l;l;vlb=0;0;0;0;vub=;A=;b=;Aeq=;beq二 口;x,fval=fmincon(,fun4,xO,A,b,Aeq,beq,vlb
27、,vub,mycon)x=fval=第 七 章 1,求 定 积 分。I=&x(2)(x-0.3)2+0 0 1-(x.0.9)2+0.04-1 1(1)format longfx=inline(sin(x)./x);|I,n=quadl(fx,0,2,le-10)1=n=19(2)fonnat longfx=inline(,l./(.A2+-L/(.A24-6,);I,n=quad(fx,0,l,le-10)1=9852,求 二 重 定 积 分。(1)Ii=/o/oe-&+y)dxdyglobal ki;ki=0;I=dblquad(fxy,O,l,OJ)ki(2)12=J:5|cos?(x+
28、y)|dxdyf=inline(abs(cos(x+y),xy);I=dblquad(f,O,pi,O,pi)3,分 别 用 矩 形,梯 形 公 式 计 算 由 表 一 中 数 据 给 出 的 定 积 分 I=/:;f(x)dx.表 一 被 积 函 数 f(x)数 据 表 k 1 2 3 4 5 6 7Xk 0,7F(Xk)x=:;F=L,J;trapz(X,F)4,设 X 由 2,2 叮 区 间 内 均 匀 分 布 的 10个 点 组 成,求 向 量 sinX的 13阶 差 分。p=0:pi/5:2*pi;for n=l:3nDX=diff(sin(p),n)end5,设SU1X用 3 种
29、不 同 的 方 法 求 函 数 f(x)的 数 值 导 数,并 在 同 一 个 坐 标 系 中 作 出 f(x)的 图 像。f=inline(*sin(x)J(x+cos(2.*x)*);g=inline(,(cos(x).*(x+cos(2*x)-sin(x).*(l-2.*sin(2*x)/(x+cos(2.*x).A2,);x=pi:pi;p=polyfit(x,f(x),5);dp=polyder(p);dpx=polyval(dp,x);%求 dp 在 假 设 点 的 函 数 值 dx=diff(f(x,)/;%直 接 对 f(x)求 数 值 导 数 gx=g(x);%求 函 数 f
30、 的 导 函 数 g 在 假 设 点 的 导 数 plot(x,dpx,x,dx,,x,gx,-);第 八 章 1,分 解 因 式。(1)xA4-yA4syms x y;s=xA4-yA4;factor(s)ans=(x-y)*(x+y)*(xA2+yA2)(2)5135factor(5135)ans=5 13 792,求 函 数 的 极 限。syms x;f=(x-2)/(xA2-4);limit(f9x,2)ans=1/4m-Jarccosx sym x;f=(sqrt(pi)-sqrt(acos(x)/sqrt(x+l);limit(f,x,-l,right,)ans=-Inf3,求 函
31、 数 的 符 号 导 数。(1)y=sinj,求 y;y”。s y m x;f=sin(l/x);diff(f;x2)ans=-cos(l/x)/xA2ans=(2*cos(l/x)/xA3-sin(l/x)/xA41-cos?(2x),x,”y=一 一,求 y,y.s y m x;f=(l-cos(2*x)/x;difT(f,k)diff(f,x,2)ans=(2*sin(2*x)/x+(cos(2*x)-l)/xA2ans=(4*cos(2*x)/x-(4*sin(2*x)/xA2-(2*(cos(2*x)-l)/xA34,求 积 分。(1)J J e,+Idxs y m x;f=sqrt
32、(exp(x)+l);int(f,x*)ans=2*(exp(x)+l)A(l/2)+2*atan(exp(x)+l)A(l/2)*i)*i 舄 dysyms x y;f=x/(x+y);int(f,y)ans=x*log(x+y)(3)/(I+ex)2dxs y m x;f=exp(x)*(l+exp(x)A2;int(f,x,0,log(2)ans=(exp(62433/9007)*(3*exp(62433/9007)+exp(62433/)+3)/3-7/3(4)/:xlnxdxsym x;f=x*log(x);int(f,x,l,exp(l)Undefined function or
33、variable*x*.5,求 下 列 级 数 之 和。(1)1-1/2+1/3-1/4+-+C-l)n+1H-sym x;s=symsum(-l)A(x+l)/x,l,Inf)s=log/C、x3 x5 x7 x+y+y+7+?symy;z=symsum(yA(2*y-l)/(2*y-l),l,Inf)z=sum(yA(2*y-1)/(2*y-1),y=I.Inf)6,求 函 数 在 x=xO处 的 泰 勒 展 开 式。、ex+ex(1)y=-,x0=0,n=5.(2)(7x3-2x+l,xO=0,n=6sym x;fl=(exp(x)+exp(-x)/2;f2=sqrt(xA3-2*x+l
34、);taylor(fl,x,5,0)taylor(f2,x,6,0)ans=xA4/24+xA2/2+1ans=-xA5/8-xA4/8-xA2/2-x+17,求 非 线 性 方 程 的 符 号 解。(1)x3+ax+1=0(2)sinx+2cosx-J x=0(3)ln|=9,ex+y=3syms x y a;x=solve(,xA3+a*x+l=0*,x)x=solve(*sin(x)+2*cos(x)-sqrt(x)=0,/x,)x y=solve(,log(x/y)=9,/exp(x+y)=3,x,/yr)x=(aA3/27+l/4)A(l/2)1/2)人(1/3)a/(3*(aA3/
35、27+1/4)八(1 1/2)八(1/3)a/(6*(aA3/27+l/4)A(l/2)-l/2)A(l/3)-(aA3/27+l/4)A(l/2)-l/2)A(l/3)/2-(3A(l/2)*(a/(3*(aA3/27+1/4)A(1/2)-l/2)A(l/3)+(aA3/27+l/4)A(l/2)-l/2)A(l/3)*i)/2a/(6*(aA3/27+1/4)八(1/2)-1/2)八(1/3)-(a八 3/27+1/4)八(1/2)1/2)八(3)/2+(3A(l/2)*(a/(3*(aA3/27+1/4)八(1/2)-1 八(1/3)+(aA3/27+l/4)A(l/2)-l/2)A(
36、l/3)*i)/2x=*i-227.X=(exp(9)*log(3)/(exp(9)+1)y=log(3)/(exp(9)+1)8,求 微 分 方 程 初 值 问 题 的 符 号 解,并 与 数 值 解 进 行 比 较。xy+(l-n)y,+y=0y(0)=y(0)=0syms n;|x,y=dsolve(,x*(D2y)+(l-n)*(Dy)+y=0,;y(0)=0*Dy(0)=0,x)Warning:possibly missing solutions solviniWarning:Explicit solution could not be found.In dsolve at 101x
37、=empty sym y=第 九 十 二 章 第 九 章 1建 立 一 个 图 形 窗 口,使 之 背 景 颜 色 为 红 色,并 在 窗 口 上 保 留 原 有 的 菜 单 项,而 且 在 按 下 鼠 标 左 键 之 后 显 示 出 Left Button Pressed字 样。hf=figure(,colorl 0 OJ/MenuBar,/figureWindowButtonDownFcn/textj/Teft ButtonPressed*)hf=12,利 用 图 形 对 象 绘 制 下 列 曲 线,要 求 先 利 用 默 认 属 性 绘 制 曲 线,然 后 通 过 图 形 句 柄 操 作
38、 来 改 变 曲 线 的 颜 色,线 型 和 线 宽,并 利 用 文 本 对 象 给 曲 线 添 加 文 字 标 注。y=;ln?(x+71+x2)clear allx=-l:5;y=*log(x+sqrt(l+x.A2);hy=plot(x,y)seUhy/color/r*,TineStyle,f/LineWidth*,1);text(3,2 八 l,y=log(x+sqrt(l+x 八 2)hy=t=0:2:100;x=t.*t;y=5*LA3;hf=plot(x,y);set(hf/color*,b,/linestyle,:,marker,p,linewidth1,;text(5000,
39、3*10A6;x=tA2 y=5tA3f)x 105 r-4.54*3.5*3x=t2 y=5t3*2.52 K1.5it求 0.51000 2000 3000 4000 5000 6000 7000 8000 9000 100003,利 用 图 形 对 象 绘 制 下 列 三 维 图 形,要 求 对 图 形 进 行 光 照 处 理。2 2z=x+y-5sin?(xy)x,y=meshgrid(-1000,10,1000);z=x.*2+y.*2-5*sin(x.*y);plot3(x,y,z);%axis equallight(color0,0,lJ;position,-1000,-1000
40、,2)shading interp;hold on;plot3(-1000,-1000,2,,P,);text(-1000,-1000,2;light,)500.0-500-light-1000-1500-2000-、500O-299-999.5-500-10001000.5-1000-1Q014.以 任 意 位 置 子 图 形 式 绘 制 出 正 弦,余 弦,正 切 和 余 切 函 数 曲 线。t=linspace(-2*pi,2*pi,100);yl=sin(t);y2=cos(t);y3=tan(t);y4=cot(t);axes(position90);axes(*position*,
41、);Iine(t,y2/Iinestyle,*,color,/b,);axes(position,1);line(t,y3,linestyle,,cok)r T);axestpositionM 0);line(t,y,linestyle.,color,g);Warning:Setting marker type based on LineStyle propertyMarker*and LineStyle have been made separate propertiesUse Marker property to remove this warning5,用 patch函 数 绘 制 一
42、个 填 充 渐 变 色 的 正 五 边 形。%clear all;clc;t=linspace(0,2*pi,6);k=t(1:Iength(t)-1);pa tch(sin(k),cos(k)9rand(19kngth(k)9 facecolor/interp1);axis equal1-0.8 0.6 0.4 0.2 0-0.2-0.4-0.6-0.8-1-0.8-0.6-0.4-0.2 0 0.2 0.4 0.6 0.8 1第 十 章 1,设 计 菜 单。菜 单 条 中 含 有 File和 H e l p两 个 菜 单 项。如 果 现 则 File中 的 N e w 命 令,则 将 显
43、示 N e w Item字 样,如 果 选 择 File中 的 O p e n 命 令,则 将 显 示 出 Open Item字 样。File中 的 Save命 令 初 始 时 处 于 禁 选 状 态,在 选 择 H e l p菜 单 项 之 后 将 此 菜 单 项 恢 复 成 可 选 状 态,如 果 选 择 File中 的 S a v e命 令,则 将 显 示 Save Item字 样。如 果 选 择 File中 的 Exit命 令,则 将 关 闭 当 前 窗 口。如 果 选 择 H e l p中 的 About 命 令,则 将 显 示 Help Item字 样,并 将 S a v e菜 单
44、 设 置 成 可 选 状 态。screen=get(O,*ScreenSize,);W=screen(3);H=screen(4);hf=figureCCok)rMl,l,lVPo itioiiMl,l,*W,*HTName;菜 单 设 计 示 例 一,;NumberTitle,off;M e n u B a r,none,);hfile=uimenu(hf,label,File*);hhelp=uimenu(hf,label*,Help*);uimenu(hfile,label1,*&New*,call*,disp(uNevv Item);uimenu(hfile,laber;&O p e
45、n,cair;disp(,Op e n Item);hsave=uimenu(hfile,UaberZ&Save,Enable*,*ofT);uimenu(hsave,*laber/Text file,/call,kl=0;k2=l;file01;*);uimenu(hsave/laber,Graphics file,call,kl=0;k2=0;fHel0;);uimenu(hfiIe/laber/Save&A s,call,disp(MSave As Item);uimenu(hfile,label1,*&Exif,Separator,*011,cair,close(hf),);uimen
46、u(hhelp,label,*About.J/cairjldispfHelp Item,);,set(hsave,uEnable,non,)*);4,设 计 一 个 用 户 界 面,其 中 有 一 个 静 态 文 本 框,一 个 编 辑 框,两 个 复 选 框,一 组 单 选 按 钮。在 编 辑 框 中 输 入 一 个 数,单 击 按 钮,可 以 设 置 静 态 文 本 框 中 的 文 字 大 小;复 选 框 用 于 设 置 文 字 是 否 为 粗 体,倾 斜;单 选 按 钮 用 于 设 置 字 体 的 颜 色。(1)打 开 G U I设 计 窗 口,添 加 有 关 图 形 对 象。在 命 令
47、 输 入 窗 口 输 入 命 令 guide”,打 开 设 计 窗 口。在 窗 口 上 添 加 一 个 按 钮、一 个 静 态 文 本 框、一 个 编 辑 框、两 个 复 选 框、一 组 3 个 单 选 按 钮。(2)设 置 图 形 对 象 的 属 性。设 置 按 钮 的 Siring属 性 为“设 置 字 体 大 小”;设 置 静 态 文 本 框 的 Backgroundcolor属 性 为 白 色、String属 性 为“字 体 大 小 颜 色 粗 体 倾 斜”;设 置 编 辑 框 的 String属 性 为“设 置 第 一 个 单 选 按 钮 的 String属 性 为“红 色”;设 置
48、 第 二 个 单 选 按 钮 的 Siring属 性 为“蓝 色”;设 置 第 三 个 单 选 按 钮 的 String属 性 为“绿 色”;设 置 第 一 个 复 选 框 的 String属 性 为“粗 体”;设 置 第 二 个 复 选 框 的 String属 性 为“倾 斜”。(3)利 用 位 置 调 整 工 具,调 整 控 件 未 知。选 中 3 个 单 选 按 钮、2 个 复 选 框,设 置 为 左 对 齐。(4)保 存 图 形 用 户 界 面。(5)编 写 程 序 点 击 按 钮 的“Callback”,添 加 程 序:a=eval(get,String)set,FontSize,a
49、)点 击“红 色”单 选 按 钮 的“Callback”,添 加 程 序:set,ForegroundColor,red);set,Value,0);set,Value,0);点 击“蓝 色”单 选 按 钮 的“Callback”,添 加 程 序:set,ForegroundColor,blue);set;ValueO);set,Value;。);点 击“绿 色”单 选 按 钮 的“Callback”,添 加 程 序:set/ForegroundColorVg*);set,Value,0);set/Value,0);点 击“粗 体”复 选 框 的“Callback”,添 加 程 序:if(ge
50、t,Value1)=get/Max)set/FontWeight/bold*)elseset/FontWeight/normar)end点 击“倾 斜”复 选 框 的“Callback”,添 加 程 序:if(get(hObject,Value)=get(hObject/Maxr)set,FontAngle,Italic);elseset,FontAngle1,normal1);enduntitled-I E x|体 小 色 体 斜 字 大 颜 粗 侦 r 红 色 6 蓝 色 r 级 色 V 您 住 r 银 斜 12 设 置 字 体 大 小 untitled|_ 口 X第 十 二 章 6 红