tr-1441.p65

Size: px
Start display at page:

Download "tr-1441.p65"

Transcription

1 NAL TR-1441 NAL TR-1441 ISSN UDC TECHNICAL REPORT OF NATIONAL AEROSPACE ACE LABORATOR TORY TR-1441 NATIONAL AEROSPACE ACE LABORATOR TORY OF JAPAN AN

2 1 1 1 Self-illumination Correction for Pressure Sensitive Paints on Airplane Models in Wind Tunnel Testing Hiroki SUGIURA 1, Yoshimi IIJIMA 1, Keisuke ASAI 1 ABSTRACT One of the potential problems in applying pressure sensitive paints to a wing-body model is selfillumination of the model. A theoretical model of self-illumination correction by W. M. Ruyten was applied to two, square, painted samples and a wind tunnel test of the MU-300 model. The spatially averaged reflection coefficient of 0.75 was determined for a pressure sensitive paint PtOEP-GP197 by measurement on the two, square, painted samples at a range of included angles. As a result, the applicability of the self-illumination correction method was validated. The self-illumination effect on the MU-300 model was so small that the effect of the self-illumination correction was little. Keywords :Pressure Sensitive Paint, Self-Illumination Correction, Reflection Coefficient Ruyten 2 PtOEP-GP MU-300 I i J j A D Bi i,j A ij j M S j x i N C B i p n P received 27 November, 2001 (Fluid Science Research Center)

3 R S S T Bi x S i i C B i x ij α γ θ x ij x j x i B C n j x ji x j S i B 2 B C 2 B C Ruyten 6) corr i 2 j ref MU ),2) % 5) Ruyten 6) 2 4 MU nm 3),4) self-illumination 1 Lambert Dornier Alpha jet R

4 3 A ( ni xij)( ni x ji) Ruyten Aij = S 4 π xij BSDF (Bidirectional Scattering Distribution Function) i=j A ij 0 n i x ij 0 n j x ji 0 A ij 0 6)Le Sant x i x j A ij 0 BSDF BSDF 3 7) 3.1 a) 2 b) I p Stern-Volmer 8) I ref p 50mm 2mm 2 = A+ B (2.1) I p ref A,B ref 2 B 545mm 25mm B convergent distortion B 2 90 B B 605mm 16 (diffusive filter) j (2.3) S S i i x i i S corr i corr S i = S i Rred AS ij j (2.2) 6) A A ij j S j x i A ij 2

5 (3.5) (2.1) =100, 120, 140, nm B C B 0 S B S C PtOEP- B C GP197 (2.2) SB( γ) RredA( γ) SC( γ) = S 14 (180 ) B (3.1) CCD 60,000 =180 2 electrons nm B C ij DBi() γ = SBi() γ S (180 ) Bi (3.2) f22 T () γ = () γ () γ Bi Aij SCi (3.3) j (3.1) CCD CCD 1 5 A D Bi T Bi CCD CCD I p 3

6 5 () () γ D γ = R T Bi red Bi T Bi C B i D Bi C B i C =100, 120, 140, : Q = DBi( γ) RredTBi( γ) (3.4) γ i =180 DBi( γ) TBi( γ) =100, 120, 140, 160 γ i Rred = =100, 120, 140, 160 T 2 Bi() γ (3.5) γ i 5 =100, 120, 140 D Bi T Bi = =100, 120, 140,

7 (deg) R Q 100 1/ / / / / / / / ,120,140 1/ ,120,140 1/ ,120,140,160 1/ =100, 120, 140, ,120,140, ,120,140,160 1/ =100, 120, 140 B C BC PtOEP-GP ,120,140,160 C (3.5) R R 75% FIB7 R B B /16 B 6 (a) ,120,140, (b) 0 1/41/16 8 3% Q 30% Q C B 7 (3.4) B =180 B 6 (c) 100,120, /4 Q= m2m

8 7 6

9 B =100, 120, 140, 160 B =180 8MU-300 8% on 2 8(a) off 8(b) 2 2 PtOEP-GP197 MU-3008%

10 m 80 8, 8 3mm CCD 9 14 CCD nm on off f22 M= =0 4 80kPa 308K CCD in situ 1.2m CCD 1 10 CCD In situ

11 M=0.75 =2.5 on 2 C P 11M=0.75 =2.5 on 2 C P

12 11 12M=0.75 =2.5 off 2 C P 13

13 M=0.75 =2.5 on 2 C P C P 11 C P RGB off 2 C P MU-300 C P 12 13% 19% C P M=0.75 =2.5 19% 9 5 2% Ruyten C P 0.03 off 2 PtOEP- GP MU-300 ) 18 4 pp m=0.75 =2.5 19% )

14 13 ( ni xij) Si 31 Ω= 2 xij xij pp ) Y. Shimbo, K. Asai, H. Kanda, Y. Iijima & N. Komatsu: Evaluation of Several Calibration Techniques for Pressure Sensitive Paint in Transonic ( nj x ji) ( ni xij) Si Testing, AIAA Ω IS jcosθ = S j I 2 xij xij xij ) Y. Shimbo, K. Asai, Y. Iijima, H. Sugiura, H. Kanda, (A3) N. Komatsu, S. Kita & M. Ishiguro:Pressure Sensitive Paint Application to a Business Jet Model in B4 I J B Transonic Testing, 7th Pressure Sensitive Paint Workshop, I = J /π (A4) ) W. M. Ruyten:Correcting Luminescent Paint (A3) (A4) S j S i Measurements for Self-illumination, ICIASF '97 ( ni xij)( nj x ji) Record, International Congress on Instrumentation S 4 j SiJ π xij in Aerospace Simulation Facilities, pp.3-9, (A5) ) W. M. Ruyten: Self-illumination calibration technique for luminescent paint measurements, Review (A5) J S j S i S i of Scientific Instruments, vol.68, no.9, pp ( ni xij)( nj x ji) 3457, Aij = S 4 j π xij ) Y. Le Sant:Overview of the Self-Illumination Effect (A6) Applied to Pressure Sensitive Paint Applications, ICIASF 2001 Record, pp , ) R.C. Crites: Pressure Sensitive Paint Technique, Von Karmann Institue for Fluid Dynamics, Lecture Series , Measurement Techniques, pp.1-69, B I J j S j 15 x j S j r P i S j P i I P i ) P = θ θ φ = 2 i rd rsin d r sinθdθdφ 63 (B1) ) 32 pp A A x j x i x ji S j x ji n j x ji ( nj xji) Sjcosθ = S x ji j (A1) i S i S j x ji x j S i (A2) S j J S i I S j S i (A1) 15 S j r P i

15 S j x ji S j cos S j x j P i P i /r 2 S j P i (B1) Pi Sjcosθ I = S sinθcosθ θ φ 2 ji d d (B2) J Sj = Iπ S r I = J/ π S j S j J j (B4) 2 π π/2 Sj Isinθcosθdθdφ = Iπ S j (B3) 0 0 C MATLAB program:calculation OF REFLEX COEFFICIENT OF PSP % by Hiroki Sugiura(NAL),99/12/28. % based on W.M.Ruyten:Rev.Sci.Instrum.68(9),1997/9. clear;tic;xmax=512;ymax=512;nrm=1/4000; sby1=211;sby2=329;;sbx1=215;sbx2=335; B xbnrm=sbx2-sbx1+1;ynrm=sby2-sby1+1; Rup=0;Rdown=0;Dbsqr=0;Zdbi=0;Ztbi=0;scf=1; xbn=xbnrm;yn=ynrm; %Read Dark Image fid=fopen(' matlab reflex T dark1.drk','r'); drk1=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T dark2.drk','r'); drk2=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T dark3.drk','r'); drk3=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T dark4.drk','r'); drk4=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T dark5.drk','r'); drk5=fread(fid,[ymax,xmax],'short');status=fclose(fid); drk=(drk1+drk2+drk3+drk4+drk5)*0.2; clear drk1 drk2 drk3 drk4 drk5; % Read Image SB(180) =180 fid=fopen(' matlab reflex T t180f1.img','r'); temp1=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t180f2.img','r'); temp2=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t180f3.img','r'); temp3=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t180f4.img','r'); temp4=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t180f5.img','r'); temp5=fread(fid,[ymax,xmax],'short');status=fclose(fid); temp=(temp1+temp2+temp3+temp4+temp5)*0.2;

