untitled

Size: px
Start display at page:

Download "untitled"

Transcription

1 24 MATLAB 1 17

2 MATLAB 18

3 MATLAB download Lena f /l / Z: matlab matlab (color.zip ) Matlab current Directry z: matlab Img=imread( z: matlab Parrots.bmp ); t b workspace 19

4 MATLAB figure(1); imagesc(img); figure(2), subplot(1,3,1), imagesc(img(:, :, 1)); subplot(1,3,2), imagesc(img(:, :, 2)); subplot(1,3,3), imagesc(img(:, :, 3)); colormap(gray); R, G, B ImgG=uint8((0.299*Img(:,:,1)+0.687*Img(:,:,2)+0.114*Img(:,:,3))/3); ( ( 2)+0 114*I ( Imwrite(ImgG, z: matlab ParrotsG.bmp, BMP ); NTSCRGB 20

5 MATLAB ( ) HSV y x 21

6 MATLAB RImgG=255-ImgG; g figure(2), subplot(1,2,1), imagesc(imgg); subplot(1,2,2),, imagesc(rimgg); g colormap( gray ); ImgG16=fix(ImgG./16).*16+16; ImgG8=fix(ImgG./32).*32+32; figure(2), subplot(1,3,1), imagesc(imgg); subplot(1,3,2), imagesc(imgg16); subplot(1,3,3), imagesc(imgg8); colormap( gray ); 22

7 MATLAB (GOG ) = Gain Offset Gamma model Gain:a ImgGA=uint8(ImgG*5); g Offset:b ImgGB=uint8(ImgG-15); Gamma:γ aveimgg=mean(mean(imgg)); ImgGR=uint8(((double(ImgG)/aveImgG).^3)*aveImgG); figure(2), subplot(2,2,1), imagesc(imgg); subplot(2,2,2), 2 2) imagesc(imgga); subplot(2,2,3), imagesc(imggb); subplot(2,2,4), imagesc(imggr); colormap( gray ); 23

8 MATLAB [xsize ysize]=size(imgg); hg=zeros(256,1); for i=1:xsize for j=1:ysize hg(imgg(i, j)+1)=hg(imgg(i, j)+1)+1; figure(2), subplot(1,2,1), imagesc(imgg); g subplot(1,2,2), bar(hg); colormap(gray); 24

9 MATLAB 25

10 MATLAB [xsize,ysize]=size(imgg); len=xsize*ysize; ImgG2=reshape(ImgG,len,1); histogram = hist(imgg2,[0:256-1]); cumhist(1)= histogram(1); for i=2:256 cumhist(i)= cumhist(i-1)+ histogram(i); for i=1:256 trans(i)=round(((256-1)/(len))*cumhist(i)); finalimage=reshape(newimage, xsize, ysize); finalhistogram=hist(newimage,[0:256-1]); figure(1), subplot(3,1,1),stem([0:256-1],histogram) subplot(3,1,2),stem([0:256-1],trans) subplot(3,1,3),stem([0:256-1],finalhistogram) figure(2), subplot(1,2,1),image(imgg),colormap('gray') subplot(1,2,2),image(finalimage),colormap('gray') newimage=zeros(len,1); for i=2:256 transx(i-1)=trans(i); for i=1:len if ImgG2(i)~=0 newimage(i)=transx(imgg2(i)); else newimage(i)=trans(1); 26

11 MATLAB [xsize ysize]=size(imgg); for i=1:xsize for j=1:ysize hg((i-1)*i+j)=double(imgg(i, j)); figure(2), subplot(1,2,1), imagesc(imgg); subplot(1,2,2), hist(hg, 16); colormap(gray); 16 27

12 MATLAB 2 [xsize ysize]=size(imgg); avg= mean(mean(imgg(i, ( ( j))); for i=1:xsize for j=1:ysize if (ImgG(i, j) > avg) ImgBw(i, j)=0; else ImgBw(i, j)=255; figure(2), subplot(1,2,1), imagesc(imgg); subplot(1,2,2), imagesc(imgbw); colormap( gray ); 28

13 MATLAB 29

14 MATLAB ( ) clear all Img=imread( z: MATLAB Parrots.bmp'); ImgG=uint8((0.299*Img(:,:,1)+0.687*Img(:,:,2)+0.114*Img(:,:,3))/3); 30

15 MATLAB Mag=0.25; [xsize,ysize]=size(imgg); X=zeros(2*xsize,2*ysize); for x=1:2*xsize for y=1:2*ysize oldx=round((x-xsize)/mag); xsize)/mag); oldy=round((y-ysize)/mag); oldx=oldx+round(xsize/2); oldy=oldy+round(ysize/2); if oldx>0 & oldx<=xsize & oldy>0 & oldy<=ysize; y X(x,y)=ImgG(oldx,oldy); ScaledImgG=X(xsize-round(xsize*Mag/2)+1: xsize+round(xsize*mag/2),... ysize-round(ysize*mag/2)+1: ysize+round(ysize*mag/2)); figure(2), subplot(1,2,1),, imagesc(imgg);axis([0 g ([ xsize 0 ysize]); subplot(1,2,2), imagesc(uint8(scaledimgg));axis([0 xsize 0 ysize]); colormap('gray'); figure(3), subplot(1,2,1), imagesc(imgg);axis([0 xsize 0 ysize]); subplot(1,2,2), imagesc(uint8(scaledimgg));axis([0 xsize*mag 0 ysize*mag]); colormap('gray'); 31

16 MATLAB 'linear' 'nearest' 'spline' 'cubic' ' [x,y] = meshgrid(-3:1:3); z = peaks(x,y); figure(1), surf(x,y,z) [xi,yi] = meshgrid(-3:0.25:3); zi1=interp2(xyzxiyi'linear'); interp2(x,y,z,xi,yi, zi2 = interp2(x,y,z,xi,yi,'nearest'); zi3 = interp2(x,y,z,xi,yi,'spline'); zi4 = interp2(xyzxiyi'cubic'); interp2(x,y,z,xi,yi, figure(2), subplot(2,2,1), surf(xi,yi,zi1) subplot(2,2,2), surf(xi,yi,zi2) subplot(2,2,3), surf(xi,yi,zi3) subplot(2,2,4), surf(xi,yi,zi4) 32

17 MATLAB 'linear' 'nearest' 'spline' 'cubic' ' [xsize ysize]=size(scaledimgg); [x,y] = meshgrid(1:1:xsize); z = ScaledImgG; figure(1), imagesc(z); colormap('gray'); [xi,yi] = meshgrid(1:0.25:xsize); zi1 = interp2(x,y,z,xi,yi,'linear'); i i ') zi2 = interp2(x,y,z,xi,yi,'nearest'); zi3 = interp2(x,y,z,xi,yi,'spline'); zi4 = interp2(x,y,z,xi,yi,'cubic'); figure(2), subplot(2,2,1),imagesc(zi1); subplot(2,2,2),imagesc(zi2); subplot(2,2,3),imagesc(zi3); subplot(2,2,4),imagesc(zi4); t(224)i ( i4) colormap('gray'); 33

18 MATLAB, 34

19 MATLAB angle=pi/3; xmove=0; ymove=0; [xsize, ysize]=size(imgg); matrix=[cos(angle) sin(angle) xmove; -sin(angle) cos(angle) ymove; 0 0 1]; AffineImgG = zeros(xsize); for i=1:xsize for j=1:ysize A=[i-round(xsize/2) j-round(ysize/2) 1]; origx=matrix*a'; origx=round(origx+[xsize/2; ysize/2; 1]); if origx(1)>0 & origx(1)<=xsize & origx(2)>0 & origx(2)<=ysize; AffineImgG(i, j)=imgg(origx(1), origx(2)); figure(2), subplot(1,2,1), imagesc(imgg); subplot(1,2,2), imagesc(affineimgg); colormap('gray'); 35

20 MATLAB Tsai pin-hole f - k - Cx, Cy Sx, Sy Xu,Yu Xd, Yd X 2 2 d = X u ( 1+ kr ) Yd = Yu ( 1+ kr ) 2 2 r = X d + Y d dx, d dy S X X + f x d y d = Cx Y f = + C y d x d y S Y 36

21 MATLAB ImgG=repmat([ones(32) zeros(32); zeros(32) ones(32)],4,4); cx=0.5; cy=0.5; k=-0.2; [xsize,ysize]=size(imgg); distimgg=zeros(xsize*2,ysize*2); rmax=sqrt((xsize/2)^2+(ysize/2)^2); 2+(ysize/2) 2); for x=1:2*xsize for y=1:2*ysize rdd=((x-cx*ysize)^2+(y-cy*xsize)^2); rd2_norm=rdd/rmax^2; kr2=k*rd2_norm; oldx=round(x+(x-cx*ysize)*kr2); oldy=round(y+(y-cy*xsize)*kr2); if oldx>0 & oldx<=xsize & oldy>0 & oldy<=ysize; distimgg(x,y)=imgg(oldx,oldy); figure(2), subplot(1,2,1), 1) imagesc(imgg);axis([0 xsize 0 ysize]); subplot(1,2,2), imagesc(uint8(distimgg));axis([0 xsize 0 ysize]); colormap('gray'); 37

22 MATLAB y/pukiwiki/index.php?%b8%a6%b5%e6%2fgml%20c%2b%2b%20camera%20php?%b8%a6%b5%e6%2fgml%20c%2b%2b%20camera%20 Calibration%20Toolbox 38

untitled

untitled 25 MATLAB 2 39 ( ) HSV y x 40 4 f (t) f ( t) = t f ( t + t) tt f ( t) 2 f ( t) g( t) = f ( τ ) g( t τ ) dτ f (τ ) t g(t τ ) 0 42 t f (τ ) f (τ ) g( t τ ) 0 g( t τ ) t 3 3 3 High pass filter Low pass filter

More information

ERATO100913

ERATO100913 ERATO September 13, 2010, DC2 1/25 1. 2 2. 2/25 3/25 3/25 2 3/25 2 3/25 1 1 0.5 0.5 0 0 0.5 1 0 0 0.5 1 4/25 1 1 0.5 0.5 0 0 0.5 1 (0, 0) 0 0 0.5 1 4/25 1 1 0.5 0.5 0 0 0.5 1 (0, 0) ( 1, 0) 0 0 0.5 1 4/25

More information

untitled

untitled 20 7 1 22 7 1 1 2 3 7 8 9 10 11 13 14 15 17 18 19 21 22 - 1 - - 2 - - 3 - - 4 - 50 200 50 200-5 - 50 200 50 200 50 200 - 6 - - 7 - () - 8 - (XY) - 9 - 112-10 - - 11 - - 12 - - 13 - - 14 - - 15 - - 16 -

More information

untitled

untitled 19 1 19 19 3 8 1 19 1 61 2 479 1965 64 1237 148 1272 58 183 X 1 X 2 12 2 15 A B 5 18 B 29 X 1 12 10 31 A 1 58 Y B 14 1 25 3 31 1 5 5 15 Y B 1 232 Y B 1 4235 14 11 8 5350 2409 X 1 15 10 10 B Y Y 2 X 1 X

More information

Processingをはじめよう

Processingをはじめよう Processing をはじめよう 第 7 章 動きその 2 目次 フレームレート スピードと方向 移動 回転 拡大 縮小 2 点間の移動 乱数 タイマー 円運動 今回はここまで 2 2 点間の移動 Example 7-6 (EX_08_06) 始点 (startx, starty) から終点 (stopx, stopy) まで移動する 座標更新の計算方法は後述 始点と終点を変更しても動作する 変更して確認

More information

II - ( 02 ) 1,,,, 2, 3. ( ) HP,. 2 MATLAB MATLAB, C Java,,., MATLAB, Workspace, Workspace. Workspace who. whos. MATLAB, MATLAB Workspace. 2.1 Workspac

II - ( 02 ) 1,,,, 2, 3. ( ) HP,. 2 MATLAB MATLAB, C Java,,., MATLAB, Workspace, Workspace. Workspace who. whos. MATLAB, MATLAB Workspace. 2.1 Workspac II - ( 02 ) 1,,,, 2, 3 ( ) HP, 2 MATLAB MATLAB, C Java,,, MATLAB, Workspace, Workspace Workspace who whos MATLAB, MATLAB Workspace 21 Workspace 211 Workspace save, Workspace, MATLAB MAT, load, MAT Workspace

More information

応力とひずみ.ppt

応力とひずみ.ppt in yukawa@numse.nagoya-u.ac.jp 2 3 4 5 x 2 6 Continuum) 7 8 9 F F 10 F L L F L 1 L F L F L F 11 F L F F L F L L L 1 L 2 12 F L F! A A! S! = F S 13 F L L F F n = F " cos# F t = F " sin# S $ = S cos# S S

