大学物理学-计算物理-编程练习题答案.docx

clear allclose allfor t=1:50 radius(t)=t/10; r=0.1*t; dcm=D:Q因子; rds=num2str(r); sav_name=dcm rds um.txt; fid(1) = fopen(sav_name,rt); for w=1:50 Q_real(w,t)=fscanf(fid(1),%e,1); end fclose(fid(1);endQ=Q_real;for t=1:50 radius(t)=t/10; r=0.1*t; dcm=D:P因子; rds=num2str(r); sav_name=dcm rds um.txt; fid(1) = fopen(sav_name,rt); for w=1:180 Q_real(w,t)=fscanf(fid(1),%e,1); end fclose(fid(1);endP=Q_real;%角度与波长的选择%P(r,theta) Q(r,lamda);%-初始化种群for g=1:120 %120个种群 angle_num(g)=ceil(50*rand(1); wave_num(g)=50-angle_num(g); for a=1:angle_num(g) try_1(g,a)=ceil(180*rand(1); end for a=1:wave_num(g) try_1(g,a+angle_num(g)=ceil(50*rand(1)+200; endendminn=1000000000;%-for u=1:800for g=1:120 for a=1:50 if try_1(g,a)50 corss_times=cross_times-50; end for i=1:cross_times cross_serial=ceil(50*2*rand(1); if cross_serial50 cross_serial=cross_serial-50; end carrier=try_1(g*2-1,cross_serial); try_1(g*2-1 ,cross_serial)=try_1(g*2,cross_serial); try_1(g*2, cross_serial)=carrier; end end%改变组与组之间的顺序 for s=1:50 carrier_arry(s)=try_1(1,s); end for g=1:120 if g120 for s=1:50 try_1(g,s)=try_1(g+1 ,s); end end if g=120 for s=1:50 try_1(g,s)= carrier_arry(s); end end enduif min_cond_tryminn minn=min_cond_try; for s=1:50 num_need(s)=try_1(min_cond_try_order,s); end end min_cond=minn max_cond= max_cond_try if min_cond=max_cond ballode break endendnum_needmin_condOf course, there are many things to learn if I were enrolled into my ideal research field. First of all, I hope I can build up a systematic view of physics, especially in the field of optics. On the other hand, Im dreaming of the participation of some relevant projects of my major. In that way, I can get more experiences in practice. And what is more, if it is possible, I hope I can continue to study for a doctorate degree. In a word, Im looking forward to making up a solid foundation for my future profession.光学当中我认为最前沿的一块就是我的