matlab旋转+双线性插值

Matlab中插值函数汇总和使用说明MATLAB中的插值函数为interp1，其调用格式为： yi= interp1(x,y,xi,'method')其中x，y为插值点，yi为在被插值点xi处的插值结果；x,y为向量， 'method'表示采用的插值方法，MATLAB提供的插值方法有几种： 'method'是最邻近插值， 'linear'线性插值； 'spline'三次样条插值； 'cubic'立方插值．缺省时表示线性插值注意：所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。

例如：在一天24小时内，从零点开始每间隔2小时测得的环境温度数据分别为12，9，9，10，18 ，24，28，27，25，20，18，15，13，推测中午12点（即13点）时的温度．x=0:2:24;y=[12 9 9 10 18 24 28 27 25 20 18 15 13];a=13;y1=interp1(x,y,a,'spline')结果为： 27.8725若要得到一天24小时的温度曲线，则：xi=0:1/3600:24;yi=interp1(x,y,xi, 'spline');plot(x,y,'o' ,xi,yi)命令1 interp1功能一维数据插值（表格查找）。

该命令对数据点之间计算内插值。

它找出一元函数f(x)在中间点的数值。

其中函数f(x)由所给数据决定。

x：原始数据点Y：原始数据点xi：插值点Yi：插值点格式(1)yi = interp1(x,Y,xi)返回插值向量yi，每一元素对应于参量xi，同时由向量x 与Y 的内插值决定。

参量x 指定数据Y 的点。

若Y 为一矩阵，则按Y 的每列计算。

yi是阶数为length(xi)*size(Y,2)的输出矩阵。

自己写的Matlab 旋转+双线性插值图像函数效果图：源码：clear all;I = imread('original.jpg');[Height,Width,RGB] = size(I);II = I;%当角度为0时直接输出%本程序是以左上角为坐标原点%angle_j是旋转角度，正值是按顺时针旋转，负值时按逆时针旋转angle_j = 181;%angle是弧度angle = 2*pi*angle_j/360;%将angle转成正值while(angle < 0)angle = 2 * pi + angle;end%约束在0-2π内while(angle > 2 * pi)angle = angle - 2 * pi;end%tag是判断下面的while循环有没有执行过tag = 0;while(angle > 0)%超过90度的旋转，都先旋转90度，直到角度在0°-90°之间%原理是旋转90度整数倍时，信息是不丢失的if angle >= pi/2a = pi/2;angle = angle - pi/2;elseif 0 < angle < pi/2a = angle;angle = 0;endif tag == 0tag = 1;elseI = II;[Height,Width,RGB] = size(I);%在旋转后的图像上继续旋转，从而实现大于90°的旋转end%正向变换用sina = sin(a);cosa = cos(a);%逆向变换用_m == _minussina_m = sin(-a);cosa_m = cos(-a);%旋转后图像的长度和宽度II_height = round(sina * Width + cosa * Height);II_width = round(sina * Height + cosa * Width);II = ones(II_height,II_width,3);%先转成unit8。

Matlab中插值函数汇总和使用说明MATLAB中的插值函数为interp1，其调用格式为：yi= interp1(x,y,xi,'method')其中x，y为插值点，yi为在被插值点xi处的插值结果；x,y为向量， 'method'表示采用的插值方法，MATLAB提供的插值方法有几种： 'method'是最邻近插值， 'linear'线性插值； 'spline'三次样条插值； 'cubic'立方插值．缺省时表示线性插值注意：所有的插值方法都要求x是单调的，并且xi不能够超过x的范围。

例如：在一天24小时内，从零点开始每间隔2小时测得的环境温度数据分别为12，9，9，10，18 ，24，28，27，25，20，18，15，13，推测中午12点（即13点）时的温度．x=0:2:24;y=[12 9 9 10 18 24 28 27 25 20 18 15 13];a=13;y1=interp1(x,y,a,'spline')结果为： 27.8725若要得到一天24小时的温度曲线，则：xi=0:1/3600:24;yi=interp1(x,y,xi, 'spline');plot(x,y,'o' ,xi,yi)命令1 interp1功能一维数据插值（表格查找）。