More information

2 T ax 2 + 2bxy + cy 2 + dx + ey + f = 0 a + b + c > 0 a, b, c A xy ( ) ( ) ( ) ( ) u = u 0 + a cos θ, v = v 0 + b sin θ 0 θ 2π u = u 0 ± a

2 T ax 2 + 2bxy + cy 2 + dx + ey + f = 0 a + b + c > 0 a, b, c A xy ( ) ( ) ( ) ( ) u = u 0 + a cos θ, v = v 0 + b sin θ 0 θ 2π u = u 0 ± a 2 T140073 1 2 ax 2 + 2bxy + cy 2 + dx + ey + f = 0 a + b + c > 0 a, b, c A xy u = u 0 + a cos θ, v = v 0 + b sin θ 0 θ 2π u = u 0 ± a cos θ, v = v 0 + b tan θ π 2 < θ < π 2 u = u 0 + 2pt, v = v 0 + pt

More information

HPC146

HPC146 2 3 4 5 6 int array[16]; #pragma xmp nodes p(4) #pragma xmp template t(0:15) #pragma xmp distribute t(block) on p #pragma xmp align array[i] with t(i) array[16] 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Node

More information

1.. 1 ll a ii. 1i. i f 1 1 a. a. i. t. 1 fi fi. t i j fj i. j ;i 1. i. aa a

