[응용수학]수치해석 실습과제

소개글

중앙대학교 정태상교수님 2005년 2학기 수치해석 실습과제

목차

1. 3D sinc 함수 구하기

2. 확률적 방법으로 pi 값 구하기

3. 파일 입출력하기

본문내용

1. 3D sinc 함수 구하기

교수님!
문제 본문에 보면 R=X^2 + Y^2일때 라고 되어 있는데 문제처럼 R값으로 r=(x.^2 + y.^2); 라고 정의하면 다음과 같이 나옵니다.


[x,y] = meshgrid(-2:0.1:2,-2:0.1:2);
r=sqrt(x.^2 + y.^2);
if x==0 & y==0
f=1;
elseif x~=0 & y~=0
f = sin(pi*r)./(pi*r);
end
mesh(x,y,f), title(`3-D SinC`)
xlabel (`x`), ylabel (`y`), zlabel (`f`)
grid


교수님께서 올려주신 캡쳐된 그림과 동일한 모양이 나와야 하기 위해서는 R값이 이 되어야 합니다. 다음과 같이 M-file을 해주면 다음과 같이 나옵니다.

※ 소스코드의 주석은 첨부한 파일에 모두 있습니다.
2. 확률적 방법으로 pi 값 구하기
N = input(`Total number of random points: `);
NInside = 0;
for nloops=1:N
Xrand = rand;
Yrand = rand; % Generate Random XY Point
Rrand = Xrand^2 + Yrand^2; % Find its distance from (0,0)
if (Rrand <= 1) % Rrand <=1 means that point exist in the circle
NInside = NInside + 1; % count inside points
end
end

참고 자료

없음