16 15 clear temp1 temp2 temp3 temp4 temp5; temp=temp-drk; at=temp(1:ymax/2+20,2:xmax);bt=temp(ymax/2+21:ymax,1:xmax-1); temp=[bt; at];clear at;clear bt; sb18=temp(sby1:sby2,sbx1:sbx2);clear temp; % Read Image SB(100) and SC(100) =100 v=100;disp(v) fid=fopen(' matlab reflex T t100f1.img','r'); temp1=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t100f2.img','r'); temp2=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t100f3.img','r'); temp3=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t100f4.img','r'); temp4=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t100f5.img','r'); temp5=fread(fid,[ymax,xmax],'short');status=fclose(fid); temp=(temp1+temp2+temp3+temp4+temp5)*0.2; clear temp1 temp2 temp3 temp4 temp5;temp=temp-drk; at=temp(1:ymax/2+20,2:xmax);bt=temp(ymax/2+21:ymax,1:xmax-1); temp=[bt; at];clear at;clear bt; sb10=temp(sby1:sby2,sbx1:sbx2);db10=sb10-sb18;db10=medfilt2(db10); % Resize SC(100) in Y-axis C Y dymax1=5;dymax2=12;sc10y1=sby1-dymax1;sc10y2=sby2+dymax2; djx=-11; sc10x1=round(sbx1-cos((180-v)*pi/180)*(sbx2-sbx1))+djx;sc10x2=sbx1-3; ynrm=sby2-sby1+1;xc10nrm=sc10x2-sc10x1+1; i=0;sc10y=zeros(ynrm,xc10nrm); %imshow(nrm*temp(sc10y1:sc10y2,sc10x1:sc10x2),64); for ix=sc10x1:sc10x2; dy1=round((sc10x2-ix)*dymax1/(xc10nrm-1)); dy2=round((sc10x2-ix)*dymax2/(xc10nrm-1)); i=i+1;sc10y(1:ynrm,i)=imresize(temp(sby1-dy1:sby2+dy2-1,ix),[ynrm,1]); sc10=sc10y;tic;clear temp; % Compress Image for Faster Calculation { } xbnrm=round(xbnrm/scf);ynrm=round(ynrm/scf);xc10nrm=round(xc10nrm/scf); db10=imresize(db10,[ynrm xbnrm]);sc10=imresize(sc10,[ynrm xc10nrm]); % Calculate Matrix A(100),DBi(100),and TBi(100) { A D Bi T Bi } dbi10=reshape(db10,[1 xbnrm*ynrm]);scj10=reshape(sc10,[xc10nrm*ynrm 1]); tan2v=(tan(v*pi/180))^2;tanv=tan(v*pi/180);tan2vp=tan2v/pi; a=zeros(1,ynrm*xc10nrm); tbi10=zeros(1,ynrm*xbnrm); sbrate=(xc10nrm*scf+djx)/(xc10nrm*scf)

17 for ix=1:xbnrm; ixix5t2p=(ix-0.5)*tan2vp; ixyn=(ix-1)*ynrm; for iy=1:ynrm; for jx=1:xc10nrm; jxm=jx*sbrate; ijx2jix2=(ix+jxm-1)^2+(jxm-0.5)^2*tan2v; bunsi=ixix5t2p*(jxm-0.5); jxyn=(jx-1)*ynrm; for jy=1:ynrm; a(jxyn+jy)=bunsi/(ijx2jix2+(iy-jy)^2)^2; tbi10(ixyn+iy)=a*scj10; % Calculate Denominator and Numerator of Reflex Coefficient (100) Rup=Rup+dbi10*tbi10';Rdown=Rdown+tbi10*tbi10'; Dbsqr=Dbsqr+dbi10*dbi10'; clear dbi10 tbi10 xc10nrm sb10 db10 sc10 sc10y scj10 tanv tan2vtan2vp a ix iy jx jy i j temp ; clear dy1 dy2 sc10y1 sc10y2 dymax1 dymax2; % Reset to Initial Value (*scf) xbnrm=xbn;ynrm=yn; % Read Image SB(120) and SC(120) =120 v=120;disp(v) fid=fopen(' matlab reflex T t120f1.img','r'); temp1=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t120f2.img','r'); temp2=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t120f3.img','r'); temp3=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t120f4.img','r'); temp4=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t120f5.img','r'); temp5=fread(fid,[ymax,xmax],'short');status=fclose(fid); temp=(temp1+temp2+temp3+temp4+temp5)*0.2; clear temp1 temp2 temp3 temp4 temp5;temp=temp-drk; at=temp(1:ymax/2+20,2:xmax);bt=temp(ymax/2+21:ymax,1:xmax-1); temp=[bt; at];clear at;clear bt; sb10=temp(sby1:sby2,sbx1:sbx2);db10=sb10-sb18;db10=medfilt2(db10); % Resize SC(120) in Y-axis C Y dymax1=5;dymax2=11;sc10y1=sby1-dymax1;sc10y2=sby2+dymax2; djx=-16; sc10x1=round(sbx1-cos((180-v)*pi/180)*(sbx2-sbx1))+djx;sc10x2=sbx1-3; ynrm=sby2-sby1+1;xc10nrm=sc10x2-sc10x1+1;

18 17 i=0;sc10y=zeros(ynrm,xc10nrm); for ix=sc10x1:sc10x2; dy1=round((sc10x2-ix)*dymax1/(xc10nrm-1)); dy2=round((sc10x2-ix)*dymax2/(xc10nrm-1)); i=i+1;sc10y(1:ynrm,i)=imresize(temp(sby1-dy1:sby2+dy2-1,ix),[ynrm,1]); sc10=sc10y;tic;clear temp; % Compress Image for Faster Calculation{ } xbnrm=round(xbnrm/scf);ynrm=round(ynrm/scf);xc10nrm=round(xc10nrm/scf); db10=imresize(db10,[ynrm xbnrm]);sc10=imresize(sc10,[ynrm xc10nrm]); % Calculate Matrix A(120),DBi(120),and TBi(120) { A D Bi T Bi } dbi10=reshape(db10,[1 xbnrm*ynrm]);scj10=reshape(sc10,[xc10nrm*ynrm 1]); tan2v=(tan(v*pi/180))^2;tanv=tan(v*pi/180);tan2vp=tan2v/pi; a=zeros(1,ynrm*xc10nrm); tbi10=zeros(1,ynrm*xbnrm); sbrate=(xc10nrm*scf+djx)/(xc10nrm*scf) for ix=1:xbnrm; ixix5t2p=(ix-0.5)*tan2vp; ixyn=(ix-1)*ynrm; for iy=1:ynrm; for jx=1:xc10nrm; jxm=jx*sbrate; ijx2jix2=(ix+jxm-1)^2+(jxm-0.5)^2*tan2v; bunsi=ixix5t2p*(jxm-0.5); jxyn=(jx-1)*ynrm; for jy=1:ynrm; a(jxyn+jy)=bunsi/(ijx2jix2+(iy-jy)^2)^2; tbi10(ixyn+iy)=a*scj10; % Calculate Denominator and Numerator of Reflex Coefficient (100,120) Rup=Rup+dbi10*tbi10';Rdown=Rdown+tbi10*tbi10'; Dbsqr=Dbsqr+dbi10*dbi10'; clear dbi10 tbi10 xc10nrm sb10 db10 sc10 sc10y scj10 tanv tan2vtan2vp a ix iy jx jy i j temp ; clear dy1 dy2 sc10y1 sc10y2 dymax1 dymax2; % Reset to Initial Value (*scf) xbnrm=xbn;ynrm=yn; % Read Image SB(140) and SC(140) =140 v=140;disp(v) fid=fopen(' matlab reflex T t140f1.img','r'); temp1=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t140f2.img','r'); temp2=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t140f3.img','r');