1.. 1 ll a ii. 1i. i f 1 1 a. a. i. t. 1 fi fi. t i j fj i. j ;i 1. i. aa a 1.. 1 ll a ii. 1i. i f 1 1 a. a. i. t. 1 fi fi. t i j fj i. j ;i 1. i. aa 1 111 0 0 0 0 a I E l21 1fi i L < i i;i1=t ii 111 1; ai i ti a t T ;,, l 1i.... E 11fi i 1t l l t2 1i i1 t Ea li )2 0 u 0 1f )2

More information

Introduction to Matlab Programming with Applications

Introduction to Matlab Programming with Applications Histogram ˆ The function hist(x) or histogram creates a histogram of x. ˆ By default, elements in x are organized into equally spaced bins between min(x) and max(x). ˆ You can specify the number of bins.

More information

1 2 2 (POC)

1 2 2 (POC) 2007 Taichi Kawasaki Kinki University School of Engineerings 1 2 2 (POC) 6 3 7 4 12 5 13 6 18 7 24 8 26 9 27 10 31 1 1 [1]. 1 [2] P.Chazal [3] 2 3 2 1: 3 2: 4 3: 5 2 (POC) f(n 1, n 2 ) g(n 1, n 2 ) 2 F

More information

mugensho.dvi

mugensho.dvi 1 1 f (t) lim t a f (t) = 0 f (t) t a 1.1 (1) lim(t 1) 2 = 0 t 1 (t 1) 2 t 1 (2) lim(t 1) 3 = 0 t 1 (t 1) 3 t 1 2 f (t), g(t) t a lim t a f (t) g(t) g(t) f (t) = o(g(t)) (t a) = 0 f (t) (t 1) 3 1.2 lim

More information

cpp4.dvi

