测控软件技术基础周鑫资料.docx
- 文档编号:8999222
- 上传时间:2023-02-02
- 格式:DOCX
- 页数:15
- 大小:66.53KB
测控软件技术基础周鑫资料.docx
《测控软件技术基础周鑫资料.docx》由会员分享,可在线阅读,更多相关《测控软件技术基础周鑫资料.docx(15页珍藏版)》请在冰豆网上搜索。
测控软件技术基础周鑫资料
测控软件技术基础
姓名:
周鑫
班级:
232131
学号:
20131000982
指导老师:
黄玉金
BasicSyntaxandCommand-LineExercises
1.Createavectoroftheevenwholenumbersbetween31and75.
代码:
x=32:
2:
75
2.Letx=[2516].
a.Add16toeachelement
b.Add3tojusttheodd-indexelements
c.Computethesquarerootofeachelement
d.Computethesquareofeachelement
代码:
a=x+16;
b=x(1:
2:
end)+3
c=sqrt(x)
d=x.^2
5.Createavectorxwiththeelements...
代码:
a=2:
2:
20
b=10:
-2:
-4
c1=1:
5,c2=1./c1,rats(c2)
d1=0:
4,d2=1:
5,d3=d1./d2,rats(d3)
6.Createavectorxwiththeelements,
xn=(-1)n+1/(2n-1)
Adduptheelementsoftheversionofthisvectorthathas100elements.
代码:
n=1:
100;
x=((-1).^(n+1))./(2*n-1);
y=sum(x)
10.Makeagoodplot(i.e.,anon-choppyplot)ofthefunction
f(x)=sin(1/x)
for0.01 good? 代码: x=0.01: 0.0001: 0.1;y=sin(1./x);plot(x,y) BasicArraySyntaxandManipulations 1.Givenx=[3157926],explainwhatthefollowingcommands"mean"by bysummarizingthenetresultofthecommand. a.x(3) b.x(1: 7) c.x(1: end) d.x(1: end-1) e.x(6: -2: 1) f.x([16211]) g.sum(x) 答: x(3)表示x行向量中的第三个元素。 x(1: 7)表示x行向量中第1到7个元素。 x(1: end)表示x行向量中第1到最后一个元素。 x(1: end-1)表示x行向量中第1到倒数第二个元素。 x(6: -2: 1)表示x行向量中第6个元素到第1个元素中的偶数项元素。 x([16211]): 方括号中表示x行向量中元素的序号,整个表示由这些元素组成的新行向量。 sum(x): 表示x中所有元素求和。 2.GiventhearrayA=[241;672;359],providethecommandsneededto a.assignthefirstrowofAtoavectorcalledx1 b.assignthelast2rowsofAtoanarraycalledy c.computethesumoverthecolumnsofA d.computethesumovertherowsofA e.computethestandarderrorofthemeanofeachcolumnofA(NB.thestandard errorofthemeanisdefinedasthestandarddeviationdividedbythe squarerootofthenumberofelementsusedtocomputethemean.) 代码: A=[241;672;359] x1=A(1,: ) y=A(end-1: end,: ) c=sum(A) d=sum(A,2) N=size(A,1),e=std(A)/sqrt(N) 6.GivethefollowingcommandstocreateanarraycalledF: >>randn('seed',123456789) >>F=randn(5,10); a.Computethemeanofeachcolumnandassigntheresultstotheelementsofa vectorcalledavg. b.Computethestandarddeviationofeachcolumnandassigntheresultstothe elementsofavectorcalleds. c.Computethevectoroft-scoresthattestthehypothesisthatthemeanofeach columnisnodifferentfromzero. d.IfPr(|t|>2.132)=0.1with4degreesoffreedom,areanyofthemeanvalues inthevectoravgstatisticallydifferentfrom0? 代码: randn('seed',123456789) F=randn(5,10); N=size(F,1) avg=mean(F) s=std(F) tscore=(avg-0)./(s/sqrt(N)) ExercisesonRelationalandLogicalOperations 2.Theexerciseshereshowthetechniquesoflogical-indexing(indexingwith 0-1vectors).Givenx=1: 10andy=[3156829470],executeand interprettheresultsofthefollowingcommands: a.(x>3)&(x<8) b.x(x>5) c.y(x<=4) d.x((x<2)|(x>=8)) e.y((x<2)|(x>=8)) f.x(y<0) 输入代码: x=1: 10;y=[3156829470]; 运行结果: a: ans= 0001111000 b: ans= 678910 c: ans= 3156 d: ans= 18910 e: ans= 3470 f: ans= Emptymatrix: 1-by-0 4.Givenx=[315912-10-12961],providethecommand(s)thatwill a....setthevaluesofxthatarepositivetozero b....setvaluesthataremultiplesof3to3(remwillhelphere) c....multiplythevaluesofxthatareevenby5 d....extractthevaluesofxthataregreaterthan10intoavectorcalledy e....setthevaluesinxthatarelessthanthemeantozero f....setthevaluesinxthatareabovethemeantotheirdifferencefromthemean 代码: x=[315912-10-12961] a=x,idxa=x>0,a(idxa)=0 b=x,idxb=~rem(x,3),b(idxb)=3 c=x,idxc=~rem(x,2),c(idxc)=5*c(idxc) 5.Createthevectorx=randperm(35)andthenevaluatethefollowingfunctionusing onlylogicalindexing: y(x)=2ifx<6 =x-4if6<=x<20 =36-xif20<=x<=35 Youcancheckyouranswerbyplottingyvs.xwithsymbols.Thecurveshouldbe atriangularshape,alwaysabovezeroandwithamaximumof16.Itmightalsobe usefultotrysettingxto1: 35.Usingmultiplesteps(orasimpleMfile)is recommendedforthisproblem. 代码: x=1: 35; y=zeros(size(x)); idx1=x<6; idx2=(x>=6)&(x<20); idx3=(x>=20)&(x<=35); y(idx1)=2; y(idx2)=x(idx2)-4; y(idx3)=36-x(idx3); disp([x(: )idx1(: )idx2(: )idx3(: )y(: )]) plot(x,y,'o') ControlofFlow: if-blocks 1.ifn>1a.n=7m=? m=n+1b.n=0m=? elsec.n=-10m=? m=n-1 end 结果: a: m=8; b: m=-1; c: m=-11 2.ifz<5a.z=1w=? w=2*zb.z=9w=? elseifz<10c.z=60w=? w=9-zd.z=200w=? elseifz<100 w=sqrt(z) else w=z end 结果: a: w=2; b: w=0; c: w=7.7460; d: w=200 3.ifT<30a.T=50h=? h=2*T+1b.T=15h=? elseifT<10c.T=0h=? h=T-2 else h=0 end 结果: a: h=0; b: h=31; c: h=1; 5.h(T)=T-10when0 =0.45T+900whenT>100 Testcases: a.T=5,h=-5 b.T=110,h=949.5 结果: 代码: if(T>0&&T<100) h(T)=T-10 elseif(T>100) h(T)=0.45*T+900 end a: ans= -5 b: ans= 949.5000 6.f(x)=-1ifx<0 =0ifx=0 =1ifx>0 CompareyourresultstotheMATLABfunctionsign. 代码: 定义函数: functiony=ckyjjs(x) [m,n]=size(x); y=zeros(size(x)); fora=1: m forb=1: n if(x(a,b)>0) y(a,b)=1; elseif(x(a,b)==0) y(a,b)=0; elseif(x(a,b)<0) y(a,b)=-1; end end end end LoopConstructs 1.Giventhevectorx=[183901],createashortsetofcommandsthatwill a.Addupthevaluesoftheelements(Checkwithsum.) b.Computestherunningsum(forelementj,therunningsumisthesumofthe elementsfrom1toj,inclusive.Checkwithcumsum.) c.computesthesineofthegivenx-values(shouldbeavector) 代码: 1.x=[183901] a: total=0; forj=1: length(x) total=total+x(j); end b: runningTotal=zeros(size(x)); runningTotal (1)=x (1); forj=2: length(x) runningTotal(j)=runningTotal(j-1)+x(j); end c: s=zeros(size(x)); forj=1: length(x) s(j)=sin(x(j)); end 2.CreateanM-by-Narrayofrandomnumbers(userand).Movethroughthe array,elementbyelement,andsetanyvaluethatislessthan0.2to0andany valuethatisgreaterthan(orequalto)0.2to1. 代码: A=rand(4,7); [M,N]=size(A); forj=1: M fork=1: N ifA(j,k)<0.2 A(j,k)=0; else A(j,k)=1; end end end 3.Givenx=[416]andy=[627],computethefollowingarrays a.aij=xiyj b.bij=xi/yj c.ci=xiyi,thenadduptheelementsofc. d.dij=xi/(2+xi+yj) e.eij=reciprocalofthelesserofxiandyj 代码: x=[416],y=[627] N=length(x); forj=1: N c(j)=x(j)*y(j); fork=1: N a(j,k)=x(j)*y(k); b(j,k)=x(j)/y(k); d(j,k)=x(j)/(2+x(j)+y(k)); e(j,k)=1/min(x(j),y(k)); end end c=sum(c) ProgrammingExercises 2.TheFibonaccinumbersarecomutedaccordingtothefollowingrelation: Fn=Fn-1+Fn-2 withF0=F1=1. a.Computethefirst10Fibonaccinumbers. b.Forthefirst50Fibonaccinumbers,computetheratio Fn/Fn-1 Itisclaimedthatthisratioapproachesthevalueofthegoldenmean ((1+sqrt(5))/2).Whatdoyourresultsshow? 代码: a: F (1)=1; F (2)=1; fora=3: 10 F(a)=F(a-1)+F(a-2); end F b: F (1)=1; F (2)=1; fora=3: 50 F(a)=F(a-1)+F(a-2); end F(50)/F(49) 结果: ans= 1.6180 3.TheLegendrepolynomials(Pn(x))aredefinedbythefollowingrecurrancerelation (n+1)Pn+1(x)-(2n+1)xPn(x)+nPn-1(x)=0 withP0(x)=1,P1(x)=xandP2(x)=(3x2-1)/2.Computethenextthree Legendrepolynomialsandplotall6overtheinterval[-1,1]. FormoreinformationonLegendrepolynomials,seetheMathworldsite. 代码: x=-1: 0.1: 1; [a,b]=size(x) fori=1: b P0(i)=1;P1(i)=x(i);P2(i)=(3*x(i).^2-1)/2;P3(i)=(5*x(i).*P2(i)-2.*P1(i))/3;P4(i)=(7*x(i).*P3(i)-3.*P2(i))/4;P5(i)=(9*x(i).*P4(i)-4.*P3(i))/5; end plot(x,P0,x,P1,x,P2,x,P3,x,P4,x,P5) 结果: 12.FollowthedirectionsforExercise11butcreateafunctionthatcomputesthe cumulativesumoftheelementsofavector.Theelementsofthecumulativesum vectoraredefinedby sj=x1+x2+...+xj forj=1: lengthofthevectorx. Thebuilt-infunctionssumandcumsumshouldbeusedinsteadofprodandcumprod,respectively. 代码: functiony=cumulative_sum(x) y=0; a=length(x); s (1)=x (1); fori=2: a s(i)=s(i-1)+x(i); end forj=1: a y=y+s(j); end 举例: x=[111111];cumulative_sum(x) ans= 21
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 测控 软件技术 基础 资料