19 temp3=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t140f4.img','r'); temp4=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t140f5.img','r'); temp5=fread(fid,[ymax,xmax],'short');status=fclose(fid); temp=(temp1+temp2+temp3+temp4+temp5)*0.2; clear temp1 temp2 temp3 temp4 temp5;temp=temp-drk; at=temp(1:ymax/2+20,2:xmax);bt=temp(ymax/2+21:ymax,1:xmax-1); temp=[bt; at];clear at;clear bt; sb10=temp(sby1:sby2,sbx1:sbx2);db10=sb10-sb18;db10=medfilt2(db10); % Resize SC(140) in Y-axis C Y dymax1=4;dymax2=8;sc10y1=sby1-dymax1;sc10y2=sby2+dymax2; djx=-14; sc10x1=round(sbx1-cos((180-v)*pi/180)*(sbx2-sbx1))+djx;sc10x2=sbx1-3; ynrm=sby2-sby1+1;xc10nrm=sc10x2-sc10x1+1; i=0;sc10y=zeros(ynrm,xc10nrm); for ix=sc10x1:sc10x2; dy1=round((sc10x2-ix)*dymax1/(xc10nrm-1)); dy2=round((sc10x2-ix)*dymax2/(xc10nrm-1)); i=i+1;sc10y(1:ynrm,i)=imresize(temp(sby1-dy1:sby2+dy2-1,ix),[ynrm,1]); sc10=sc10y;tic;clear temp; % Compress Image for Faster Calculation{ } xbnrm=round(xbnrm/scf);ynrm=round(ynrm/scf);xc10nrm=round(xc10nrm/scf); db10=imresize(db10,[ynrm xbnrm]);sc10=imresize(sc10,[ynrm xc10nrm]); % Calculate Matrix A(140),DBi(140),and TBi(140) { A D Bi T Bi } dbi10=reshape(db10,[1 xbnrm*ynrm]);scj10=reshape(sc10,[xc10nrm*ynrm 1]); tan2v=(tan(v*pi/180))^2;tanv=tan(v*pi/180);tan2vp=tan2v/pi; a=zeros(1,ynrm*xc10nrm); tbi10=zeros(1,ynrm*xbnrm); sbrate=(xc10nrm*scf+djx)/(xc10nrm*scf) for ix=1:xbnrm; ixix5t2p=(ix-0.5)*tan2vp; ixyn=(ix-1)*ynrm; for iy=1:ynrm; for jx=1:xc10nrm; jxm=jx*sbrate; ijx2jix2=(ix+jxm-1)^2+(jxm-0.5)^2*tan2v; bunsi=ixix5t2p*(jxm-0.5); jxyn=(jx-1)*ynrm; for jy=1:ynrm; a(jxyn+jy)=bunsi/(ijx2jix2+(iy-jy)^2)^2; tbi10(ixyn+iy)=a*scj10;

20 19 % Calculate Denominator and Numerator of Reflex Coefficient (100,120,140) Rup=Rup+dbi10*tbi10';Rdown=Rdown+tbi10*tbi10'; Dbsqr=Dbsqr+dbi10*dbi10'; clear dbi10 tbi10 xc10nrm sb10 db10 sc10 sc10y scj10 tanv tan2vtan2vp a ix iy jx jy i j temp ; clear dy1 dy2 sc10y1 sc10y2 dymax1 dymax2; % Reset to Initial Value (*scf) xbnrm=xbn;ynrm=yn; % Read Image SB(160) and SC(160) =160 v=160;disp(v) fid=fopen(' matlab reflex T t160f1.img','r'); temp1=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t160f2.img','r'); temp2=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t160f3.img','r'); temp3=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t160f4.img','r'); temp4=fread(fid,[ymax,xmax],'short');status=fclose(fid); fid=fopen(' matlab reflex T t160f5.img','r'); temp5=fread(fid,[ymax,xmax],'short');status=fclose(fid); temp=(temp1+temp2+temp3+temp4+temp5)*0.2; clear temp1 temp2 temp3 temp4 temp5; temp=temp-drk; at=temp(1:ymax/2+20,2:xmax);bt=temp(ymax/2+21:ymax,1:xmax-1); temp=[bt; at];clear at;clear bt; sb10=temp(sby1:sby2,sbx1:sbx2); db10=sb10-sb18;db10=medfilt2(db10); % Resize SC(160) in Y-axis C Y dymax1=2;dymax2=4;sc10y1=sby1-dymax1;sc10y2=sby2+dymax2; djx=-9; sc10x1=round(sbx1-cos((180-v)*pi/180)*(sbx2-sbx1))+djx;sc10x2=sbx1-3; ynrm=sby2-sby1+1;xc10nrm=sc10x2-sc10x1+1; i=0;sc10y=zeros(ynrm,xc10nrm); for ix=sc10x1:sc10x2; dy1=round((sc10x2-ix)*dymax1/(xc10nrm-1)); dy2=round((sc10x2-ix)*dymax2/(xc10nrm-1)); i=i+1;sc10y(1:ynrm,i)=imresize(temp(sby1-dy1:sby2+dy2-1,ix),[ynrm,1]); sc10=sc10y;tic;clear temp; % Compress Image for Faster Calculation { } xbnrm=round(xbnrm/scf);ynrm=round(ynrm/scf);xc10nrm=round(xc10nrm/scf); db10=imresize(db10,[ynrm xbnrm]);sc10=imresize(sc10,[ynrm xc10nrm]);

21 % Calculate Matrix A(160),DBi(160),and TBi(160) { A D Bi T Bi } dbi10=reshape(db10,[1 xbnrm*ynrm]);scj10=reshape(sc10,[xc10nrm*ynrm 1]); tan2v=(tan(v*pi/180))^2;tanv=tan(v*pi/180);tan2vp=tan2v/pi; a=zeros(1,ynrm*xc10nrm); tbi10=zeros(1,ynrm*xbnrm); sbrate=(xc10nrm*scf+djx)/(xc10nrm*scf) for ix=1:xbnrm; ixix5t2p=(ix-0.5)*tan2vp; ixyn=(ix-1)*ynrm; for iy=1:ynrm; for jx=1:xc10nrm; jxm=jx*sbrate; ijx2jix2=(ix+jxm-1)^2+(jxm-0.5)^2*tan2v; bunsi=ixix5t2p*(jxm-0.5); jxyn=(jx-1)*ynrm; for jy=1:ynrm; a(jxyn+jy)=bunsi/(ijx2jix2+(iy-jy)^2)^2; tbi10(ixyn+iy)=a*scj10; % Calculate Denominator and Numerator of Reflex Coefficient (100,120,140,160) Rup=Rup+dbi10*tbi10';Rdown=Rdown+tbi10*tbi10'; Dbsqr=Dbsqr+dbi10*dbi10'; clear dbi10 tbi10 xc10nrm sb10 db10 sc10 sc10y scj10 tanv tan2vtan2vp a ix iy jx jy i j temp ; clear dy1 dy2 sc10y1 sc10y2 dymax1 dymax2;clear drk; % Calculate Reflex Coefficient { } R=Rup/Rdown;display(R) toc % Calculate Uncertainty of R { } qq=dbsqr-2*r*rup+r^2*rdown; urred=(qq/rdown)^0.5;display(urred)

Gmech08.dvi

Gmech08.dvi 145 13 13.1 13.1.1 0 m mg S 13.1 F 13.1 F /m S F F 13.1 F mg S F F mg 13.1: m d2 r 2 = F + F = 0 (13.1) 146 13 F = F (13.2) S S S S S P r S P r r = r 0 + r (13.3) r 0 S S m d2 r 2 = F (13.4) (13.3) d 2

More information

28 Horizontal angle correction using straight line detection in an equirectangular image

28 Horizontal angle correction using straight line detection in an equirectangular image 28 Horizontal angle correction using straight line detection in an equirectangular image 1170283 2017 3 1 2 i Abstract Horizontal angle correction using straight line detection in an equirectangular image

More information

(a) 1 (b) 3. Gilbert Pernicka[2] Treibitz Schechner[3] Narasimhan [4] Kim [5] Nayar [6] [7][8][9] 2. X X X [10] [11] L L t L s L = L t + L s

(a) 1 (b) 3. Gilbert Pernicka[2] Treibitz Schechner[3] Narasimhan [4] Kim [5] Nayar [6] [7][8][9] 2. X X X [10] [11] L L t L s L = L t + L s 1 1 1, Extraction of Transmitted Light using Parallel High-frequency Illumination Kenichiro Tanaka 1 Yasuhiro Mukaigawa 1 Yasushi Yagi 1 Abstract: We propose a new sharpening method of transmitted scene

More information

SP qxd

SP qxd JAXA-SP-06-020 JAXA/WINTE PSP Research and Development of Low-Speed PSP Measurement System in Wind Tunnel Technology enter of JAXA Kazunori MITSUO, Mitsuru KURITA, Kazuyuki NAKAKITA Keisuke FUJII and Shigeya

