数学建模社区-数学中国
标题:
Bernstein基构建函数数据MATLAB代码
[打印本页]
作者:
建不了的模。
时间:
2015-1-20 13:50
标题:
Bernstein基构建函数数据MATLAB代码
%================Bernstein基构建函数数据(二维)==========================
4 m7 {2 S: h6 V l
%
! K: I( {2 h6 o, @8 s" l& D* ?
%see also
http://www.matlabsky.com
0 o# S0 x! W$ J. i0 s7 H8 v# P4 o; o
%
' m" R o- a. B+ O
%===得到Bernsein基的次数================
" R& t3 h8 I! { Y) W
clear;
5 a" e: ^9 p( B! i" @5 a8 u
X=load('data.txt');
" V" |6 t& }7 @ G1 K
%画出原图形的
" E _1 u- Q" F& {9 G
plot(X(:,1),X(:,2));
1 p9 f9 a: X q, \( i+ ?
title('原始数据');
* d. _, I; s/ j2 w; v3 P* O; N
%数据参数化,将X(:,1)变换为[0,1]间数据,
' d# ^+ w3 {. n* y9 H8 n, Y5 z
minX=min(X(:,1));
" Y. T) F+ @/ u' E& E- G" ^& @. A0 q
maxX=max(X(:,1));
2 J5 {- Y' c: D
paraX=0;
- q o! q1 `8 V5 N6 {# |/ u, X- s
for i=1:length(X)
9 E1 Z; f/ r5 M3 k8 `# {$ Y |4 g
paraX(i)=(X(i,1)-minX)/(maxX-minX);
, i' g# c" C' W. `0 J3 v: t% m
end
5 h/ t7 K6 P: P
m=input('请输入Bernsein基的次数m(即m+1阶):');
) ] k% a: S9 p' x& n7 g& A
b=OLS(X,paraX,m);
7 v0 D E( w% C' r' C. J& Q) |, V4 @- z
residu=residual(X,m,paraX,b);
' S9 B, P5 h( K w' Y+ u
%循环得读取命令
3 f |% \0 ~% Z( P& b0 M9 e2 k
while 1
0 p% B1 w* |# W `
fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
& V! e0 y5 H- W# B$ _; v
fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
! I* S" V1 h. ~6 s6 ~
fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
4 H. e" I* Q8 ]# V
fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
. ~* A/ O- {% v8 H- {+ G
cmd=input('请按照以上提示选择操作:');
( a' g1 P% o: T4 F
! A- v( O+ b, q) d# N
switch cmd
4 W: |( `; [2 U5 {2 a/ j- ?
case 1
5 f4 m1 C8 H: ^- ^' Z* o6 E* X/ }! k4 U4 S
; A9 U- A* k7 O- N" P9 `
6 B8 k/ O; P5 w& N- H% h- v
# |, T' I/ d! [& o4 h
作者:
宇仲
时间:
2015-1-20 14:28
楼主辛苦了，继续加油啊！
& S" B3 m$ K" [$ d x( k0 b9 z3 E
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5