数学建模社区-数学中国
标题:
求助，matlab second derivative approximation,导数逼近
[打印本页]
作者:
splelvis
时间:
2010-2-19 05:39
标题:
求助，matlab second derivative approximation,导数逼近
clear all;
# [# c, d# |: q3 Y6 ~4 W
nx=50;
. B9 M/ t# l# C) y
ny=30;
: `5 F2 e2 R6 k
hx=pi/nx;
$ H/ O+ D7 u4 A% A
x=linspace(0,pi,nx+1);
. d5 f y; {" C% m9 g
y=linspace(0,pi,ny+1);
5 \' z% O, g& _
x_plus_h=x+hx.*ones(1,nx+1);
4 C/ z3 m$ G5 [" O1 C* |
x_minus_h=x-hx.*ones(1,nx+1);
' C# @8 v, X3 x% Z
for i=1:nx+1
6 E- e& \: k5 Z
for j=1:ny+1
: Q) ^- a' f9 p8 i; l: Q
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)
8 G$ S0 N5 z. C3 Q5 `6 B! b V
end;
/ f1 [# p( T; X, C; I! A
end;
8 |& X- B1 P8 y' z% b& o
[xx,yy]=meshgrid(x,y);
+ r) }) a! n% w2 n4 ?* {
surf(xx,yy,f_xx);
& D# t" B% ]) T0 W7 q2 O
title('partial derivative w.r.t.xx')
" h% A* U) q9 h" ]$ N
，
; ~2 g1 B! Z* W* B3 S0 @
, @6 l Z6 B. @4 o
另开一个f9文件，
}: u8 C5 v8 N1 P% {0 y. n# o/ `% l
function l=f9(x,y)
8 H( W) h! E& |! K
l=exp(sin(x-y));
6 q* [: ]) B' y& e
7 C2 J4 o% o" B. g2 e& W4 s+ W
/ n1 M) X% J. ?7 H* S' r: t! ?4 y& d
然后出来说错误了，
* w5 [7 F$ u5 T! E7 ^: J, ]
??? Input argument "x" is undefined.
3 X+ `$ S* }& [+ r6 ?/ e1 U7 K- y" z
* m/ P' c7 j& N( ~
Error in ==> f9 at 2
% K& `. a. k; \( c. B$ C- y
l=exp(sin(x-y));
& W) |8 K) A% W
Error in ==> partial_derivxx at 11
" l* n5 w/ E4 ^& u; ?; a" G) s' M
f_xx(j,i)=(f9(x_plus_h(i),y(j))-2*f9+f9(x_minus_h(i),y(j)))./(hx.^2);
, A) |! @ `3 P' K/ M
, S) R" ?. ?4 F% I. \3 F& H
各位要怎么改才对呢
作者:
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