회원가입 ID/PW 찾기
AA

문제는 

-xu"-u'+u*x^2=cosx

u(0)=u(1)=0

0<x<1

일 때 u 의 해에 대해 plotting하는 건데요...

local 행렬을 구해서 u에 대한 global 행렬로 assembling해서 푸는 형식입니다.

저는 이 문제에서 전체 term을 x로 나눠서 

-u"-(1/x)*u'+x*u=cosx/x

이렇게 해서 boundary condition을 더 간단하게 만들어서 푸는 중인데요

도저히 그 중의 한 term인 (cosx/x)((-x+b)/(b-a))부분과 (cosx/x)((x-a)/(b-a))부분의 적분결과를 설정할 수 가 없어서요..

cosx/x를 적분하면 element로 표현하기 어려워서 function으로 해보았는데,,,아직 초보라 제가 function을 제대로 설정하였는지도 모르겠어요..

여기서 a와 b는 0부터 1까지를 n등분한 것 중의 한 개의 element part이구요

일단 첫번째는

function lF1=int(N1,N2);

lF1=(cos(x)/x)*((-x+N2)/(N2-N1)); 이걸 엠파일로 만들어서


lF(1,1)=  int(sym('int'),N1,N2); 이렇게 실행을 시켰는데.......

실행은 되요...맞는지는 모르겠지만..ㅠㅠ 이런식으로 두번째 적분도 했구요...


그다음엔..정확한 정답하고 제가 코딩한 거하고 비교해야하는데 

여기서도 좀 잘못 되지 않았을까 싶어요.


x_1=eCoord;


% compare the solution with that obtained using Matlab built-in functions.

solinit = bvpinit(linspace(0.00000001,1,100), [0 1]);

sol = bvp4c(@problem_,@bc,solinit);


alpha=linspace(0.00000001,1,100);

beta=deval(sol,alpha);


title=sprintf('comparison of FEM solution=%3d',Nelem);

figure('Name',title);

plot(x_1,sU, alpha, beta(1,:),'-ro');

xlabel('x')

ylabel('y')    


여기서 


problem_.m은 

function dydx = problem_(x,y)


dydx = [y(2) ; x*y(1)-(1/x)*y(2)-cos(x)/x];


bc.m은 

function res = bc(ya,yb)

res = [ya(1)-0   yb(1)-0];로 했구요...

뭐가 문제일까요....

ㅠㅠ 그래프는 점점 산으로 가네요..도와주세요

댓글 0

하드웨어 설계 및 개발에 대하여 개발자들이 자유롭게 토론하는 공간입니다.
- Q&A, 자유주재 토론, 관련 정보 공유
- 분야 : 마이크로프로세서 응용, 전기/전자(아날로그/디지털) 회로 설계, C/C++ 프로그래밍, 펌웨어,
         PCB Artwork, 트러블슈팅 등 하드웨어 설계에 관한 전반인 내용
※ 게시글에 맞는 분류를 선택하여 글을 작성해 주시면 쾌적한 사이트 운영에 많은 도움이 됩니다.
※ 하드웨어 인사이트는 회원들간의 거래정보를 게재할 뿐이지, 그 어떤 책임과 의무도 가지지 않습니다.

search
번호 분류 제목 글쓴이 조회 수 날짜
3 머신러닝, AI & 알고리즘 HOT오목 게임 알고리즘3 새로운하늘 3457 2010.03.29
2 머신러닝, AI & 알고리즘 HOT지하철 최단거리 알고리즘은?2 지워나 2859 2008.06.20
1 머신러닝, AI & 알고리즘 HOT순차 프로그램에 대한 간단한 알고리즘과 재귀 알고리즘의 속도 테스트 - by 아크마7 아크마 3056 2007.08.08
  • 인생은 선을 실행하기 위하여 만들어졌다.
    - 칸트
  • * 납포인트 정보 *
  • 글 작성 : 3
  • 댓글 작성 : 1
  • 내 글이 추천받음 : 1
저작권법에 위배되는 콘텐츠는 등록 불가하며, 저작물에 대한 권리는 저작자에게 있습니다.
Copyright 2006-2021 © hardwareis.com, All rights reserved.