cpp4.dvi 2017 c 4 C++ (4) C++, 41, 42, 1, 43,, 44 45, 41 (inheritance),, C++,, 100, 50, PCMCIA,,,,,,,,, 42 1 421 ( ), car 1 [List 41] 1: class car { 2: private: 3: std::string m model; // 4: std::string m maker;

More information

2004 2005 2 2 1G01P038-0 1 2 1.1.............................. 2 1.2......................... 2 1.3......................... 3 2 4 2.1............................ 4 2.2....................... 4 2.3.......................

More information

数値計算:常微分方程式

数値計算:常微分方程式 ( ) 1 / 82 1 2 3 4 5 6 ( ) 2 / 82 ( ) 3 / 82 C θ l y m O x mg λ ( ) 4 / 82 θ t C J = ml 2 C mgl sin θ θ C J θ = mgl sin θ = θ ( ) 5 / 82 ω = θ J ω = mgl sin θ ω J = ml 2 θ = ω, ω = g l sin θ = θ ω ( )

More information

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2 II No.1 [n/] [1]H n x) H n x) = 1) r n! r!n r)! x)n r r= []H n x) n,, H n x) = 1) n H n x) [3] H n x) = 1) n dn x e dx n e x [4] H n+1 x) = xh n x) nh n 1 x) ) d dx x H n x) = H n+1 x) d dx H nx) = nh

More information

Sample Input Output for the Sample Input Sample Input Output for the Sample Input Sample Input 4-1 Output fo

Sample Input Output for the Sample Input Sample Input Output for the Sample Input Sample Input 4-1 Output fo Problem A Input: Standard Input Japan Alumni Group Summer Camp 2014 Day 4, AtCoder 15 Sep 2014 A A a b c A 60 A 60t + c ( t ) A -1 1 a b c Input 3 a, b, c a b c a, b, c 0 < a, b, c < 60 Constraints Output

More information

O1-1 O1-2 O1-3 O1-4 O1-5 O1-6

O1-1 O1-2 O1-3 O1-4 O1-5 O1-6 O1-1 O1-2 O1-3 O1-4 O1-5 O1-6 O1-7 O1-8 O1-9 O1-10 O1-11 O1-12 O1-13 O1-14 O1-15 O1-16 O1-17 O1-18 O1-19 O1-20 O1-21 O1-22 O1-23 O1-24 O1-25 O1-26 O1-27 O1-28 O1-29 O1-30 O1-31 O1-32 O1-33 O1-34 O1-35

More information

1 matlab 1 1.1.................................................... 1 1.1.1............................................... 2 1.2.......................

1 matlab 1 1.1.................................................... 1 1.1.1............................................... 2 1.2....................... : : 3 (2004/11/2) : MATLAB (2004/11/9) : (2004/11/16) : 2004/11/30 : 16 10 26 : 15 12 29 : : 1 matlab 1 1.1.................................................... 1 1.1.1...............................................

More information

[] Tle () i ( ) e . () [].....[], ..i.et.. N Z, I Q R C N Z Q R C R i R {,,} N A B X N Z Q R,, c,,, c, A, B, C, L, y, z,, X, L pq p q def f () lim f ( ) f ( ) ( ), p p q r q r p q r p q r c c,, f ( )

More information

untitled

untitled 16 Edge emphasis filter of images that uses Histogram 1050307 2005 2 18 RGB 9,, i Abstract Edge emphasis filter of images that uses Histogram Noriaki Okamoto Theedgeemphasisfilter is usually used to recognize

More information

A

A A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................

More information

web05.dvi

web05.dvi 5 3 copyright c Tatsuya Kitamura / All rights reserved. 57 5 5. 3 MATLAB 3 line plot plot3 MAT plot Z Octave gnuplot splot gsplot OCT 3 splot plot3 MAT plot 5. 5. 5.3 plot 5. >> t=:pi/:*pi; >> plot3(t.*cos(t),t.*sin(t),t)

More information

dx dt = f x,t ( ) t

dx dt = f x,t ( ) t MATLAB 1. 2. Runge-Kutta 3. 1. 2. 3. dx dt = f x,t ( ) t dx( t) dt = lim Δt x( t + Δt) x( t) Δt t = nδt n = 0,1,2,3,4,5, x = x( nδt) n Δt dx ( ) dt = f x,t x n +1 x n Δt = f ( x,nδt) n 1 x n = x 0 n =

More information

) Binary Cubic Forms / 25

) Binary Cubic Forms / 25 2016 5 2 ) Binary Cubic Forms 2016 5 2 1 / 25 1 2 2 2 3 2 3 ) Binary Cubic Forms 2016 5 2 2 / 25 1.1 ( ) 4 2 12 = 5+7, 16 = 5+11, 36 = 7+29, 1.2 ( ) p p+2 3 5 5 7 11 13 17 19, 29 31 41 43 ) Binary Cubic

More information

1 1 1 1 1 1 2 f z 2 C 1, C 2 f 2 C 1, C 2 f(c 2 ) C 2 f(c 1 ) z C 1 f f(z) xy uv ( u v ) = ( a b c d ) ( x y ) + ( p q ) (p + b, q + d) 1 (p + a, q + c) 1 (p, q) 1 1 (b, d) (a, c) 2 3 2 3 a = d, c = b

More information

x ( ) x dx = ax

x ( ) x dx = ax x ( ) x dx = ax 1 dx = a x log x = at + c x(t) = e at C (C = e c ) a > 0 t a < 0 t 0 (at + b ) h dx = lim x(t + h) x(t) h 0 h x(t + h) x(t) h x(t) t x(t + h) x(t) ax(t) h x(t + h) x(t) + ahx(t) 0, h, 2h,

More information

ax 2 + bx + c = n 8 (n ) a n x n + a n 1 x n a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 4

