matlab csapi代码,【matlab】樣條函數:csapi、spap2和augknt的使用
可以調用樣條工具箱(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的使用相关推荐
- 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); %畫訓練誤差 ...
- imageset matlab,如何以imageSet或imageDataStore的形式向MATLAB中的BagOfFeatures()函數提供輸入?...
我想使用MATLAB的bagOfFeatures()函數.但它需要以imageSet或imageDataStore的形式輸入.我想運行的代碼如下:如何以imageSet或imageDataStore的 ...
- matlab工具函数newp,常用的 MATLAB 神經網絡工具箱函數
列表中所列出的函數適用於 MATLAB5.3 以上版本,若需要進一步的了解,請參閱 MATLAB 的幫助文檔. 1. 網絡創建函數 newp 創建感知器網絡 newlind 設計一線性層 ne ...
- matlab lbm 代码,Matlab实现格子玻尔兹曼方法(Lattice Boltzmann Method,LBM)模拟
%LBM的matlab代码 %Matlab实现格子玻尔兹曼方法(Lattice Boltzmann Method,LBM)模拟 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % c ...
- matlab 五子棋代码,matlab编程(五子棋)
带有详细中文注释的经典Matlab五子棋程序代码 break end if y+i>10 break; end if board(x+i,y+i)~=flag 判断在此过程中,棋子与刚下的棋的颜 ...
- matlab bwlable代码,matlab的bwlabel
精品文档 clear; clc; matlab 标注连通域 f=imread('c:\\1.jpg'); gray_level=graythresh(f); f=im2bw(f,gray_level) ...
- matlab振动函数代码,Matlab振动程序-代码作业.docx
一.课题任务要求 随着机械工业不断向自动化.高精度.智能化等方向的发展,在机械设备运行及生产过程中进行参量测试.分析与诊断等处理过程已成为必要环节,许多信号处理方法如时域统计分析.相关分析.相干分析. ...
- matlab 实验代码,matlab实验2代码
1.画出对数和指数函数曲线,并分别加上标题.轴标记和曲线说明(这里可采用多种方法来标注曲线) 解: x=0.01:0.1:10; y1=log10(x); y2=exp(x); figure(1) s ...
- matlab+awgn和wgn,噪聲強度(噪聲功率) 噪聲方差到底有什么關系? matlab中的awgn函數...
以matlab中awgn函數為例說明: 在matlab中無論是wgn還是awgn函數,實質都是由randn函數產生的噪聲.即:wgn函數中調用了randn函數,而awgn函數中調用了wgn函數. 根據 ...
最新文章
- ApplicationStartedEvent 与 ContextStartedEvent 有区别吗?
- 《Windows Server 2012 Hyper-V虚拟化管理实践》——3.2 Hyper-V主机日常管理
- ELK学习6_Kafka-Logstash-Elasticsearch数据流操作
- 取sql数据乱码_不基于备份和表,生产系统数据误删就能完全恢复
- OpenMap教程3 –使用MapHandler构建基本的地图应用程序–第2部分
- ScrumChina 2008 上海之行
- RxJava 的基本使用
- vim文件时自动添加作者、时间、版权等信息
- 190407每日一句
- paip.android 读取docx总结
- java工程师什么城市就业_热门城市的Java薪资情况
- hdoj1290切球形蛋糕(递推和划分问题)
- 中南大学c语言课程设计报告,中南大学 C++课程设计报告.doc
- 七彩安卓影视APP源码独立解析接口
- 京东、搜狗“带狗”都好好的,为何司机对快狗打车不满意?
- 多线程学习(九)-可阻塞的队列
- linux ora-39002,ORA-39002:invalidoperation
- workbench设置单元坐标系_ansys workbench中新建坐标系的问题,求高人指点。
- 在c51语言的循环语句中 用作循环结束,在C51语言的循环语句中,用作循环结束条件判断的表达式为()...
- 百度地图经纬度和像素坐标互转
热门文章
- java jvisualvm linux,从Linux JDK中发出jvisualvm时出现乱码
- callback函数_Nodejs 源码解析 util.promisify 如何将 Callback 转为 Promise
- java select 不是date,iPhonedateselect器,而不是键盘?
- mysql json类型存储list_list存储json字符串
- HP-UX Oracle Ioctl ASYNC_CONFIG error=1
- 在linux下运行qiime,科学网—QIIME2在Linux环境下使用遇到的问题 - 肖斌的博文
- 自定义Hibernate Validator规则注解
- Spring DI(依赖注入)构造器注入篇
- springMVC后台验证信息不显示的原因初探
- Apache httpd 安装