该命令对数据点之间计算内插值。

它找出一元函数f(x)在中间点的数值。

其中函数f(x)由所给数据决定。

x：原始数据点Y：原始数据点xi：插值点Yi：插值点格式(1)yi = interp1(x,Y,xi)返回插值向量yi，每一元素对应于参量xi，同时由向量x 与Y 的内插值决定。

参量x 指定数据Y 的点。

若Y 为一矩阵，则按Y 的每列计算。

yi 是阶数为length(xi)*size(Y,2)的输出矩阵。

matlab多项式插值函数
在MATLAB中，多项式插值函数可以通过使用interp1、
interp2和interp3等函数来实现。

其中，interp1用于一维插值，interp2用于二维插值，interp3用于三维插值。

这些函数可以用来在给定的数据点上进行插值，生成一个平滑的曲线或曲面。

以interp1为例，可以使用以下语法进行多项式插值：
matlab.
interp1(x, y, xi, 'method')。

其中，x和y分别是原始数据点的横坐标和纵坐标，xi是要进行插值的点的横坐标，'method'是插值方法，可以选择 'linear'（线性插值）、'nearest'（最近邻插值）、'spline'（样条插值）或 'pchip'（分段三次 Hermite 插值）等。

除了interp1之外，还可以使用polyfit和polyval函数进行多项式插值。

polyfit用于拟合多项式系数，polyval用于计算插值点的函数值。

在实际使用中，可以根据具体的数据特点和需要选择合适的插值方法，以及适当的多项式次数，来实现多项式插值函数的应用。

总之，MATLAB提供了丰富的函数和工具，可以帮助我们实现多项式插值，并且可以根据具体情况选择合适的方法和参数进行插值计算。

希望以上信息能够帮助到你。

插值方法晚上做一个曲线拟合，结果才开始用最小二乘法拟合时，拟合出来的东西太难看了！于是尝试用其他方法。

经过一番按图索骥，终于发现做曲线拟合的话，采用插值法是比较理想的方法。

尤其是样条插值，插完后线条十分光滑。

方法付后，最关键的问题是求解时要积分，放这里想要的时候就可以直接过来拿，不用死去搜索啦。