ax 2 + bx + c = n 8 (n ) a n x n + a n 1 x n a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 4 20 20.0 ( ) 8 y = ax 2 + bx + c 443 ax 2 + bx + c = 0 20.1 20.1.1 n 8 (n ) a n x n + a n 1 x n 1 + + a 1 x + a 0 = 0 ( a n, a n 1,, a 1, a 0 a n 0) n n ( ) ( ) ax 3 + bx 2 + cx + d = 0 444 ( a, b, c, d

More information

II (10 4 ) 1. p (x, y) (a, b) ε(x, y; a, b) 0 f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a 2 x a 0 A = f x (

II (10 4 ) 1. p (x, y) (a, b) ε(x, y; a, b) 0 f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a 2 x a 0 A = f x ( II (1 4 ) 1. p.13 1 (x, y) (a, b) ε(x, y; a, b) f (x, y) f (a, b) A, B (6.5) y = b f (x, b) f (a, b) x a = A + ε(x, b; a, b) x a x a A = f x (a, b) y x 3 3y 3 (x, y) (, ) f (x, y) = x + y (x, y) = (, )

More information

(pdf) (cdf) Matlab χ ( ) F t

(pdf) (cdf) Matlab χ ( ) F t (, ) (univariate) (bivariate) (multi-variate) Matlab Octave Matlab Matlab/Octave --...............3. (pdf) (cdf)...3.4....4.5....4.6....7.7. Matlab...8.7.....9.7.. χ ( )...0.7.3.....7.4. F....7.5. t-...3.8....4.8.....4.8.....5.8.3....6.8.4....8.8.5....8.8.6....8.9....9.9.....9.9.....0.9.3....0.9.4.....9.5.....0....3

More information

M3 x y f(x, y) (= x) (= y) x + y f(x, y) = x + y + *. f(x, y) π y f(x, y) x f(x + x, y) f(x, y) lim x x () f(x,y) x 3 -

M3 x y f(x, y) (= x) (= y) x + y f(x, y) = x + y + *. f(x, y) π y f(x, y) x f(x + x, y) f(x, y) lim x x () f(x,y) x 3 - M3............................................................................................ 3.3................................................... 3 6........................................... 6..........................................

More information

2008 DS T050049

2008 DS T050049 DS T050049. PSP DS DS DS RPG DS OS Windows XP DevkiPro OS DS CPU ARM devkitarm MSYS MinGW MSYS MinGW Unix OS C++ C++ make nds nds DS DS micro SD Card nds DS DS DS nds C Java C++ nds nds DS 2008 DS T050049

More information

XACCの概要

XACCの概要 2 global void kernel(int a[max], int llimit, int ulimit) {... } : int main(int argc, char *argv[]){ MPI_Int(&argc, &argc); MPI_Comm_rank(MPI_COMM_WORLD, &rank); MPI_Comm_size(MPI_COMM_WORLD, &size); dx

More information

II 2014 2 (1) log(1 + r/100) n = log 2 n log(1 + r/100) = log 2 n = log 2 log(1 + r/100) (2) y = f(x) = log(1 + x) x = 0 1 f (x) = 1/(1 + x) f (0) = 1

II 2014 2 (1) log(1 + r/100) n = log 2 n log(1 + r/100) = log 2 n = log 2 log(1 + r/100) (2) y = f(x) = log(1 + x) x = 0 1 f (x) = 1/(1 + x) f (0) = 1 II 2014 1 1 I 1.1 72 r 2 72 8 72/8 = 9 9 2 a 0 1 a 1 a 1 = a 0 (1+r/100) 2 a 2 a 2 = a 1 (1 + r/100) = a 0 (1 + r/100) 2 n a n = a 0 (1 + r/100) n a n a 0 2 n a 0 (1 + r/100) n = 2a 0 (1 + r/100) n = 2

More information

1 1 Pixel 0 n 1 n=8 56 R G B RGB M RGB (1) M = 0.99R G B (1) () 4 π d 4 B = L cos φ () 4 ID B L d ID φ d / ID F R φ (3) R

1 1 Pixel 0 n 1 n=8 56 R G B RGB M RGB (1) M = 0.99R G B (1) () 4 π d 4 B = L cos φ () 4 ID B L d ID φ d / ID F R φ (3) R I-07 404 137, Email:dsusuki@ipc.shizuoka.ac.jp When preparing the manuscript, read and observe carefully this sample as well as the instruction manual for the manuscript (1) OS of the Transaction of Japan

More information

9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x

9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x 2009 9 6 16 7 1 7.1 1 1 1 9 2 1 f(x, y) = xy sin x cos y x y cos y y x sin x d (x, y) = y cos y (x sin x) = y cos y(sin x + x cos x) x dx d (x, y) = x sin x (y cos y) = x sin x(cos y y sin y) y dy 1 sin

More information

31 4 MATLAB A, B R 3 3 A = , B = mat_a, mat_b >> mat_a = [-1, -2, -3; -4, -5, -6; -7, -8, -9] mat_a =

31 4 MATLAB A, B R 3 3 A = , B = mat_a, mat_b >> mat_a = [-1, -2, -3; -4, -5, -6; -7, -8, -9] mat_a = 3 4 MATLAB 3 4. A, B R 3 3 2 3 4 5 6 7 8 9, B = mat_a, mat_b >> mat_a = [-, -2, -3; -4, -5, -6; -7, -8, -9] 9 8 7 6 5 4 3 2 mat_a = - -2-3 -4-5 -6-7 -8-9 >> mat_b = [-9, -8, -7; -6, -5, -4; -3, -2, -]

More information

2

2 16 1050026 1050042 1 2 1 1.1 3 1.2 3 1.3 3 2 2.1 4 2.2 4 2.2.1 5 2.2.2 5 2.3 7 2.3.1 1Basic 7 2.3.2 2 8 2.3.3 3 9 2.3.4 4window size 10 2.3.5 5 11 3 3.1 12 3.2 CCF 1 13 3.3 14 3.4 2 15 3.5 3 17 20 20 20

More information

ARMSによる化学反応のシミュレーション

ARMSによる化学反応のシミュレーション 2001 3 Brusselator BrusselatorMulti Multi-set ARMSAbstruct Rewriting Multiset System BZ BrusselatorARMS Belousov-Zhabotinsky BZ Belousov-Zhabotinsky BZ B.P.Belousov 1951 A.M.Zhabotinsky Belousov 2 Belousov-Zhabotinsky

More information

2 MATLABについて n MATLABとは 科 学 技 術 計 算 のための 高 性 能 プログラミング 言 語 n 特 徴 配 列 が 基 本 データ 型 ベクトル(1 次 元 配 列 ) 行 列 (2 次 元 配 列 ) 対 話 的 システム 豊 富 な 関 数 ライブラリとグラフィックツー

2 MATLABについて n MATLABとは 科 学 技 術 計 算 のための 高 性 能 プログラミング 言 語 n 特 徴 配 列 が 基 本 データ 型 ベクトル(1 次 元 配 列 ) 行 列 (2 次 元 配 列 ) 対 話 的 システム 豊 富 な 関 数 ライブラリとグラフィックツー 2016 年 5 月 31 日 版 MATLABの 基 本 的 な 使 い 方 担 当 : 荻 田 武 史 2 MATLABについて n MATLABとは 科 学 技 術 計 算 のための 高 性 能 プログラミング 言 語 n 特 徴 配 列 が 基 本 データ 型 ベクトル(1 次 元 配 列 ) 行 列 (2 次 元 配 列 ) 対 話 的 システム 豊 富 な 関 数 ライブラリとグラフィックツール

More information

R R 16 ( 3 )

R R 16   ( 3 ) (017 ) 9 4 7 ( ) ( 3 ) ( 010 ) 1 (P3) 1 11 (P4) 1 1 (P4) 1 (P15) 1 (P16) (P0) 3 (P18) 3 4 (P3) 4 3 4 31 1 5 3 5 4 6 5 9 51 9 5 9 6 9 61 9 6 α β 9 63 û 11 64 R 1 65 13 66 14 7 14 71 15 7 R R 16 http://wwwecoosaka-uacjp/~tazak/class/017

More information

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

untitled

untitled yoshi@image.med.osaka-u.ac.jp http://www.image.med.osaka-u.ac.jp/member/yoshi/ II Excel, Mathematica Mathematica Osaka Electro-Communication University (2007 Apr) 09849-31503-64015-30704-18799-390 http://www.image.med.osaka-u.ac.jp/member/yoshi/

More information

第13回:交差項を含む回帰・弾力性の推定

第13回:交差項を含む回帰・弾力性の推定 13 2018 7 27 1 / 31 1. 2. 2 / 31 y i = β 0 + β X x i + β Z z i + β XZ x i z i + u i, E(u i x i, z i ) = 0, E(u i u j x i, z i ) = 0 (i j), V(u i x i, z i ) = σ 2, i = 1, 2,, n x i z i 1 3 / 31 y i = β

More information

Microsoft Word - 計算力学2007有限要素法.doc

Microsoft Word - 計算力学2007有限要素法.doc 95 2 x y yz = zx = yz = zx = { } T = { x y z xy } () {} T { } T = { x y z xy } = u u x y u z u x x y z y + u y (2) x u x u y x y x y z xy E( ) = ( + )( 2) 2 2( ) x y z xy (3) E x y z z = z = (3) z x y

More information

2 2.1 ( ) >> x=linspace(255,0,12) 2 x =

2 2.1 ( ) >> x=linspace(255,0,12) 2 x = 2 2.1 ( ) 2.1 1 >> x=linspace(255,0,12) 2 x = 3 1 5 4 255.0000 231.8182 208.6364 185.4545 162.2727 5 6 10 6 139.0909 115.9091 92.7273 69.5455 46.3636 7 11 12 8 23.1818 0 9 >> x=uint8(x) 10 x = 11 1 12

More information

70の法則

70の法則 70 70 1 / 27 70 1 2 3 4 5 6 2 / 27 70 70 70 X r % = 70 2 r r r 10 72 70 72 70 : 1, 2, 5, 7, 10, 14, 35, 70 72 : 1, 2, 3, 4, 6, 8, 9, 12, 18, 24, 36, 72 3 / 27 r = 10 70 r = 10 70 1 : X, X 10 = ( X + X

More information

L P y P y + ɛ, ɛ y P y I P y,, y P y + I P y, 3 ŷ β 0 β y β 0 β y β β 0, β y x x, x,, x, y y, y,, y x x y y x x, y y, x x y y {}}{,,, / / L P / / y, P

L P y P y + ɛ, ɛ y P y I P y,, y P y + I P y, 3 ŷ β 0 β y β 0 β y β β 0, β y x x, x,, x, y y, y,, y x x y y x x, y y, x x y y {}}{,,, / / L P / / y, P 005 5 6 y β + ɛ {x, x,, x p } y, {x, x,, x p }, β, ɛ E ɛ 0 V ɛ σ I 3 rak p 4 ɛ i N 0, σ ɛ ɛ y β y β y y β y + β β, ɛ β y + β 0, β y β y ɛ ɛ β ɛ y β mi L y y ŷ β y β y β β L P y P y + ɛ, ɛ y P y I P y,,

More information

1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2

1.2 y + P (x)y + Q(x)y = 0 (1) y 1 (x), y 2 (x) y 1 (x), y 2 (x) (1) y(x) c 1, c 2 y(x) = c 1 y 1 (x) + c 2 y 2 (x) 3 y 1 (x) y 1 (x) e R P (x)dx y 2 1 1.1 R(x) = 0 y + P (x)y + Q(x)y = R(x)...(1) y + P (x)y + Q(x)y = 0...(2) 1 2 u(x) v(x) c 1 u(x)+ c 2 v(x) = 0 c 1 = c 2 = 0 c 1 = c 2 = 0 2 0 2 u(x) v(x) u(x) u (x) W (u, v)(x) = v(x) v (x) 0 1 1.2

More information

°ÌÁê¿ô³ØII

°ÌÁê¿ô³ØII July 14, 2007 Brouwer f f(x) = x x f(z) = 0 2 f : S 2 R 2 f(x) = f( x) x S 2 3 3 2 - - - 1. X x X U(x) U(x) x U = {U(x) x X} X 1. U(x) A U(x) x 2. A U(x), A B B U(x) 3. A, B U(x) A B U(x) 4. A U(x),

More information

ユニセフ表紙_CS6_三.indd

ユニセフ表紙_CS6_三.indd 16 179 97 101 94 121 70 36 30,552 1,042 100 700 61 32 110 41 15 16 13 35 13 7 3,173 41 1 4,700 77 97 81 47 25 26 24 40 22 14 39,208 952 25 5,290 71 73 x 99 185 9 3 3 3 8 2 1 79 0 d 1 226 167 175 159 133

More information

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x [ ] IC. f(x) = e x () f(x) f (x) () lim f(x) lim f(x) x + x (3) lim f(x) lim f(x) x + x (4) y = f(x) ( ) ( s46). < a < () a () lim a log xdx a log xdx ( ) n (3) lim log k log n n n k=.3 z = log(x + y ),

More information

Assignment_.java 課題 : 転置行列 / class Assignment_ public static void main(string[] args) int i,j; int[][] array = 1,,,,,,,,,,,,,1,1,; 行 列行列 i

Assignment_.java 課題 : 転置行列 / class Assignment_ public static void main(string[] args) int i,j; int[][] array = 1,,,,,,,,,,,,,1,1,; 行 列行列 i 1 1 0 1 Assignment_1.java 課題 1: チェッカー / class Assignment_1 public static void main(string[] args) int i,j; チェッカー用の 次元配列 int[][] checker=new int[][]; チェッカーパターンを書き込む for(i=0;i

More information

IPSJ SIG Technical Report Vol.2014-DPS-158 No.27 Vol.2014-CSEC-64 No /3/6 1,a) 2,b) 3,c) 1,d) 3 Cappelli Bazen Cappelli Bazen Cappelli 1.,,.,.,

IPSJ SIG Technical Report Vol.2014-DPS-158 No.27 Vol.2014-CSEC-64 No /3/6 1,a) 2,b) 3,c) 1,d) 3 Cappelli Bazen Cappelli Bazen Cappelli 1.,,.,., 1,a),b) 3,c) 1,d) 3 Cappelli Bazen Cappelli Bazen Cappelli 1.,,,,,.,,,,.,,.,,,,.,, 1 Department of Electrical Electronic and Communication Engineering Faculty of Science and Engineering Chuo University

