可以調用樣條工具箱(Spline ToolBox)

csapi

clear;clc

z0=double(rgb2gray(imread('1.jpg'))); %讀取數據

x0=1:size(z0,1);y0=1:size(z0,2);

sp=csapi({x0 y0},z0);%得到每片三次樣條曲面的系數

fnplt(sp); %可得到擬合曲面的效果圖

vals = fnval(sp,{x0 y0});%計算擬合后的數據

---------------------------------------------------------------------------------------------------------------

spap2

clc; clear; close all;

z=[0.2 0.24 0.25 0.26 0.25 0.25 0.25 0.26 0.26 0.29 0.25 0.29;

0.27 0.31 0.3 0.3 0.26 0.28 0.29 0.26 0.26 0.26 0.26 0.29;

0.41 0.41 0.37 0.37 0.38 0.35 0.34 0.35 0.35 0.34 0.35 0.35;

0.41 0.42 0.42 0.41 0.4 0.39 0.39 0.38 0.36 0.36 0.36 0.36;

0.3 0.36 0.4 0.43 0.45 0.45 0.51 0.42 0.4 0.37 0.37 0.37];

% z=double(rgb2gray(imread('1.jpg')));

x=1:size(z,2);%長為12

y=1:size(z,1);%寬為5

[x0, y0] = meshgrid(x, y);%產生采樣網格

surf(x,y,z);%畫出原始數據的三維圖

title('Original data Plot');

xlabel('X'), ylabel('Y'), zlabel('Z')

% colormap; colorbar;

% axis([0 15 0 6 0.15 0.55]);

% hold on;

% plot3(x0, y0, z, 'p');

figure;

ky = 3; knotsy = augknt([1,2.5,13],ky);%設置節點序列1 1 1 2.5 13 13 13

%augknt(knots,k)意義是產生一個逐增的節點序列,並且頭和尾節點重復為k個,

%如augknt([1 2 3 3 3 4 4 4 4],2)和augknt([1 2 3 3 3 4],2)得到結果一樣

%1 1 2 3 3 3 4 4

sp = spap2(knotsy,ky,y,z');%根據最小二乘准則的B樣條擬合

yy = 1:size(z,1); vals = fnval(sp,yy);%計算擬合的數據

fnplt(sp); %畫出擬合效果圖

matlab csapi代码,【matlab】樣條函數:csapi、spap2和augknt的使用相关推荐

  1. newff matlab使用例子,matlab神經網絡newff函數的用法

    設[P,T]是訓練樣本,[X,Y]是測試樣本: net=newrb(P,T,err_goal,spread); %建立網絡 q=sim(net,p); e=q-T; plot(p,q); %畫訓練誤差 ...

  2. imageset matlab,如何以imageSet或imageDataStore的形式向MATLAB中的BagOfFeatures()函數提供輸入?...

    我想使用MATLAB的bagOfFeatures()函數.但它需要以imageSet或imageDataStore的形式輸入.我想運行的代碼如下:如何以imageSet或imageDataStore的 ...

  3. matlab工具函数newp,常用的 MATLAB 神經網絡工具箱函數

    列表中所列出的函數適用於 MATLAB5.3 以上版本,若需要進一步的了解,請參閱 MATLAB 的幫助文檔. 1. 網絡創建函數 newp   創建感知器網絡 newlind   設計一線性層 ne ...

  4. matlab lbm 代码,Matlab实现格子玻尔兹曼方法(Lattice Boltzmann Method,LBM)模拟

    %LBM的matlab代码 %Matlab实现格子玻尔兹曼方法(Lattice Boltzmann Method,LBM)模拟 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % c ...

  5. matlab 五子棋代码,matlab编程(五子棋)

    带有详细中文注释的经典Matlab五子棋程序代码 break end if y+i>10 break; end if board(x+i,y+i)~=flag 判断在此过程中,棋子与刚下的棋的颜 ...

  6. matlab bwlable代码,matlab的bwlabel

    精品文档 clear; clc; matlab 标注连通域 f=imread('c:\\1.jpg'); gray_level=graythresh(f); f=im2bw(f,gray_level) ...

  7. matlab振动函数代码,Matlab振动程序-代码作业.docx

    一.课题任务要求 随着机械工业不断向自动化.高精度.智能化等方向的发展,在机械设备运行及生产过程中进行参量测试.分析与诊断等处理过程已成为必要环节,许多信号处理方法如时域统计分析.相关分析.相干分析. ...

  8. matlab 实验代码,matlab实验2代码

    1.画出对数和指数函数曲线,并分别加上标题.轴标记和曲线说明(这里可采用多种方法来标注曲线) 解: x=0.01:0.1:10; y1=log10(x); y2=exp(x); figure(1) s ...

  9. matlab+awgn和wgn,噪聲強度(噪聲功率) 噪聲方差到底有什么關系? matlab中的awgn函數...

    以matlab中awgn函數為例說明: 在matlab中無論是wgn還是awgn函數,實質都是由randn函數產生的噪聲.即:wgn函數中調用了randn函數,而awgn函數中調用了wgn函數. 根據 ...

最新文章

  1. ApplicationStartedEvent 与 ContextStartedEvent 有区别吗?
  2. 《Windows Server 2012 Hyper-V虚拟化管理实践》——3.2 Hyper-V主机日常管理
  3. ELK学习6_Kafka-Logstash-Elasticsearch数据流操作
  4. 取sql数据乱码_不基于备份和表,生产系统数据误删就能完全恢复
  5. OpenMap教程3 –使用MapHandler构建基本的地图应用程序–第2部分
  6. ScrumChina 2008 上海之行
  7. RxJava 的基本使用
  8. vim文件时自动添加作者、时间、版权等信息
  9. 190407每日一句
  10. paip.android 读取docx总结
  11. java工程师什么城市就业_热门城市的Java薪资情况
  12. hdoj1290切球形蛋糕(递推和划分问题)
  13. 中南大学c语言课程设计报告,中南大学 C++课程设计报告.doc
  14. 七彩安卓影视APP源码独立解析接口
  15. 京东、搜狗“带狗”都好好的,为何司机对快狗打车不满意?
  16. 多线程学习(九)-可阻塞的队列
  17. linux ora-39002,ORA-39002:invalidoperation
  18. workbench设置单元坐标系_ansys workbench中新建坐标系的问题,求高人指点。
  19. 在c51语言的循环语句中 用作循环结束,在C51语言的循环语句中,用作循环结束条件判断的表达式为()...
  20. 百度地图经纬度和像素坐标互转

热门文章

  1. java jvisualvm linux,从Linux JDK中发出jvisualvm时出现乱码
  2. callback函数_Nodejs 源码解析 util.promisify 如何将 Callback 转为 Promise
  3. java select 不是date,iPhonedateselect器,而不是键盘?
  4. mysql json类型存储list_list存储json字符串
  5. HP-UX Oracle Ioctl ASYNC_CONFIG error=1
  6. 在linux下运行qiime,科学网—QIIME2在Linux环境下使用遇到的问题 - 肖斌的博文
  7. 自定义Hibernate Validator规则注解
  8. Spring DI(依赖注入)构造器注入篇
  9. springMVC后台验证信息不显示的原因初探
  10. Apache httpd 安装