呵呵插值方法的Matlab实现一维数据插值MATLAB中用函数interp1来拟合一维数据，语法是YI = INTERP1(X,Y,XI,方法)其中（X，Y）是已给的数据点，XI 是插值点，其中方法主要有'linear' -线性插值，默认'pchip' -逐段三次Hermite插值'spline' -逐段三次样条函数插值其中最后一种插值的曲线比较平滑例：x=0:.12:1; x1=0:.02:1;%（其中x=0:.12:1表示显示的插值点，x1=0:.02:1表示插值的步长）y=(x.^2-3*x+5).*exp(-5*x).*sin(x);plot(x,y,'o'); hold on;y1=interp1(x,y,x1,'spline');plot(x1,y1,':')如果要根据样本点求函数的定积分，而函数又是比较光滑的，则可以用样条函数进行插值后再积分，在MATLAB中可以编写如下程序：function y=quadspln(x0,y0,a,b)f=inline('interp1(x0,y0,x,''spline'')','x','x0','y0');y=quadl(f,a,b,1e-8,[],x0,y0);现求sin(x)在区间[0，pi]上的定积分，只取5点x0=[0,0.4,1,2,pi];y0=sin(x0);I=quadspln(x0,y0,0,pi)结果得到的值为2.01905，精确值为2求一段matlab插值程序悬赏分：20 - 解决时间：2009-12-26 19:57已知5个数据点：x=[0.25 0.5 0.75 1] y=[0 0.3104 0.6177 0.7886 1] ，求一段matlab插值程序，求过这5个数据点的插值多项式，并在x-y坐标中画出y=f（x）图形，并且求出f （x）与x轴围成图形的面积（积分），不胜感激！使用Lagrange 插值多项式的方法：首先把下面的代码复制到M文件中，保存成lagranfunction [C,L]=lagran(X,Y)% input - X is a vector that contains a list of abscissas% - Y is a vector that contains a list of ordinates% output - C is a matrix that contains the coefficients of the lagrange interpolatory polynomial%- L is a matrix that contains the lagrange coefficients polynomialw=length(X);n=w-1;L=zeros(w,w);for k=1:n+1V=1;for j=1:n+1if k~=jV=conv(V,poly(X(j)))/(X(k)-X(j));endendL(k,:)=V;endC=Y*L;然后在命令窗口中输入以下内容：x=[0 0.25 0.5 0.75 1];y=[0 0.3104 0.6177 0.7886 1];lagran(x,y)ans =3.3088 -6.3851 3.3164 0.7599 0得到的数据就是多项式各项的系数，注意最后一个是常数项，即x^0，所以表达式为:f=3.3088*x.^4-6.3851*x.^3+3.3164*x.^2 +0.7599*x求面积就是积分求解>> f=@(x)3.3088*x.^4-6.3851*x.^3+3.3164*x.^2 +0.7599*x;>> quad(f,0,1)ans =0.5509这些点肯定是通过这个多项式的！MATLAB插值与拟合§1曲线拟合实例：温度曲线问题气象部门观测到一天某些时刻的温度变化数据为：试描绘出温度变化曲线。

matlab 插值法MATLAB 插值法是数据处理和信号处理中常用的一种算法。

在数据采集或数据处理中，通常会遇到数据缺失或者采样点不足的情况，这时候就需要用到插值法来对数据进行补充或者重构。

插值法的基本思想是，给定一些离散的数据点，通过一种数学方法，构造出一个连续的函数，使得在已知数据点处，该函数与原数据点一致。

常见的插值方法有线性插值、多项式插值、样条插值等。

线性插值法是最简单的一种插值方法。

在采样点之间的区域内，采用一次多项式函数去逼近该区域内的某个未知函数。

其公式如下所示：f(x) = f(x0)(1 - t) + f(x1)t其中，x0 和 x1 是相邻两个采样点，t 是一个权重系数，表示该点在两个采样点之间的位置。

多项式插值法是用一个 n 次多项式函数逼近原函数 f(x)。

在采样点处，两个函数的取值相同，同时也能保证一定的光滑性。

其公式如下所示：f(x) = a0 + a1x + a2x^2 + ... + anxnS''(x) = M0(x - x0) + N0, x0 ≤ x ≤ x1其中，M 和 N 是未知的系数，通过计算两个相邻区间中的连续性和光滑性来解出系数。

除了以上三种插值方法，还有其他的插值算法，例如离散傅里叶插值法、拉格朗日插值法等。

总之，MATLAB 中的插值函数为 interp1，它的语法格式如下：yi = interp1(x, y, xi, method)其中，x 和 y 为已知函数的取值点，xi 为要进行插值的点的位置，method 是采用的插值方式。

例如，method = 'linear' 表示采用线性插值法。

MATLAB 中还提供了很多其他的 method 选项，用户可以根据实际情况选择适合的方法。

MATLAB 插值算法在信号处理和图像处理中广泛应用，例如，图像的放大缩小、色彩调整、去噪等都可以用插值算法实现。

因此，掌握 MATLAB 插值算法可以帮助我们更好地进行数据处理和信号处理。

% 本程序实现双线性插值% 对于取整% there is a decimal a, the value of round and ceil is the nearest integer >a;% on the contrast, floor and fix is the nearest integer <a;% round = ceil ；floor = fix% 双线性插值：% 输出像素值是它在输入图像中2*2邻域采样点的平均值，% 它根据某像素周围4个像素的灰度值在水平和垂直两个方法上对其插值。

