数学建模社区-数学中国
标题:
求助，matlab second derivative approximation,导数逼近
[打印本页]
作者:
splelvis
时间:
2010-2-19 05:39
标题:
求助，matlab second derivative approximation,导数逼近
clear all;
8 S t7 C" h& V
nx=50;
( J) Z3 Q/ U0 g8 P4 f
ny=30;
* B! M6 p2 C9 g: M. \
hx=pi/nx;
/ F; |$ J" L0 o: Z+ ~
x=linspace(0,pi,nx+1);
( \4 L4 T4 C# X) u0 W! r4 @
y=linspace(0,pi,ny+1);
& C3 r7 x& B& a) t
x_plus_h=x+hx.*ones(1,nx+1);
9 P% `' L- P' h: j( J5 l
x_minus_h=x-hx.*ones(1,nx+1);
, b& I2 O/ o( L! C# Y
for i=1:nx+1
! `; a. ]. J5 T1 x
for j=1:ny+1
- S- Y3 T* J" Q+ c9 M
f_xx(j,i)=(f9(x_plus_h(i),y(j))-2*f9+f9(x_minus_h(i),y(j)))./(hx.^2); %原式f(x+hx,y)-2f(x,y)+f(x-hx,y)./(hx^2)
' J, X/ I7 {) D% c. o
end;
9 F% U) ~1 w6 M3 L& c* R" I+ w
end;
% K$ D. ?* K& [" I/ K
[xx,yy]=meshgrid(x,y);
: W1 J/ A/ [) o
surf(xx,yy,f_xx);
3 R0 d+ T, {0 `- x9 O
title('partial derivative w.r.t.xx')
' M6 ]$ r5 f7 k, i1 p
，
9 W# x% E4 Z- k: {* b
# c$ c& k8 K f# L+ B, y" G8 J, J, j
另开一个f9文件，
& y" j& }, P- l' v s
function l=f9(x,y)
; j# k7 Z B* Y1 Y6 Z+ E, u$ L
l=exp(sin(x-y));
* x& F* k/ I( \ J
7 f2 N7 z0 m, D
G' m4 Z% Y* ~: r3 B3 s4 F
然后出来说错误了，
, P; z0 `( k, X2 ?: u* Q
??? Input argument "x" is undefined.
6 W: H5 B1 o. F2 I& d
8 b! p @1 }6 k6 F) P* x
Error in ==> f9 at 2
+ ^/ x# ?* T, F" L$ r* b4 l) [
l=exp(sin(x-y));
8 ]1 A. P4 C' H3 X+ V
Error in ==> partial_derivxx at 11
$ K5 j B, L+ U! [. {
f_xx(j,i)=(f9(x_plus_h(i),y(j))-2*f9+f9(x_minus_h(i),y(j)))./(hx.^2);
% s8 z2 S6 m5 ^, a
: t& M* S( B) i% [* |5 d
各位要怎么改才对呢
作者:
splelvis
时间:
2010-2-19 07:35
我自己找到了问题，f9应该是 f9(x(i),y(j))
作者:
douglas8j
时间:
2010-2-19 07:35
abc~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
作者:
splelvis
时间:
2010-2-19 07:49
不好意思，我也觉得自己太粗心了
作者:
douglas8j
时间:
2010-2-19 07:59
henhaodeyipianwenzhang~~~~~~~~~~~~~~~`
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5