More information

JFE.dvi

JFE.dvi ,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho

More information

SO(2)

SO(2) TOP URL http://amonphys.web.fc2.com/ 1 12 3 12.1.................................. 3 12.2.......................... 4 12.3............................. 5 12.4 SO(2).................................. 6

More information

[2] OCR [3], [4] [5] [6] [4], [7] [8], [9] 1 [10] Fig. 1 Current arrangement and size of ruby. 2 Fig. 2 Typography combined with printing

[2] OCR [3], [4] [5] [6] [4], [7] [8], [9] 1 [10] Fig. 1 Current arrangement and size of ruby. 2 Fig. 2 Typography combined with printing 1,a) 1,b) 1,c) 2012 11 8 2012 12 18, 2013 1 27 WEB Ruby Removal Filters Using Genetic Programming for Early-modern Japanese Printed Books Taeka Awazu 1,a) Masami Takata 1,b) Kazuki Joe 1,c) Received: November

More information

xx/xx Vol. Jxx A No. xx 1 Fig. 1 PAL(Panoramic Annular Lens) PAL(Panoramic Annular Lens) PAL (2) PAL PAL 2 PAL 3 2 PAL 1 PAL 3 PAL PAL 2. 1 PAL

xx/xx Vol. Jxx A No. xx 1 Fig. 1 PAL(Panoramic Annular Lens) PAL(Panoramic Annular Lens) PAL (2) PAL PAL 2 PAL 3 2 PAL 1 PAL 3 PAL PAL 2. 1 PAL PAL On the Precision of 3D Measurement by Stereo PAL Images Hiroyuki HASE,HirofumiKAWAI,FrankEKPAR, Masaaki YONEDA,andJien KATO PAL 3 PAL Panoramic Annular Lens 1985 Greguss PAL 1 PAL PAL 2 3 2 PAL DP

More information

mt_4.dvi

mt_4.dvi ( ) 2006 1 PI 1 1 1.1................................. 1 1.2................................... 1 2 2 2.1...................................... 2 2.1.1.......................... 2 2.1.2..............................

More information

2000年度『数学展望 I』講義録

2000年度『数学展望 I』講義録 2000 I I IV I II 2000 I I IV I-IV. i ii 3.10 (http://www.math.nagoya-u.ac.jp/ kanai/) 2000 A....1 B....4 C....10 D....13 E....17 Brouwer A....21 B....26 C....33 D....39 E. Sperner...45 F....48 A....53

More information

chap1.dvi

chap1.dvi 1 1 007 1 e iθ = cos θ + isin θ 1) θ = π e iπ + 1 = 0 1 ) 3 11 f 0 r 1 1 ) k f k = 1 + r) k f 0 f k k = 01) f k+1 = 1 + r)f k ) f k+1 f k = rf k 3) 1 ) ) ) 1+r/)f 0 1 1 + r/) f 0 = 1 + r + r /4)f 0 1 f

More information

TM

TM NALTR-1390 TR-1390 ISSN 0452-2982 UDC 533.6.013.1 533.6.013.4 533.6.69.048 NAL TECHNICAL REPORT OF NATIONAL AEROSPACE LABORATORY TR-1390 e N 1999 11 NATIONAL AEROSPACE LABORATORY ... 1 e N... 2 Orr-Sommerfeld...

More information

撮 影

撮 影 DC cathode ray tube, 2.2 log log log + log log / / / A method determining tone conversion characteristics of digital still camera from two pictorial images Tone conversion characteristic Luminance

More information

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting

More information

kiyo5_1-masuzawa.indd

kiyo5_1-masuzawa.indd .pp. A Study on Wind Forecast using Self-Organizing Map FUJIMATSU Seiichiro, SUMI Yasuaki, UETA Takuya, KOBAYASHI Asuka, TSUKUTANI Takao, FUKUI Yutaka SOM SOM Elman SOM SOM Elman SOM Abstract : Now a small

More information

29 1 6 1 1 1.1 1.1 1.1( ) 1.1( ) 1.1: 2 1.2 1.2( ) 4 4 1 2,3,4 1 2 1 2 1.2: 1,2,3,4 a 1 2a 6 2 2,3,4 1,2,3,4 1.2( ) 4 1.2( ) 3 1.2( ) 1.3 1.3 1.3: 4 1.4 1.4 1.4: 1.5 1.5 1 2 1 a a R = l a l 5 R = l a +

More information

塗装深み感の要因解析

塗装深み感の要因解析 17 Analysis of Factors for Paint Depth Feeling Takashi Wada, Mikiko Kawasumi, Taka-aki Suzuki ( ) ( ) ( ) The appearance and quality of objects are controlled by paint coatings on the surfaces of the objects.

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc2.com/ 1 30 3 30.1.............. 3 30.2........................... 4 30.3...................... 5 30.4........................ 6 30.5.................................. 8 30.6...............................

More information

(MIRU2008) HOG Histograms of Oriented Gradients (HOG)

(MIRU2008) HOG Histograms of Oriented Gradients (HOG) (MIRU2008) 2008 7 HOG - - E-mail: katsu0920@me.cs.scitec.kobe-u.ac.jp, {takigu,ariki}@kobe-u.ac.jp Histograms of Oriented Gradients (HOG) HOG Shape Contexts HOG 5.5 Histograms of Oriented Gradients D Human

More information

Classic HD:ŠŸŠp”Ò:Discovery:‘‚ŠÞ:‚²„¤:Ÿ_Ł¶:Simulation_for_HRO.dvi

Classic HD:ŠŸŠp”Ò:Discovery:‘‚ŠÞ:‚²„¤:Ÿ_Ł¶:Simulation_for_HRO.dvi HRO 52 2002 7 19 HRO (Ham Radio Observation) HRO HRO 1 2 100km Microsoft Excel HRO 1999 12 12 16 HRO HRO 1 2 1.1........................................... 2 1.2 HRO............................ 2 1.3 HRO...................................

More information

Gmech08.dvi

Gmech08.dvi 51 5 5.1 5.1.1 P r P z θ P P P z e r e, z ) r, θ, ) 5.1 z r e θ,, z r, θ, = r sin θ cos = r sin θ sin 5.1) e θ e z = r cos θ r, θ, 5.1: 0 r

More information

JAB RL504:201 JAB NOTE 4 2 January 25, A B

JAB RL504:201 JAB NOTE 4 2 January 25, A B JAB RL504:201 JAB NOTE 4 JAB RL504:201 2201125 1200025 200-0-25-1/- 2 201-01-25 JAB RL504:201 JAB NOTE 4 2 January 25, 201............... 6 5.1... 6 5.2... 6 5. A... 7 5.4 B... 7 5.5... 8 5.6... 9 5.7...

More information

Undulator.dvi

Undulator.dvi X X 1 1 2 Free Electron Laser: FEL 2.1 2 2 3 SACLA 4 SACLA [1]-[6] [7] 1: S N λ [9] XFEL OHO 13 X [8] 2 2.1 2(a) (c) z y y (a) S N 90 λ u 4 [10, 11] Halbach (b) 2: (a) (b) (c) (c) 1 2 [11] B y = n=1 B

More information

(iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y = 0., y x, y = x. (v) 1x = x. (vii) (α + β)x = αx + βx. (viii) (αβ)x = α(βx)., V, C.,,., (1)

(iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y = 0., y x, y = x. (v) 1x = x. (vii) (α + β)x = αx + βx. (viii) (αβ)x = α(βx)., V, C.,,., (1) 1. 1.1...,. 1.1.1 V, V x, y, x y x + y x + y V,, V x α, αx αx V,, (i) (viii) : x, y, z V, α, β C, (i) x + y = y + x. (ii) (x + y) + z = x + (y + z). 1 (iii) 0 V, x V, x + 0 = x. 0. (iv) x V, y V, x + y

More information

zz + 3i(z z) + 5 = 0 + i z + i = z 2i z z z y zz + 3i (z z) + 5 = 0 (z 3i) (z + 3i) = 9 5 = 4 z 3i = 2 (3i) zz i (z z) + 1 = a 2 {

zz + 3i(z z) + 5 = 0 + i z + i = z 2i z z z y zz + 3i (z z) + 5 = 0 (z 3i) (z + 3i) = 9 5 = 4 z 3i = 2 (3i) zz i (z z) + 1 = a 2 { 04 zz + iz z) + 5 = 0 + i z + i = z i z z z 970 0 y zz + i z z) + 5 = 0 z i) z + i) = 9 5 = 4 z i = i) zz i z z) + = a {zz + i z z) + 4} a ) zz + a + ) z z) + 4a = 0 4a a = 5 a = x i) i) : c Darumafactory