为了方便理解，先考虑一维情况下的线性插值：% 对于一个数列c，% 假设c[a]到c[a+1]之间是线性变化的，那么对于浮点数% %-------------------第一种表达% x( a <= x < a+1)% c(x) = c[a+1]*( x - a) + c[a]*( 1 + a - x);% c(x) = c[a] + [ (c[a+1]-c[a])/(b-a)]*( x - a);% b = a+1;% c(x) = c[a] + [ c[a+1] - c[a]]*( x - a);% 把这种插值方式扩展到二维情况：对于一个二维数组c，% 我们假设对于任意一个浮点数i,c(a,i)到c(a+1,i)之间是线性变化的,c(i,b)到c(i,b+1)之间也是线性变化的(a,b都是整数) ，% 那么对于浮点数的坐标(x,y)满足(a <= x < a+1, b <= y < b+1)，可以先分别求出c(x,b)和c(x,b+1):% c(x,b) = c[a+1][b]*( x - a) + c[a][b]*( 1 + a - x);% c(x,b+1) = c[a+1][b+1]*( x - a) + c[a][b+1]*( 1 + a - x);% 好，现在已经知道c(x,b)和c(x,b+1)了，而根据假设c(x,b)到c(x,b+1)也是线性变化的，所以:% c(x,y) = c(x,b+1)*( y - b) + c(x,b)*( 1 + b - y)% %------------第二种表达% c(x) = c[a+1]*( x - a) + c[a]*( 1 + a - x);% = c[a] + c[a+1]*( x - a) - c[a]( x - a);% c(x) = c[a] + ( c[a+1] - c[a])*( x - a);% 对于二维情况：% c(x,b) = c[a,b] + ( c[a+1,b] - c[a,b])*( x - a);% c(x,b+1) = c[a,b+1] + ( c[a+1,b+1] - c[a,b+1])*( x - a);% c(x,y) = c(x,b) + ( c(x,b+1) - c(x,b))*( y - b);%---------------------------------------------------------------------clc;clear all;Image = imread('1218.bmp');grayImage = rgb2gray(Image);figure,imshow(Image);figure,imshow(grayImage);% grayImage = [1 2 3 4 5;1 2 3 4 5;1 2 3 4 5;1 2 3 4 5;1 2 3 4 5];% grayImage = [1 2 3;1 2 3;1 2 3];rotation = [0.5 0.5;-0.5 0.5];[rows cols] = size(grayImage);grayImage = double(grayImage);for r = 1:rowsfor c = 1:colstemp = rotation*[r-rows/2;c-cols/2]+[rows/2;cols/2]; % 旋转% temp = [rows/2;cols/2];if( temp(1)>1 && temp(2)>1 && temp(1)<rows && temp(2)<cols) % 双线性插值x = floor(temp(1));y = floor(temp(2));% 插值公式1% v1 = grayImage( x+1,y)*( temp(1)-x) + grayImage( x,y)*( 1+x-temp(1));% v2 = grayImage( x+1,y+1)*( temp(1)-x) + grayImage( x,y+1)*( 1+x-temp(1)); % v = v2*( temp(2)-y) + v1*( 1+y-temp(2));% 插值公式2v1 = grayImage(x,y) + ( grayImage(x+1,y) - grayImage(x,y))*( temp(1) - x);v2 = grayImage(x,y+1) + ( grayImage(x+1,y+1) - grayImage(x,y+1))*( temp(1) - x);v = v1 + ( v2 - v1)*( temp(2) - y);grayImageN(r,c) = round(v);endendendgrayImage = uint8(grayImageN);figure,imshow(grayImage);本文来自: 高校自动化网() 详细出处参考(转载请保留本链接)：/html/matlab/15410_2.html。

实验三图像的几何操作一、图像的基本运算1、加法运算图像相加一般用于对同一场景的多幅图像求平均效果，以便有效地降低具有叠加性质的随机噪声。