More information

( )$("canvas").drawarc({strokestyle:"red", x:100, y:100, radius:20, start:0, end:360); drawline(x1:, y1:,... xn:, yn:) drawline n 2 n 3 x1: y1: xn: yn

( )$(canvas).drawarc({strokestyle:red, x:100, y:100, radius:20, start:0, end:360); drawline(x1:, y1:,... xn:, yn:) drawline n 2 n 3 x1: y1: xn: yn 0.1. jquaery jcanvas jquery John Resig OSS JavaScript Web JavaScript jquery jquery 1 0.2. jcanvas jcanvas 0.3. jcanvas HTML5 Canvas Canvas jcanvas jcanvas jquery Canvas API jcanvas Grouping Layer jcanvas

More information

untitled

untitled . 96. 99. ( 000 SIC SIC N88 SIC for Windows95 6 6 3 0 . amano No.008 6. 6.. z σ v σ v γ z (6. σ 0 (a (b 6. (b 0 0 0 6. σ σ v σ σ 0 / v σ v γ z σ σ 0 σ v 0γ z σ / σ ν /( ν, ν ( 0 0.5 0.0 0 v sinφ, φ 0 (6.

More information

Acrobat Distiller, Job 128

Acrobat Distiller, Job 128 (2 ) 2 < > ( ) f x (x, y) 2x 3+y f y (x, y) x 2y +2 f(3, 2) f x (3, 2) 5 f y (3, 2) L y 2 z 5x 5 ` x 3 z y 2 2 2 < > (2 ) f(, 2) 7 f x (x, y) 2x y f x (, 2),f y (x, y) x +4y,f y (, 2) 7 z (x ) + 7(y 2)

More information

分散分析・2次元正規分布

分散分析・2次元正規分布 2 II L10(2016-06-30 Thu) : Time-stamp: 2016-06-30 Thu 13:55 JST hig F 2.. http://hig3.net ( ) L10 2 II(2016) 1 / 24 F 2 F L09-Q1 Quiz :F 1 α = 0.05, 2 F 3 H 0, : σ 2 1 /σ2 2 = 1., H 1, σ 2 1 /σ2 2 1. 4

More information

Copyright c 2006 Zhenjiang Hu, All Right Reserved.

Copyright c 2006 Zhenjiang Hu, All Right Reserved. 1 2006 Copyright c 2006 Zhenjiang Hu, All Right Reserved. 2 ( ) 3 (T 1, T 2 ) T 1 T 2 (17.3, 3) :: (Float, Int) (3, 6) :: (Int, Int) (True, (+)) :: (Bool, Int Int Int) 4 (, ) (, ) :: a b (a, b) (,) x y

More information

26 No.36 2001.1

26 No.36 2001.1 No.36 2001.1 25 26 No.36 2001.1 No.36 2001.1 27 28 No.36 2001.1 No.36 2001.1 29 30 No.36 2001.1 No.36 2001.1 31 32 No.36 2001.1 No.36 2001.1 33 Mycobacterium bovis 34 No.36 2001.1 No.36 2001.1 35 36 No.36

More information

plotwsx PLOT-WSX X11 X spp[2]./a.out PLOT-WSX PLOT IPENS -13 IPENS 0 IPENS 0 2 s p PLOT MRI 3.2 PLOT-WSX PLOT IPE- NS -13 IPENS 0 1 Continue PLO

plotwsx PLOT-WSX X11 X spp[2]./a.out PLOT-WSX PLOT IPENS -13 IPENS 0 IPENS 0 2 s p PLOT MRI 3.2 PLOT-WSX PLOT IPE- NS -13 IPENS 0 1 Continue PLO PLOT-WSX 3 1 EWS FORTRAN PLOT-WSX FORTRAN spp 2 2.1 1 1 1 1 Y (Xmax,Ymax) (0,Ymax) 2.2 PLOT-WSX NEWPEN 1 1. NEWPEN 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 2.3 PLOT-WSX PLOTS 4.1.1 PLOT-WSX outle.ps 3 X (0,0)

More information

³ÎΨÏÀ

³ÎΨÏÀ 2017 12 12 Makoto Nakashima 2017 12 12 1 / 22 2.1. C, D π- C, D. A 1, A 2 C A 1 A 2 C A 3, A 4 D A 1 A 2 D Makoto Nakashima 2017 12 12 2 / 22 . (,, L p - ). Makoto Nakashima 2017 12 12 3 / 22 . (,, L p

More information

1/1 lim f(x, y) (x,y) (a,b) ( ) ( ) lim limf(x, y) lim lim f(x, y) x a y b y b x a ( ) ( ) xy x lim lim lim lim x y x y x + y y x x + y x x lim x x 1

1/1 lim f(x, y) (x,y) (a,b) ( ) ( ) lim limf(x, y) lim lim f(x, y) x a y b y b x a ( ) ( ) xy x lim lim lim lim x y x y x + y y x x + y x x lim x x 1 1/5 ( ) Taylor ( 7.1) (x, y) f(x, y) f(x, y) x + y, xy, e x y,... 1 R {(x, y) x, y R} f(x, y) x y,xy e y log x,... R {(x, y, z) (x, y),z f(x, y)} R 3 z 1 (x + y ) z ax + by + c x 1 z ax + by + c y x +

More information

Underlying mechanisms of biochemical oscillations

Underlying mechanisms of biochemical oscillations Underlying mechanisms of biochemical oscillations (YUGI, Katsuyuki) Kuroda Lab., The University of Tokyo Electrocardiograph Borisuk and Tyson (1998) k1.2x10 2 sec 1 Hes1 (Notch signalling

More information

課題

課題 int starttime_msec; boolean counting = false; size(400,200); smooth(); //font は各自のものに変更してください font = loadfont("serif-48.vlw"); void mouseclicked(){ counting = true; starttime_msec = millis(); int t=0;

More information

2

2 CONTENT S 01 02 04 06 08 10 20 25 25 30 38 49 53 58 60 64 74 8 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 2 6 1 2 2002,000 200 7 Web 20 20 13 13 4 13 13F 13F 20 20 13 13 4 13 13F 13F 13 6 A B C A B 13 20 13 20

More information

a x x x x 1 x 2 Ý; x. x = x 1 + x 2 + Ý + x = 10 1; 1; 3; 3; 4; 5; 8; 8; 8; 9 1 + 1 + 3 + 3 + 4 + 5 + 8 + 8 + 8 + 9 10 = 50 10 = 5 . 1 1 Ý Ý # 2 2 Ý Ý & 7 7; 9; 15; 21; 33; 44; 56 21 8 7; 9; 15; 20; 22;

More information

DVIOUT-MTT元原

DVIOUT-MTT元原 TI-92 -MTT-Mathematics Thinking with Technology MTT ACTIVITY Discussion 1 1 1.1 v t h h = vt 1 2 gt2 (1.1) xy (5, 0) 20m/s [1] Mode Graph Parametric [2] Y= [3] Window [4] Graph 1.1: Discussion 2 Window

More information

II I Riemann 2003

II I Riemann 2003 II I Remann 2003 Dfferental Geometry II Dfferental Geometry II I Dfferental Geometry I 1 1 1.1.............................. 1 1.2................................. 10 1.3.......................... 16 1.4.........................

More information

a q q y y a xp p q y a xp y a xp y a x p p y a xp q y x yaxp x y a xp q x p y q p x y a x p p p p x p

a q q y y a xp p q y a xp y a xp y a x p p y a xp q y x yaxp x y a xp q x p y q p x y a x p p p p x p a a a a y y ax q y ax q q y ax y ax a a a q q y y a xp p q y a xp y a xp y a x p p y a xp q y x yaxp x y a xp q x p y q p x y a x p p p p x p y a xp q y a x p q p p x p p q p q y a x xy xy a a a y a x

More information