More information

December 28, 2018

December 28, 2018 e-mail : kigami@i.kyoto-u.ac.jp December 28, 28 Contents 2............................. 3.2......................... 7.3..................... 9.4................ 4.5............. 2.6.... 22 2 36 2..........................

More information

ライトカーブ観測から何がわかるか

ライトカーブ観測から何がわかるか 2004.7.2 M 2.5log F F 10 0.4*M M F S= (abc)[sin 2 A(sin 2 ( )/a 2 +cos 2 ( )/b 2 )+cos 2 (A)/c 2 ] 1/2 a,b,c A aspect angle S M = (abc)[sin 2 (A)/b 2 +cos 2 (A)/c 2 ] 1/2 S m = (abc)[sin 2 (A)/a 2 +cos

More information

29

29 9 .,,, 3 () C k k C k C + C + C + + C 8 + C 9 + C k C + C + C + C 3 + C 4 + C 5 + + 45 + + + 5 + + 9 + 4 + 4 + 5 4 C k k k ( + ) 4 C k k ( k) 3 n( ) n n n ( ) n ( ) n 3 ( ) 3 3 3 n 4 ( ) 4 4 4 ( ) n n

More information

1 Kinect for Windows M = [X Y Z] T M = [X Y Z ] T f (u,v) w 3.2 [11] [7] u = f X +u Z 0 δ u (X,Y,Z ) (5) v = f Y Z +v 0 δ v (X,Y,Z ) (6) w = Z +

1 Kinect for Windows M = [X Y Z] T M = [X Y Z ] T f (u,v) w 3.2 [11] [7] u = f X +u Z 0 δ u (X,Y,Z ) (5) v = f Y Z +v 0 δ v (X,Y,Z ) (6) w = Z + 3 3D 1,a) 1 1 Kinect (X, Y) 3D 3D 1. 2010 Microsoft Kinect for Windows SDK( (Kinect) SDK ) 3D [1], [2] [3] [4] [5] [10] 30fps [10] 3 Kinect 3 Kinect Kinect for Windows SDK 3 Microsoft 3 Kinect for Windows

More information

IHI Robust Path Planning against Position Error for UGVs in Rough Terrain Yuki DOI, Yonghoon JI, Yusuke TAMURA(University of Tokyo), Yuki IKEDA, Atsus

IHI Robust Path Planning against Position Error for UGVs in Rough Terrain Yuki DOI, Yonghoon JI, Yusuke TAMURA(University of Tokyo), Yuki IKEDA, Atsus IHI Robust Path Planning against Position Error for UGVs in Rough Terrain Yuki DOI, Yonghoon JI, Yusuke TAMURA(University of Tokyo), Yuki IKEDA, Atsushi UMEMURA, Yoshiharu KANESHIMA, Hiroki MURAKAMI(IHI

More information

1 [1, 2, 3, 4, 5, 8, 9, 10, 12, 15] The Boston Public Schools system, BPS (Deferred Acceptance system, DA) (Top Trading Cycles system, TTC) cf. [13] [

1 [1, 2, 3, 4, 5, 8, 9, 10, 12, 15] The Boston Public Schools system, BPS (Deferred Acceptance system, DA) (Top Trading Cycles system, TTC) cf. [13] [ Vol.2, No.x, April 2015, pp.xx-xx ISSN xxxx-xxxx 2015 4 30 2015 5 25 253-8550 1100 Tel 0467-53-2111( ) Fax 0467-54-3734 http://www.bunkyo.ac.jp/faculty/business/ 1 [1, 2, 3, 4, 5, 8, 9, 10, 12, 15] The

More information

AHPを用いた大相撲の新しい番付編成

AHPを用いた大相撲の新しい番付編成 5304050 2008/2/15 1 2008/2/15 2 42 2008/2/15 3 2008/2/15 4 195 2008/2/15 5 2008/2/15 6 i j ij >1 ij ij1/>1 i j i 1 ji 1/ j ij 2008/2/15 7 1 =2.01/=0.5 =1.51/=0.67 2008/2/15 8 1 2008/2/15 9 () u ) i i i

More information

Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Enginee

Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Enginee Study on Application of the cos a Method to Neutron Stress Measurement Toshihiko SASAKI*3 and Yukio HIROSE Department of Materials Science and Engineering, Kanazawa University, Kakuma-machi, Kanazawa-shi,

More information

1. 4cm 16 cm 4cm 20cm 18 cm L λ(x)=ax [kg/m] A x 4cm A 4cm 12 cm h h Y 0 a G 0.38h a b x r(x) x y = 1 h 0.38h G b h X x r(x) 1 S(x) = πr(x) 2 a,b, h,π

1. 4cm 16 cm 4cm 20cm 18 cm L λ(x)=ax [kg/m] A x 4cm A 4cm 12 cm h h Y 0 a G 0.38h a b x r(x) x y = 1 h 0.38h G b h X x r(x) 1 S(x) = πr(x) 2 a,b, h,π . 4cm 6 cm 4cm cm 8 cm λ()=a [kg/m] A 4cm A 4cm cm h h Y a G.38h a b () y = h.38h G b h X () S() = π() a,b, h,π V = ρ M = ρv G = M h S() 3 d a,b, h 4 G = 5 h a b a b = 6 ω() s v m θ() m v () θ() ω() dθ()

More information

AC Modeling and Control of AC Motors Seiji Kondo, Member 1. q q (1) PM (a) N d q Dept. of E&E, Nagaoka Unive

AC Modeling and Control of AC Motors Seiji Kondo, Member 1. q q (1) PM (a) N d q Dept. of E&E, Nagaoka Unive AC Moeling an Control of AC Motors Seiji Kono, Member 1. (1) PM 33 54 64. 1 11 1(a) N 94 188 163 1 Dept. of E&E, Nagaoka University of Technology 163 1, Kamitomioka-cho, Nagaoka, Niigata 94 188 (a) 巻数

More information

x x x 2, A 4 2 Ax.4 A A A A λ λ 4 λ 2 A λe λ λ2 5λ + 6 0,...λ 2, λ 2 3 E 0 E 0 p p Ap λp λ 2 p 4 2 p p 2 p { 4p 2 2p p + 2 p, p 2 λ {

x x x 2, A 4 2 Ax.4 A A A A λ λ 4 λ 2 A λe λ λ2 5λ + 6 0,...λ 2, λ 2 3 E 0 E 0 p p Ap λp λ 2 p 4 2 p p 2 p { 4p 2 2p p + 2 p, p 2 λ { K E N Z OU 2008 8. 4x 2x 2 2 2 x + x 2. x 2 2x 2, 2 2 d 2 x 2 2.2 2 3x 2... d 2 x 2 5 + 6x 0 2 2 d 2 x 2 + P t + P 2tx Qx x x, x 2 2 2 x 2 P 2 tx P tx 2 + Qx x, x 2. d x 4 2 x 2 x x 2.3 x x x 2, A 4 2

More information

I

I I 6 4 10 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

5 c P 5 kn n t π (.5 P 7 MP π (.5 n t n cos π. MP 6 4 t sin π 6 cos π 6.7 MP 4 P P N i i i i N i j F j ii N i i ii F j i i N ii li i F j i ij li i i i

5 c P 5 kn n t π (.5 P 7 MP π (.5 n t n cos π. MP 6 4 t sin π 6 cos π 6.7 MP 4 P P N i i i i N i j F j ii N i i ii F j i i N ii li i F j i ij li i i i i j ij i j ii,, i j ij ij ij (, P P P P θ N θ P P cosθ N F N P cosθ F Psinθ P P F P P θ N P cos θ cos θ cosθ F P sinθ cosθ sinθ cosθ sinθ 5 c P 5 kn n t π (.5 P 7 MP π (.5 n t n cos π. MP 6 4 t sin π 6

More information

1).1-5) - 9 -

1).1-5) - 9 - - 8 - 1).1-5) - 9 - ε = ε xx 0 0 0 ε xx 0 0 0 ε xx (.1 ) z z 1 z ε = ε xx ε x y 0 - ε x y ε xx 0 0 0 ε zz (. ) 3 xy ) ε xx, ε zz» ε x y (.3 ) ε ij = ε ij ^ (.4 ) 6) xx, xy ε xx = ε xx + i ε xx ε xy = ε

More information

Abstract This paper concerns with a method of dynamic image cognition. Our image cognition method has two distinguished features. One is that the imag

Abstract This paper concerns with a method of dynamic image cognition. Our image cognition method has two distinguished features. One is that the imag 2004 RGB A STUDY OF RGB COLOR INFORMATION AND ITS APPLICATION 03R3237 Abstract This paper concerns with a method of dynamic image cognition. Our image cognition method has two distinguished features. One

More information

arctan 1 arctan arctan arctan π = = ( ) π = 4 = π = π = π = =

arctan 1 arctan arctan arctan π = = ( ) π = 4 = π = π = π = = arctan arctan arctan arctan 2 2000 π = 3 + 8 = 3.25 ( ) 2 8 650 π = 4 = 3.6049 9 550 π = 3 3 30 π = 3.622 264 π = 3.459 3 + 0 7 = 3.4085 < π < 3 + 7 = 3.4286 380 π = 3 + 77 250 = 3.46 5 3.45926 < π < 3.45927

More information

Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue

Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Date 2017 URL http://hdl.handle.net/2433/229150 Right

More information

SO(3) 49 u = Ru (6.9), i u iv i = i u iv i (C ) π π : G Hom(V, V ) : g D(g). π : R 3 V : i 1. : u u = u 1 u 2 u 3 (6.10) 6.2 i R α (1) = 0 cos α

SO(3) 49 u = Ru (6.9), i u iv i = i u iv i (C ) π π : G Hom(V, V ) : g D(g). π : R 3 V : i 1. : u u = u 1 u 2 u 3 (6.10) 6.2 i R α (1) = 0 cos α SO(3) 48 6 SO(3) t 6.1 u, v u = u 1 1 + u 2 2 + u 3 3 = u 1 e 1 + u 2 e 2 + u 3 e 3, v = v 1 1 + v 2 2 + v 3 3 = v 1 e 1 + v 2 e 2 + v 3 e 3 (6.1) i (e i ) e i e j = i j = δ ij (6.2) ( u, v ) = u v = ij

More information

<4D F736F F D DB82CC88F892A38BAD937893C190AB76355F8D5A897B8CE3325F2E646F63>

<4D F736F F D DB82CC88F892A38BAD937893C190AB76355F8D5A897B8CE3325F2E646F63> Asymmetry of Seismic Displacement Response of Highway Bridge Supported by Spread Foundation Takeshi SHIMABUKURO* Rei FUJITA** and Norihiko YAMASHITA*** The objective of this paper is to discuss the causes

More information

日立金属技報 Vol.34

日立金属技報 Vol.34 Influence of Misorientation Angle between Adjacent Grains on Magnetization Reversal in Nd-Fe-B Sintered Magnet Tomohito Maki Rintaro Ishii Mitsutoshi Natsumeda Takeshi Nishiuchi Ryo Uchikoshi Masaaki Takezawa

More information

149 (Newell [5]) Newell [5], [1], [1], [11] Li,Ryu, and Song [2], [11] Li,Ryu, and Song [2], [1] 1) 2) ( ) ( ) 3) T : 2 a : 3 a 1 :

149 (Newell [5]) Newell [5], [1], [1], [11] Li,Ryu, and Song [2], [11] Li,Ryu, and Song [2], [1] 1) 2) ( ) ( ) 3) T : 2 a : 3 a 1 : Transactions of the Operations Research Society of Japan Vol. 58, 215, pp. 148 165 c ( 215 1 2 ; 215 9 3 ) 1) 2) :,,,,, 1. [9] 3 12 Darroch,Newell, and Morris [1] Mcneil [3] Miller [4] Newell [5, 6], [1]