直接采集的图像品质一般都较好，不需要进行加法运算处理，但是对于那些经过长距离模拟通信方式传送的图像（如卫星图像），这种处理是必不可少的。

在MA TLAB7.0中，如果要进行两幅图像的加法，或者给一幅图像加上一个常数，可以调用imadd函数来实现。

imadd函数将某一幅输入图像的每一个像素值与另一幅图像相应的像素值相加，返回相应的像素值之和作为输入图像。

imadd函数的调用格式可参考图像处理的工具箱。

下面的程序可将图1两幅图像叠加在一起，叠加效果如图2所示。

imshow('rice.png');imshow('cameraman.tif');I=imread('rice.png');J=imread('cameraman.tif');K=imadd(I,J,'uint16');%大小必须一样imshow(K,[])图1 图2给图像的每一个像素加上一个常数可以使图像的亮度增加。

例如以下程序示例的处理效果如图3所示。

I=imread('rice.png');J=imadd(I,50);subplot(1,2,1),imshow(I);subplot(1,2,2),imshow(J);图32、减法运算图像减法也称为差分方法，是一种常用于检测图像变化及运动物体的图像处理方法。

图像减法可以作为许多图像处理过程的准备步骤。

例如，可以使用图像减法来检测一系列相同场景图像的差异。

图像减法与阈值化处理的综合使用通常是建立机器视觉系统最有效的方法之一。

当然，在利用图像减法处理图像时，往往需要考虑背景的更新机制，尽量补偿因天气、光线等因素对图像显示效果造成的影响。

在MATLAB7.0中，使用imsubtract函数可以将一幅图像从另一幅图像中减去，或者从一幅图像中减去一个常数。

matlab插值函数Matlab（矩阵实验室）是一种高级的数学软件，它可以帮助人们解决复杂的数学和工程问题。

其中最重要的功能之一就是插值函数。

插值函数是一种在非等间距的数据点之间通过计算曲线或其他曲面来估算未知值的技术。

本文将就插值函数在MATLAB中的发挥，以及MATLAB提供的插值函数进行详细的介绍。

一、关于MATLAB插值函数的简介MATLAB的插值函数是一种估算未知值的方法，它可以帮助人们在非等间距的数据点之间通过计算曲线或其他曲面来估算未知值。

使用插值函数，可以从已知数据中推断未知数据。

MATLAB插值函数可以用于曲线拟合、寻找特定函数极值、以及求解线性和非线性方程组。

一般来说，使用插值函数进行重新排列或插值时，数据点之间的距离应尽量小，而不是间歇的大距离。

二、MATLAB提供的插值函数MATLAB拥有多种内置的插值函数，以下是MATLAB中最常用的几种插值函数：1.性插值：Linear interpolation，它将所求点放在两个已知点之间，并使用两个已知点的线性函数值来拟合它；2.式插值：多项式插值是使用一组已知点拟合一个多项式来估算未知点的最小二乘法插值法；3.条插值：样条插值是利用一些已知点来拟合出一个连续可微的样条函数来进行插值处理的；4.数插值：指数插值是根据一些已知的指数函数拟合出曲线来做插值处理的。

三、MATLAB插值函数的应用MATLAB插值函数的应用非常广泛，它可以用来解决和处理复杂的数学和工程问题。

例如，可以使用MATLAB插值函数来拟合数据；对解决非线性方程组有很大帮助；可以将数据绘制出来，以便于观察、比较、识别出特殊的性质；还可以用来估算未知函数值；最后还可以根据求解结果来求解极限问题，等等。

四、总结本文介绍了MATLAB插值函数的基本概念，以及MATLAB提供的几种常见的插值函数，包括线性插值、多项式插值、样条插值和指数插值。

这些插值函数的应用及其广泛，可以用来拟合复杂的数据，以及帮助解决一些复杂的数学和工程问题。