안녕하세요
외국에서 공부하고 있는 유학생입니다!!
몇 년째 한글을 안 쓰다 보니 문법이나 맞춤법이나 흐름이 이상해도 이해해주세요 ㅠㅠ
대학교 때배워 본 적도 없는 매트렙을
대학원 들어오니까 많이 쓰길래 독학 중입니다
정말.. 어렵네요 ㅠㅠ
계산중에 문제가 생겨서 고수분들에게 도움을 요청하고자 합니다!
푸리에 적분을 계산하려고 하는데요 ㅠㅠ
이게 계산 식입니다
A(w)랑 B(w)는 구할 수가 있는데 xx(t)를 적분하는 과정에서 계속 문제가 생기고 어떻게 해야 할지 모르겠어요ㅠ
제가 사용하려는 방법은 A(w)와B(w)를 구한 후 xx(t)식에 대입해서 적분을 마친 후에 xx(t) 함수를 구해내고
ezplolt를 통해서 xx(t)를 그래프로 그리려고 했어요!
clear all
clc
load signal04
l1=length(time4);
syms w v t;
xx_v=1;
a_0=int(xx_v*cos(w*v),v,3,6);
b_0=int(xx_v*sin(w*v),v,3,6);
a_1=(1/pi)*a_0;
b_1=(1/pi)*b_0;
disp(a_1)
disp(b_1)
xx_f=int(a_1*cos(w*t)+b_1*sin(w*t),w)
disp(xx_f)
이게 제가 쓴 코드에요.
함수는 이렇게 생겼고요
그래서 저는 3~6 이외는 전부 값이 0이고 3~6의 값은1이어서 xx_v=1 이렇게 하고 3에서 6을 적분을 했어요
A(w)와 B(w)의 값을 구하는 데 성공했고요!
A(w)=
B(w)=
근데 xx(t)를 구하려고 적분을 시도하면
Warning: Explicit integral could not be found.
int((5734161139222659*sin(t*w)*(cos(3*w) - cos(6*w)))/(18014398509481984*w) -
(5734161139222659*cos(t*w)*(sin(3*w) - sin(6*w)))/(18014398509481984*w), w)
이런식 으로 적분을 할 수가 없다고만 나오네요 ㅠ
적분은 고의로 부정적분을 했어요. 부정적분조차 안된다면 정적분은 더 안될거 같아서..
대체 어떤 게 문제인지 알 수 없어서
손으로 적분을 계산 해 봣는데
여기서 적분을 어떻게 해야 할지 모르겠더라고요
그래서 매트랩 으로 위 식을 계산 해봤는데 pi/4가 나오더라고요
매트랩엔 a=int((sin(w)/w)*cos(w),w,0,inf) 이렇게 썼어요
그다음엔 제가
먼저 이식으로 A(w)를 구한 후에
값이 영이 나오더라고요
syms w v;
xx_v=1;
a_0=int(xx_v*cos(w*v),v);
a_1=(1/pi)*a_0;
a=int(a_1*cos(w),w,0,inf)
매트렙에선 이렇게 썻습니다
위에 제가 구하고 싶은 식이랑 t가 하나 차이나서 cos안에 t를 w 에 곱하는 식으로 썼어요 밑에 식처럼요
매트랩에선
syms w v t;
xx_v=1;
a_0=int(xx_v*cos(w*v),v,3,6);
a_1=(1/pi)*a_0;
a=int(a_1*cos(w*t),w,0,inf)
이렇게 썼습니다!
그랬더니 결과가
Warning: Explicit integral could not be found.
a =
int(-(5734161139222659*cos(t*w)*(sin(3*w) - sin(6*w)))/(18014398509481984*w), w = 0..Inf)
이렇게 뜨네요 ㅠㅠ
이 문제로 벌써 거의 일주일째 고민 중이에요
교수님은 알려주지도 않고 선배들도 모른다고만 하고요..
그래서 자력으로 해결해 보려 했지만
여기만 일주일째 풀질 못 하니까 정말 답답하더라고요
오유에서 만큼은 이런 질문 드리기 싫었는데 괜히 민폐 끼쳐서 죄송합니다 ㅠㅠ
고수님들 도와주세요! ㅠㅠ
그리고 혹시 매트랩 토론하는 곳이나 배울 수 있는 좋은 사이트 있으면 추천 좀 부탁합니다!
감사합니다!!
창작글
본인삭제금지
외부펌금지
글쓰기