More information

Input image Initialize variables Loop for period of oscillation Update height map Make shade image Change property of image Output image Change time L

Input image Initialize variables Loop for period of oscillation Update height map Make shade image Change property of image Output image Change time L 1,a) 1,b) 1/f β Generation Method of Animation from Pictures with Natural Flicker Abstract: Some methods to create animation automatically from one picture have been proposed. There is a method that gives

More information

空力騒音シミュレータの開発

空力騒音シミュレータの開発 41 COSMOS-V, an Aerodynamic Noise Simulator Nariaki Horinouchi COSMOS-V COSMOS-V COSMOS-V 3 The present and future computational problems of the aerodynamic noise analysis using COSMOS-V, our in-house

More information

udc-2.dvi

udc-2.dvi 13 0.5 2 0.5 2 1 15 2001 16 2009 12 18 14 No.39, 2010 8 2009b 2009a Web Web Q&A 2006 2007a20082009 2007b200720082009 20072008 2009 2009 15 1 2 2 2.1 18 21 1 4 2 3 1(a) 1(b) 1(c) 1(d) 1) 18 16 17 21 10

More information

Part () () Γ Part ,

Part () () Γ Part , Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35

More information

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =,

.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =, [ ] IC. r, θ r, θ π, y y = 3 3 = r cos θ r sin θ D D = {, y ; y }, y D r, θ ep y yddy D D 9 s96. d y dt + 3dy + y = cos t dt t = y = e π + e π +. t = π y =.9 s6.3 d y d + dy d + y = y =, dy d = 3 a, b

More information

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi) 0. A A = 4 IC () det A () A () x + y + z = x y z X Y Z = A x y z ( 5) ( s5590) 0. a + b + c b c () a a + b + c c a b a + b + c 0 a b c () a 0 c b b c 0 a c b a 0 0. A A = 7 5 4 5 0 ( 5) ( s5590) () A ()

More information

note1.dvi

note1.dvi (1) 1996 11 7 1 (1) 1. 1 dx dy d x τ xx x x, stress x + dx x τ xx x+dx dyd x x τ xx x dyd y τ xx x τ xx x+dx d dx y x dy 1. dx dy d x τ xy x τ x ρdxdyd x dx dy d ρdxdyd u x t = τ xx x+dx dyd τ xx x dyd

More information

IPSJ SIG Technical Report Vol.2009-BIO-17 No /5/26 DNA 1 1 DNA DNA DNA DNA Correcting read errors on DNA sequences determined by Pyrosequencing

IPSJ SIG Technical Report Vol.2009-BIO-17 No /5/26 DNA 1 1 DNA DNA DNA DNA Correcting read errors on DNA sequences determined by Pyrosequencing DNA 1 1 DNA DNA DNA DNA Correcting read errors on DNA sequences determined by Pyrosequencing Youhei Namiki 1 and Yutaka Akiyama 1 Pyrosequencing, one of the DNA sequencing technologies, allows us to determine

More information

‚æ4“ƒ.ren

‚æ4“ƒ.ren 69 1 1 13 14 70 2 3 1972 4 5 1992 6 7 1980 100 1997 71 226 2100 8 100 50 9 21 21 21 21 10 11 2 1968 12 13 72 1980 14 15 16 17 73 18 18 20 5 2002 12 9 19 1964 1980 20 21 22 74 23 100 10 10 100 101 1 101

More information

数学の基礎訓練I

数学の基礎訓練I I 9 6 13 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 3 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

+ 1 ( ) I IA i i i 1 n m a 11 a 1j a 1m A = a i1 a ij a im a n1 a nj a nm.....

+   1 ( ) I IA i i i 1 n m a 11 a 1j a 1m A = a i1 a ij a im a n1 a nj a nm..... + http://krishnathphysaitama-uacjp/joe/matrix/matrixpdf 1 ( ) I IA i i i 1 n m a 11 a 1j a 1m A = a i1 a ij a im a n1 a nj a nm (1) n m () (n, m) ( ) n m B = ( ) 3 2 4 1 (2) 2 2 ( ) (2, 2) ( ) C = ( 46

More information

_TZ_4797-haus-local

_TZ_4797-haus-local 1.1.................................... 3.3.................................. 4.4......................... 8.5... 10.6.................... 1.7... 14 3 16 3.1 ()........................... 16 3. 7... 17

More information

Microsoft Word - 信号処理3.doc

Microsoft Word - 信号処理3.doc Junji OHTSUBO 2012 FFT FFT SN sin cos x v ψ(x,t) = f (x vt) (1.1) t=0 (1.1) ψ(x,t) = A 0 cos{k(x vt) + φ} = A 0 cos(kx ωt + φ) (1.2) A 0 v=ω/k φ ω k 1.3 (1.2) (1.2) (1.2) (1.1) 1.1 c c = a + ib, a = Re[c],

More information

01-加藤 実-5.02

01-加藤 実-5.02 Bull. Natl. Mus. Nat. Sci., Ser. E, 30, pp. 1 13, December 21, 2007 1 2 3 1 169 0073 3 23 1 2 523 0058 961 3 248 0036 3 5 6 The Mechanism of the Automatic Wari-koma Dial in the Japanese Clocks and its

More information

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 第 2 版 1 刷発行時のものです.

医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます.   このサンプルページの内容は, 第 2 版 1 刷発行時のものです. 医系の統計入門第 2 版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009192 このサンプルページの内容は, 第 2 版 1 刷発行時のものです. i 2 t 1. 2. 3 2 3. 6 4. 7 5. n 2 ν 6. 2 7. 2003 ii 2 2013 10 iii 1987

More information

弾性定数の対称性について

弾性定数の対称性について () by T. oyama () ij C ij = () () C, C, C () ij ji ij ijlk ij ij () C C C C C C * C C C C C * * C C C C = * * * C C C * * * * C C * * * * * C () * P (,, ) P (,, ) lij = () P (,, ) P(,, ) (,, ) P (, 00,

More information

18 2 20 W/C W/C W/C 4-4-1 0.05 1.0 1000 1. 1 1.1 1 1.2 3 2. 4 2.1 4 (1) 4 (2) 4 2.2 5 (1) 5 (2) 5 2.3 7 3. 8 3.1 8 3.2 ( ) 11 3.3 11 (1) 12 (2) 12 4. 14 4.1 14 4.2 14 (1) 15 (2) 16 (3) 17 4.3 17 5. 19

More information

Estimation of Photovoltaic Module Temperature Rise Motonobu Yukawa, Member, Masahisa Asaoka, Non-member (Mitsubishi Electric Corp.) Keigi Takahara, Me

Estimation of Photovoltaic Module Temperature Rise Motonobu Yukawa, Member, Masahisa Asaoka, Non-member (Mitsubishi Electric Corp.) Keigi Takahara, Me Estimation of Photovoltaic Module Temperature Rise Motonobu Yukawa, Member, Masahisa Asaoka, Non-member (Mitsubishi Electric Corp.) Keigi Takahara, Member (Okinawa Electric Power Co.,Inc.) Toshimitsu Ohshiro,

More information

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K

II 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K II. () 7 F 7 = { 0,, 2, 3, 4, 5, 6 }., F 7 a, b F 7, a b, F 7,. (a) a, b,,. (b) 7., 4 5 = 20 = 2 7 + 6, 4 5 = 6 F 7., F 7,., 0 a F 7, ab = F 7 b F 7. (2) 7, 6 F 6 = { 0,, 2, 3, 4, 5 },,., F 6., 0 0 a F

More information

1., 1 COOKPAD 2, Web.,,,,,,.,, [1]., 5.,, [2].,,.,.,, 5, [3].,,,.,, [4], 33,.,,.,,.. 2.,, 3.., 4., 5., ,. 1.,,., 2.,. 1,,

1., 1 COOKPAD 2, Web.,,,,,,.,, [1]., 5.,, [2].,,.,.,, 5, [3].,,,.,, [4], 33,.,,.,,.. 2.,, 3.., 4., 5., ,. 1.,,., 2.,. 1,, THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS TECHNICAL REPORT OF IEICE.,, 464 8601 470 0393 101 464 8601 E-mail: matsunagah@murase.m.is.nagoya-u.ac.jp, {ide,murase,hirayama}@is.nagoya-u.ac.jp,

More information

4 4 θ X θ P θ 4. 0, 405 P 0 X 405 X P 4. () 60 () 45 () 40 (4) 765 (5) 40 B 60 0 P = 90, = ( ) = X

4 4 θ X θ P θ 4. 0, 405 P 0 X 405 X P 4. () 60 () 45 () 40 (4) 765 (5) 40 B 60 0 P = 90, = ( ) = X 4 4. 4.. 5 5 0 A P P P X X X X +45 45 0 45 60 70 X 60 X 0 P P 4 4 θ X θ P θ 4. 0, 405 P 0 X 405 X P 4. () 60 () 45 () 40 (4) 765 (5) 40 B 60 0 P 0 0 + 60 = 90, 0 + 60 = 750 0 + 60 ( ) = 0 90 750 0 90 0

More information

3. ( 1 ) Linear Congruential Generator:LCG 6) (Mersenne Twister:MT ), L 1 ( 2 ) 4 4 G (i,j) < G > < G 2 > < G > 2 g (ij) i= L j= N

3. ( 1 ) Linear Congruential Generator:LCG 6) (Mersenne Twister:MT ), L 1 ( 2 ) 4 4 G (i,j) < G > < G 2 > < G > 2 g (ij) i= L j= N RMT 1 1 1 N L Q=L/N (RMT), RMT,,,., Box-Muller, 3.,. Testing Randomness by Means of RMT Formula Xin Yang, 1 Ryota Itoi 1 and Mieko Tanaka-Yamawaki 1 Random matrix theory derives, at the limit of both dimension

More information

空気の屈折率変調を光学的に検出する超指向性マイクロホン

空気の屈折率変調を光学的に検出する超指向性マイクロホン 23 2 1M36268 2 2 4 5 6 7 8 13 15 2 21 2 23 2 2 3 32 34 38 38 54 57 62 63 1-1 ( 1) ( 2) 1-1 a ( sinθ ) 2J D ( θ ) = 1 (1-1) kaka sinθ ( 3) 1-2 1 Back face hole Amplifier Diaphragm Equiphase wave surface

More information

206“ƒŁ\”ƒ-fl_“H„¤‰ZŁñ

206“ƒŁ\”ƒ-fl_“H„¤‰ZŁñ 51 206 51 63 2007 GIS 51 1 60 52 2 60 1 52 3 61 2 52 61 3 58 61 4 58 Summary 63 60 20022005 2004 40km 7,10025 2002 2005 19 3 19 GIS 2005GIS 2006 2002 2004 GIS 52 2062007 1 2004 GIS Fig.1 GIS ESRIArcView

More information

all.dvi

all.dvi 38 5 Cauchy.,,,,., σ.,, 3,,. 5.1 Cauchy (a) (b) (a) (b) 5.1: 5.1. Cauchy 39 F Q Newton F F F Q F Q 5.2: n n ds df n ( 5.1). df n n df(n) df n, t n. t n = df n (5.1) ds 40 5 Cauchy t l n mds df n 5.3: t

More information

20 4 20 i 1 1 1.1............................ 1 1.2............................ 4 2 11 2.1................... 11 2.2......................... 11 2.3....................... 19 3 25 3.1.............................

More information

1: A/B/C/D Fig. 1 Modeling Based on Difference in Agitation Method artisoc[7] A D 2017 Information Processing

1: A/B/C/D Fig. 1 Modeling Based on Difference in Agitation Method artisoc[7] A D 2017 Information Processing 1,a) 2,b) 3 Modeling of Agitation Method in Automatic Mahjong Table using Multi-Agent Simulation Hiroyasu Ide 1,a) Takashi Okuda 2,b) Abstract: Automatic mahjong table refers to mahjong table which automatically

More information

85 4

85 4 85 4 86 Copright c 005 Kumanekosha 4.1 ( ) ( t ) t, t 4.1.1 t Step! (Step 1) (, 0) (Step ) ±V t (, t) I Check! P P V t π 54 t = 0 + V (, t) π θ : = θ : π ) θ = π ± sin ± cos t = 0 (, 0) = sin π V + t +V

More information

,,.,.,,.,.,.,.,,.,..,,,, i

,,.,.,,.,.,.,.,,.,..,,,, i 22 A person recognition using color information 1110372 2011 2 13 ,,.,.,,.,.,.,.,,.,..,,,, i Abstract A person recognition using color information Tatsumo HOJI Recently, for the purpose of collection of

More information

Vol. 44 No. SIG 9(CVIM 7) ) 2) 1) 1 2) 3 7) 1) 2) 3 3) 4) 5) (a) (d) (g) (b) (e) (h) No Convergence? End (f) (c) Yes * ** * ** 1

Vol. 44 No. SIG 9(CVIM 7) ) 2) 1) 1 2) 3 7) 1) 2) 3 3) 4) 5) (a) (d) (g) (b) (e) (h) No Convergence? End (f) (c) Yes * ** * ** 1 Vol. 44 No. SIG 9(CVIM 7) July 2003, Robby T. Tan, 1 Estimating Illumination Position, Color and Surface Reflectance Properties from a Single Image Kenji Hara,, Robby T. Tan, Ko Nishino, Atsushi Nakazawa,

More information

untitled

untitled 2009 57 2 393 411 c 2009 1 1 1 2009 1 15 7 21 7 22 1 1 1 1 1 1 1 1. 1 1 1 2 3 4 12 2000 147 31 1 3,941 596 1 528 1 372 1 1 1.42 350 1197 1 13 1 394 57 2 2009 1 1 19 2002 2005 4.8 1968 5 93SNA 6 12 1 7,

More information

MUFFIN3

MUFFIN3 MUFFIN - MUltiFarious FIeld simulator for Non-equilibrium system - ( ) MUFFIN WG3 - - JCII, - ( ) - ( ) - ( ) - (JSR) - - MUFFIN sec -3 msec -6 sec GOURMET SUSHI MUFFIN -9 nsec PASTA -1 psec -15 fsec COGNAC

More information

TOP URL 1

TOP URL   1 TOP URL http://amonphys.web.fc.com/ 3.............................. 3.............................. 4.3 4................... 5.4........................ 6.5........................ 8.6...........................7

More information

Hansen 1 2, Skinner 5, Augustinus 6, Harvey 7 Windle 8 Pels 9 1 Skinner 5 Augustinus 6 Pels 9 NL Harvey ML 11 NL

Hansen 1 2, Skinner 5, Augustinus 6, Harvey 7 Windle 8 Pels 9 1 Skinner 5 Augustinus 6 Pels 9 NL Harvey ML 11 NL HANAOKA, Shinya 1 3 Hansen1, 2 1 2 3 Hansen 2 3 4 5 2 2.1 002 Vol.5 No.4 2003 Winter 3 4 2.2 Hansen 1 2, Skinner 5, Augustinus 6, Harvey 7 Windle 8 Pels 9 1 Skinner 5 Augustinus 6 Pels 9 NL Harvey 10 2.3

More information

IA 2013 : :10722 : 2 : :2 :761 :1 (23-27) : : ( / ) (1 /, ) / e.g. (Taylar ) e x = 1 + x + x xn n! +... sin x = x x3 6 + x5 x2n+1 + (

IA 2013 : :10722 : 2 : :2 :761 :1 (23-27) : : ( / ) (1 /, ) / e.g. (Taylar ) e x = 1 + x + x xn n! +... sin x = x x3 6 + x5 x2n+1 + ( IA 2013 : :10722 : 2 : :2 :761 :1 23-27) : : 1 1.1 / ) 1 /, ) / e.g. Taylar ) e x = 1 + x + x2 2 +... + xn n! +... sin x = x x3 6 + x5 x2n+1 + 1)n 5! 2n + 1)! 2 2.1 = 1 e.g. 0 = 0.00..., π = 3.14..., 1

More information

IA hara@math.kyushu-u.ac.jp Last updated: January,......................................................................................................................................................................................

More information

Physical and Psychological Effects of Stressors in Female College Students Reizou Mita*1, Konosuke Tomabechi*1, Isao Yamaguchi*1, Naoko Soeno*1, Shuhe

Physical and Psychological Effects of Stressors in Female College Students Reizou Mita*1, Konosuke Tomabechi*1, Isao Yamaguchi*1, Naoko Soeno*1, Shuhe Physical and Psychological Effects of Stressors in Female College Students Reizou Mita*1, Konosuke Tomabechi*1, Isao Yamaguchi*1, Naoko Soeno*1, Shuhei Kobayashi*2, Mamoru Nishimuta*2, Michiyuki Shimizu*3,

More information

, 1.,,,.,., (Lin, 1955).,.,.,.,. f, 2,. main.tex 2011/08/13( )

, 1.,,,.,., (Lin, 1955).,.,.,.,. f, 2,. main.tex 2011/08/13( ) 81 4 2 4.1, 1.,,,.,., (Lin, 1955).,.,.,.,. f, 2,. 82 4.2. ζ t + V (ζ + βy) = 0 (4.2.1), V = 0 (4.2.2). (4.2.1), (3.3.66) R 1 Φ / Z, Γ., F 1 ( 3.2 ). 7,., ( )., (4.2.1) 500 hpa., 500 hpa (4.2.1) 1949,.,

More information

IHIMU Energy-Saving Principle of the IHIMU Semicircular Duct and Its Application to the Flow Field Around Full Scale Ships IHI GHG IHIMU CFD PIV IHI M

IHIMU Energy-Saving Principle of the IHIMU Semicircular Duct and Its Application to the Flow Field Around Full Scale Ships IHI GHG IHIMU CFD PIV IHI M IHIMU Energy-Saving Principle of the IHIMU Semicircular Duct and Its Application to the Flow Field Around Full Scale Ships IHI GHG IHIMU PIV IHI Marine United Inc. ( IHIMU ) has already developed several

More information

浜松医科大学紀要

浜松医科大学紀要 On the Statistical Bias Found in the Horse Racing Data (1) Akio NODA Mathematics Abstract: The purpose of the present paper is to report what type of statistical bias the author has found in the horse

More information

10 1 1 (1) (2) (3) 3 3 1 3 1 3 (4) 2 32 2 (1) 1 1

10 1 1 (1) (2) (3) 3 3 1 3 1 3 (4) 2 32 2 (1) 1 1 10 10 1 1 (1) (2) (3) 3 3 1 3 1 3 (4) 2 32 2 (1) 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 (2) 1 (3) JI S JI S JI S JI S 25 175 J AS 3 (1) 3 70 (2) (3) 100 4 (1)69 (2) (3) (4) (5) (6) (7) (8)70 (9) (10)2 (11)

More information

70 : 20 : A B (20 ) (30 ) 50 1

70 : 20 : A B (20 ) (30 ) 50 1 70 : 0 : A B (0 ) (30 ) 50 1 1 4 1.1................................................ 5 1. A............................................... 6 1.3 B............................................... 7 8.1 A...............................................

More information

Developement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Sum

Developement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Sum Developement of Plastic Collocation Method Extension of Plastic Node Method by Yukio Ueda, Member Masahiko Fujikubo, Member Masahiro Miura, Member Summary Previously, the authors developed the plastic

More information

ISSN NII Technical Report Patent application and industry-university cooperation: Analysis of joint applications for patent in the Universit

ISSN NII Technical Report Patent application and industry-university cooperation: Analysis of joint applications for patent in the Universit ISSN 1346-5597 NII Technical Report Patent application and industry-university cooperation: Analysis of joint applications for patent in the University of Tokyo Morio SHIBAYAMA, Masaharu YANO, Kiminori

More information

X線分析の進歩36 別刷

X線分析の進歩36 別刷 X X X-Ray Fluorescence Analysis on Environmental Standard Reference Materials with a Dry Battery X-Ray Generator Hideshi ISHII, Hiroya MIYAUCHI, Tadashi HIOKI and Jun KAWAI Copyright The Discussion Group

More information

高等学校学習指導要領解説 数学編

高等学校学習指導要領解説 数学編 5 10 15 20 25 30 35 5 1 1 10 1 1 2 4 16 15 18 18 18 19 19 20 19 19 20 1 20 2 22 25 3 23 4 24 5 26 28 28 30 28 28 1 28 2 30 3 31 35 4 33 5 34 36 36 36 40 36 1 36 2 39 3 41 4 42 45 45 45 46 5 1 46 2 48 3

More information

A consideration about Moral Education based on scientific recognition: Through an analyzing of one reading, which is mainly for the field about relationship with a community and a society, in Our morality

More information

c 2009 i

c 2009 i I 2009 c 2009 i 0 1 0.0................................... 1 0.1.............................. 3 0.2.............................. 5 1 7 1.1................................. 7 1.2..............................

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,. 24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)

More information