下面是关于C++ Sobel算子的程序代码 相关知识了解 语句看不懂 希望大神能帮我注释下 方便我理解 谢谢

//Sobel算子

BOOL SySobel(HDIB hDib)

{

if(hDib==NULL)return FALSE;

LPBYTE lpBits=FindDIBBits((LPBYTE)GlobalLock(hDib));

BITMAPINFO bmi = (BITMAPINFO )GlobalLock(hDib);

LPBITMAPINFOHEADER lpbi = (LPBITMAPINFOHEADER)bmi;

int nx,ny;

int height,width;

// get color number

WORD wNumColors = DIBNumColors((LPBYTE)bmi);

height=bmi->bmiHeaderbiHeight;

width=bmi->bmiHeaderbiWidth;

double gray;

double grayl;

double grayr;

BYTE logImg;//二维缓冲区

logImg=new BYTE [height];

for(ny=0;ny<height;ny++)

logImg[ny]=new BYTE [width];

BYTE logImg1;//二维缓冲区

logImg1=new BYTE [height];

for(ny=0;ny<height;ny++)

logImg1[ny]=new BYTE [width];

BYTE logImg2;//二维缓冲区

logImg2=new BYTE [height];

for(ny=0;ny<height;ny++)

logImg2[ny]=new BYTE [width];

lpBits=(LPBYTE)lpbi+lpbi->biSize+(int)wNumColorssizeof(RGBQUAD);//获取数据矩阵指针

int nDelta = WIDTHBYTES(lpbi->biBitCountlpbi->biWidth) - lpbi->biWidthlpbi->biBitCount/8;

for ( ny=0; ny<height; ny++)

{

for ( nx=0; nx<width; nx++)

{

gray=(lpBits);

logImg[ny][nx]=(BYTE)gray;

logImg1[ny][nx]=(BYTE)gray;

logImg2[ny][nx]=(BYTE)gray;

lpBits++;

}

lpBits += nDelta;

}

for( ny=1; ny<(height-1); ny++)

{

for( nx=1; nx<(width-1); nx++)

{

//gray=(BYTE)sqrt(graylgrayl+grayrgrayr);

grayr=fabs(grayr);

grayl=fabs(grayl);

if (grayr>255)

logImg1[ny][nx]=255;

else

logImg1[ny][nx]=(BYTE)grayr;

if (grayl>255)

logImg2[ny][nx]=255;

else

logImg2[ny][nx]=(BYTE)grayr;

}

}lpBits=(LPBYTE)lpbi+lpbi->biSize+(int)wNumColorssizeof(RGBQUAD);//获取数据矩阵指针

for( ny=1; ny<(height-1); ny++)

{

for( nx=1; nx<(width-1); nx++)

{

if(logImg2[ny][nx]>logImg1[ny][nx])

{

logImg1[ny][nx]=logImg2[ny][nx];

}

}

}

for ( ny=0; ny<lpbi->biHeight; ny++)

{

for ( nx=0; nx<lpbi->biWidth; nx++)

{

gray=logImg1[ny][nx];

(lpBits)=(BYTE)gray;

lpBits++;

}

lpBits += nDelta;

}

for(ny=0;ny<height;ny++)

delete [] logImg[ny];

delete [] logImg;

for(ny=0;ny<height;ny++)

delete [] logImg1[ny];

delete [] logImg1;

for(ny=0;ny<height;ny++)

delete [] logImg2[ny];

delete [] logImg2;

GlobalUnlock(hDib);

//WaitCursorEnd();

return TRUE;

}

位运算

在很多系统程序中常要求在位(bit)一级进行运算或处理。C语言提供了位运算的功能， 这使得C语言也能像汇编语言一样用来编写系统程序。

一、位运算符C语言提供了六种位运算符：

& 按位与

| 按位或

^ 按位异或

~ 取反

<< 左移

>> 右移

1 按位与运算 按位与运算符"&"是双目运算符。其功能是参与运算的两数各对应的二进位相与。只有对应的两个二进位均为1时，结果位才为1 ，否则为0。参与运算的数以补码方式出现。

例如：9&5可写算式如下： 00001001 (9的二进制补码)&00000101 (5的二进制补码) 00000001 (1的二进制补码)可见9&5=1。

按位与运算通常用来对某些位清0或保留某些位。例如把a 的高八位清 0 ， 保留低八位， 可作 a&255 运算 ( 255 的二进制数为0000000011111111)。

main(){

int a=9,b=5,c;

c=a&b;

printf("a=%d\nb=%d\nc=%d\n",a,b,c);

}

2 按位或运算 按位或运算符“|”是双目运算符。其功能是参与运算的两数各对应的二进位相或。只要对应的二个二进位有一个为1时，结果位就为1。参与运算的两个数均以补码出现。

例如：9|5可写算式如下： 00001001|00000101

00001101 (十进制为13)可见9|5=13

main(){

int a=9,b=5,c;

c=a|b;

printf("a=%d\nb=%d\nc=%d\n",a,b,c);

}

3 按位异或运算 按位异或运算符“^”是双目运算符。其功能是参与运算的两数各对应的二进位相异或，当两对应的二进位相异时，结果为1。参与运算数仍以补码出现，例如9^5可写成算式如下： 00001001^00000101 00001100 (十进制为12)

main(){

int a=9;

a=a^15;

printf("a=%d\n",a);

}

4 求反运算 求反运算符～为单目运算符，具有右结合性。 其功能是对参与运算的数的各二进位按位求反。例如～9的运算为： ~(0000000000001001)结果为：1111111111110110

5 左移运算 左移运算符“<<”是双目运算符。其功能把“<< ”左边的运算数的各二进位全部左移若干位，由“<<”右边的数指定移动的位数，

高位丢弃，低位补0。例如： a<<4 指把a的各二进位向左移动4位。如a=00000011(十进制3)，左移4位后为00110000(十进制48)。6 右移运算 右移运算符“>>”是双目运算符。其功能是把“>> ”左边的运算数的各二进位全部右移若干位，“>>”右边的数指定移动的位数。

例如：设 a=15，a>>2 表示把000001111右移为00000011(十进制3)。 应该说明的是，对于有符号数，在右移时，符号位将随同移动。当为正数时， 最高位补0，而为负数时，符号位为1，最高位是补0或是补1 取决于编译系统的规定。Turbo C和很多系统规定为补1。

main(){

unsigned a,b;

printf("input a number: ");

scanf("%d",&a);

b=a>>5;

b=b&15;

printf("a=%d\tb=%d\n",a,b);

}

请再看一例!

main(){

char a='a',b='b';

int p,c,d;

p=a;

p=(p<<8)|b;

d=p&0xff;

c=(p&0xff00)>>8;

printf("a=%d\nb=%d\nc=%d\nd=%d\n",a,b,c,d);

}

位域

有些信息在存储时，并不需要占用一个完整的字节， 而只需占几个或一个二进制位。例如在存放一个开关量时，只有0和1 两种状态， 用一位二进位即可。为了节省存储空间，并使处理简便，C语言又提供了一种数据结构，称为“位域”或“位段”。所谓“位域”是把一个字节中的二进位划分为几个不同的区域， 并说明每个区域的位数。每个域有一个域名，允许在程序中按域名进行 *** 作。 这样就可以把几个不同的对象用一个字节的二进制位域来表示。一、位域的定义和位域变量的说明位域定义与结构定义相仿，其形式为：

struct 位域结构名

{ 位域列表 };

其中位域列表的形式为： 类型说明符 位域名：位域长度

例如：

struct bs

{

int a:8;

int b:2;

int c:6;

};

位域变量的说明与结构变量说明的方式相同。 可采用先定义后说明，同时定义说明或者直接说明这三种方式。例如：

struct bs

{

int a:8;

int b:2;

int c:6;

}data;

说明data为bs变量，共占两个字节。其中位域a占8位，位域b占2位，位域c占6位。对于位域的定义尚有以下几点说明：

1 一个位域必须存储在同一个字节中，不能跨两个字节。如一个字节所剩空间不够存放另一位域时，应从下一单元起存放该位域。也可以有意使某位域从下一单元开始。例如：

struct bs

{

unsigned a:4

unsigned :0 /空域/

unsigned b:4 /从下一单元开始存放/

unsigned c:4

}

在这个位域定义中，a占第一字节的4位，后4位填0表示不使用，b从第二字节开始，占用4位，c占用4位。

2 由于位域不允许跨两个字节，因此位域的长度不能大于一个字节的长度，也就是说不能超过8位二进位。

3 位域可以无位域名，这时它只用来作填充或调整位置。无名的位域是不能使用的。例如：

struct k

{

int a:1

int :2 /该2位不能使用/

int b:3

int c:2

};

从以上分析可以看出，位域在本质上就是一种结构类型， 不过其成员是按二进位分配的。

二、位域的使用位域的使用和结构成员的使用相同，其一般形式为： 位域变量名·位域名 位域允许用各种格式输出。

main(){

struct bs

{

unsigned a:1;

unsigned b:3;

unsigned c:4;

} bit,pbit;

bita=1;

bitb=7;

bitc=15;

printf("%d,%d,%d\n",bita,bitb,bitc);

pbit=&bit;

pbit->a=0;

pbit->b&=3;

pbit->c|=1;

printf("%d,%d,%d\n",pbit->a,pbit->b,pbit->c);

}

上例程序中定义了位域结构bs，三个位域为a,b,c。说明了bs类型的变量bit和指向bs类型的指针变量pbit。这表示位域也是可以使用指针的。

程序的9、10、11三行分别给三个位域赋值。( 应注意赋值不能超过该位域的允许范围)程序第12行以整型量格式输出三个域的内容。第13行把位域变量bit的地址送给指针变量pbit。第14行用指针方式给位域a重新赋值，赋为0。第15行使用了复合的位运算符"&="， 该行相当于： pbit->b=pbit->b&3位域b中原有值为7，与3作按位与运算的结果为3(111&011=011,十进制值为3)。同样，程序第16行中使用了复合位运算"|="， 相当于： pbit->c=pbit->c|1其结果为15。程序第17行用指针方式输出了这三个域的值。

类型定义符typedef

C语言不仅提供了丰富的数据类型，而且还允许由用户自己定义类型说明符，也就是说允许由用户为数据类型取“别名”。 类型定义符typedef即可用来完成此功能。例如，有整型量a,b,其说明如下： int aa,b; 其中int是整型变量的类型说明符。int的完整写法为integer，

为了增加程序的可读性，可把整型说明符用typedef定义为： typedef int INTEGER 这以后就可用INTEGER来代替int作整型变量的类型说明了。 例如： INTEGER a,b;它等效于： int a,b; 用typedef定义数组、指针、结构等类型将带来很大的方便，不仅使程序书写简单而且使意义更为明确，因而增强了可读性。例如：

typedef char NAME[20]; 表示NAME是字符数组类型，数组长度为20。

然后可用NAME 说明变量，如： NAME a1,a2,s1,s2;完全等效于： char a1[20],a2[20],s1[20],s2[20]

又如：

typedef struct stu{ char name[20];

int age;

char sex;

} STU;

定义STU表示stu的结构类型，然后可用STU来说明结构变量： STU body1,body2;

typedef定义的一般形式为： typedef 原类型名 新类型名 其中原类型名中含有定义部分，新类型名一般用大写表示， 以

便于区别。在有时也可用宏定义来代替typedef的功能，但是宏定义是由预处理完成的，而typedef则是在编译时完成的，后者更为灵活方便。

我们学校的数学建模上机课也有Mathlab程序，看看下面有没有你要找的。

一 基本运算

1 求

输入(12+2(7-4))/3^2执行

2 输入x = (52+13-08)10^2/25执行

再输入y= 2x+1执行

3 执行clear命令。观察结果

4计算圆面积Area = ,半径r = 2，则可键入

r=2;area=pir^2; area

问：语句末尾加分号与不加分号有何区别？请试验之

5常用函数

名称 含义 名称 含义

sin 正弦 exp E为底的指数

cos 余弦 log 自然对数

tan 正切 log10 10为底的对数

cot 余切 log2 2为底的对数

asin 反正弦 abs 绝对值

acos 反余弦

例：1）执行y = sin(10)exp(-034^2)

2） 想计算 的值

输入y1=2sin(03pi)/(1+sqrt(5))执行之

若又想计算 ，可以简便地用 *** 作：先按á键则会出现上面输入过的指令 y1=2sin(03pi)/(1+sqrt(5)) ；然后移动光标，把y1改成y2；把 sin 改成 cos 便可。即得

y2=2cos(03pi)/(1+sqrt(5))然后执行之。

系统默认4位有效数字，若想提高精度则可如下：

digits(10);sym(y2,'d') 执行就可精确到小数点后10位，还可将10改为其它数字试验

二 矩阵运算

1要得到矩阵 ，

可输入A = [1,2,3; 4,5,6; 7,8,9] 执行，观察结果

还可分行输入

A=[1,2,3

4,5,6

7,8,9]

效果相同

2 注意 %号后的语句为注释，练习时不必输入

>>a=[1,4,6,8,10] %一维矩阵

>>a(3) % a的第三个元素

ans =

6

»x =[1 2 3 4 5 6 7 8

4 5 6 7 8 9 10 11]; %二维2x8 矩阵

执行后双击左边Workspace里的x，观察之

» x(3) % x的第三个元素

ans =

2

» x([1 2 5]) % x的第一、二、五个元素

ans =

1 4 3

如需要还可定义b=x([1 2 5])执行后结果为

b =

1 4 3

>> x(2,3) % x的第二行第三列的元素

ans =

6

x(1:5) % x的第前五个元素

ans =

1 4 2 5 3

» x(10:end) % x的第十个元素后的元素

ans =

8 6 9 7 10 8 11

执行后双击左边Workspace里的x，观察是哪十个元素

» x(10:-1:2) % x的第十个元素和第二个元素的倒排

ans =

8 5 7 4 6 3 5 2 4

» x(find(x>5)) % x中大于5的元素

ans =

6 7 8 6 9 7 10 8 11

» x(4)=100 %给x的第四个元素重新给值

x =

1 2 3 4 5 6 7 8

4 100 6 7 8 9 10 11

» x(3)=[] % 删除第三个元素（不是二维数组）

x =

Columns 1 through 12

1 4 100 3 6 4 7 5 8 6 9 7

Columns 13 through 15

10 8 11

» x(16)=1 % 加入第十六个元素

x =

Columns 1 through 12

1 4 100 3 6 4 7 5 8 6 9 7

Columns 13 through 16

10 8 11 1

3 如不需要以前的变量时，为不干扰以后计算，可执行clear清除以前的变量

当元素很多的时候，则须采用以下的方式：

» x=(1:2:121); % 以起始值为1,增量值为2,终止值为121的矩阵

» x=linspace(0,1,100); % 利用linspace，生成以0为起始值,1为终止值,元素数目为100的矩阵

»a=[] %空矩阵

a =

[]

» zeros(2,2) %全为0的矩阵

ans =

0 0

0 0

» ones(3,3) %全为1的矩阵

ans =

1 1 1

1 1 1

1 1 1

» rand(2,4); %随机矩阵

4另外一种定义矩阵的方式

»a=1:7; b=11:2:23;

»c=[b a]; %利用上面建立的阵列 a 及阵列 b ，组成新阵列c

»d=[b ; a]; %利用a及b，组成新矩阵d

执行后双击左边Workspace里的c与d，比较之

再如 已知y=[-1,6,15,7,31,2,4,5];

x=y(3:5) %x为y的第三到第五个元素组成的新向量

或 x=[y(5),y(3),y(7)] %x为y的第五、第三、第七个元素组成的新向量

或这样更简单 x=y([5,3,7])

5 输入矩阵x=[4,8,12,10,23;6,3,15,13,19;9,1,2,18,14;11,7,5,21,17]

依次输入下列命令并执行，观察结果，各命令分别有什么作用？

max(x)

min(x) （问：如何得到整个矩阵的最小值与最大值？）

[m,n]=size(x)

L=length(x)

y=x’

a=x( :,2)

b=x( :,2)’

c=x(3, :)

d=x(1 :3,3 :5)

y(2,3)=y(2,3)/2

y(2, :)=y(2, :)/2

y( :,4)=y( :,4)+y( :,2)

6 点运算 执行下列命令，指出点运算的作用

x=1 :8 （或对另外的向量或矩阵来作）

y=2^x

z=x/y

w=x^2

u=sin(x)

常用命令

min 最小值 max 最大值

mean 平均值 std 标准差

sort 排序 diff 相邻元素的差

length 个数 sum 总和

dot 内积 cross 外积

三 画图

二维图形

命 令 含 义 plot绘图函数的叁数

plot 建立向量或矩阵各队队向量的图形 字元 颜色 字元 图线型态

loglog x、y轴都取对数标度建立图形 y ** 点

semilogx x轴用于对数标度，y轴线性标度绘制图形 k 黑色 o 圆

semilogy y轴用于对数标度，x轴线性标度绘制图形 w 白色 x x

title 给图形加标题 b 蓝色 + +

xlabel 给x轴加标记 g 绿色

ylabel 给y轴加标记 r 红色 - 实线

text 在图形指定的位置上加文本字符串 c 亮青色 : 点线

gtext 在鼠标的位置上加文本字符串 m 锰紫色 - 点虚线

grid 打开网格线 -- 虚线

hold on 命令用于在已画好的图形上添加新的图形

1 x=0:0001:10; % 0到10的1000个点(每隔0001画一个点)的x座标

y=sin(x); % 对应的y座标

plot(x,y); % 绘图

注：matlab画图实际上就是描点连线，因此如果点取得不密，画出来就成了折线图，请试验之

2 Y=sin(10x);

plot(x,y,'r:',x,Y,'b') % 同时画两个函数

3 若要改变颜色，在座标对后面加上相关字串即可：

x=0:001:10;

plot(x,sin(x),'r')

4 若要同时改变颜色及图线型态（Line style），也是在坐标对后面加上相关字串即可：

plot(x,sin(x),'r')

5 用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围

axis([0,6,-15,1])

6 MATLAB也可对图形加上各种注解与处理：（见上表）

xlabel('x轴'); % x轴注解

ylabel('y轴'); % y轴注解

title('余弦函数'); % 图形标题

legend('y = cos(x)'); % 图形注解

gtext('y = cos(x)'); % 图形注解 ,用鼠标定位注解位置

grid on; % 显示格线

7画椭圆

a = [0:pi/50:2pi]'; %角度

X = cos(a)3; %参数方程

Y = sin(a)2;

plot(X,Y);

xlabel('x'), ylabel('y');

title('椭圆')

8 绘制函数 在0 ≤ x ≤ 1时的曲线。

x=0:01:1

y=xexp(-x) %为什么用点运算？若不用会怎样

plot(x,y),xlabel('x'),ylabel('y'),title('y=xexp(-x)')

9 画出衰减振荡曲线 与它的包络线 及 。t 的取值范围是[0, 4π] 。

t=0:pi/50:4pi;

y0=exp(-t/3);

y=exp(-t/3)sin(3t);

plot(t,y,'-r',t,y0,':b',t,-y0,':b') % -r表示红色实线，:b表示蓝色点线，看上表

grid

10 在同一个画面上建立几个坐标系, 用subplot(m,n,p)命令；把一个画面分成m×n个图形区域, p代表当前的区域号，在每个区域中分别画一个图,如

x=linspace(0,2pi,30); y=sin(x); z=cos(x);

u=2sin(x)cos(x); v=sin(x)/cos(x);

subplot(2,2,1),plot(x,y),axis([0 2pi -1 1]),title('sin(x)')

subplot(2,2,2),plot(x,z),axis([0 2pi -1 1]),title('cos(x)')

subplot(2,2,3),plot(x,u),axis([0 2pi -1 1]),title('2sin(x)cos(x)')

subplot(2,2,4),plot(x,v),axis([0 2pi -20 20]),title('sin(x)/cos(x)')

三维图形

11三维螺旋线：

t=0:pi/50:10pi;

plot3(sin(t),cos(t),t) %参数方程

grid %添加网格

12 t=linspace(0,20pi, 501);

plot3(tsin(t), tcos(t), t); %注意点乘

也可以同时画出两条曲线，格式与二维情况类似，兹不举例。

13用mesh命令画曲面

画出由函数 形成的立体网状图:

a=linspace(-2, 2, 25); % 在x轴上从（-2，2）取25点

b=linspace(-2, 2, 25); % 在y轴上取25点

[x,y]=meshgrid(a, b); % x和y都是21x21的矩阵

z=xexp(-x^2-y^2); % 计算函数值，z也是21x21的矩阵

mesh(x, y, z); % 画出立体网状图

14 surf和mesh的用法类似：

a=linspace(-2, 2, 25); % 在x轴上取25点

b=linspace(-2, 2, 25); % 在y轴上取25点

[x,y]=meshgrid(a, b); % x和y都是21x21的矩阵

z=xexp(-x^2-y^2); % 计算函数值，z也是21x21的矩阵

surf(x, y, z); % 画出立体曲面图

四 程序设计

1 M-文件： 上面所做的运算都是在命令窗口中输入一条或两三条命令，然后执行，再输入，再执行，以这样交谈式的方式进行。如果为了解决某一问题需要很多命令，这样做就很不方便了。这时我们把解决某一问题的所有命令集中放在一个文档里，命名、保存。然后只要在命令窗口中输入文档名，执行即可。

例:（1）编写文档:点击MATLAB指令窗口上面最左端的图标 ，即新建文件，就可打开MATLAB文件编辑器。用户即可在空白窗口中编写程序。例如输入下面的程序：

x=linspace(0,2pi,20);

y=sin(x);

plot(x,y,'r+')

title('2D plot')

(2）点击文件编辑器上面工具条中的保存 ，命名（例如将上面的程序命名为picture），然后保存。像这样在MATLAB文件编辑器中编写的文件叫M-文件（M-file）。

(3)运行：i）在命令窗口中输入文件名(如上面的picture),然后执行。

ii）或直接在文件编辑器上面的工具条中找到debug（即调试），点击，再找到run（即运行），再点击即可。

同学们可以把前面画图的一些问题放在文件编辑器里再做一下。

2 自己编写函数：我们经常用到的像sin、cos、exp这样的一些函数都是MATLAB软件自身所带的函数，因此直接应用即可，但有时我们为了解决一些问题需要自己编写函数。自己编写函数有两个基本要求i）必须在MATLAB文件编辑器中编写。ii）函数名和文件名必须相同。 例: 编写函数 , 计算f(1)f(2)+f2(3)

（1）打开MATLAB文件编辑器，输入

function Y= fun1(x) % 表示Y是x的函数，x是自变量，fun1是函数名

Y=(x^3 - 2x^2 + x - 63)/(x^2 + 005x - 314);

然后保存。

注：在自己编写的函数前都要写上function，表示这是自己定义的函数。fun1表示函数名，那么最后文件名也应命名为fun1。

（2）这样在命令窗口中就可以像应用sin、cos那样来使用函数fun1，如：在命令窗口中输入 >> fun1(1)fun1(2)+fun1(3)fun1(3) 结果为：

ans =

-126023

3 for循环语句（这里的for语句与C语言中的for语句不同，要更简单一些）

例：一个简单的for循环示例。

for i=1:10; % i依次取1,2,…10,

x(i)=2i; % 对每个i值，重复执行该指令

end; % 表示循环结束，每一个for要对应一个end

x % 要求显示运行后数组x的值。

输入后观察结果，体会for语句的作用。

注：在MATLAB里（在C语言中也一样）， 的作用表示把等号右边的值送给左边的变量，这和数学中相等的意思不同。下面的例子中都要这样理解，否则就不能明白程序的含义。

4 while循环语句

例: Fibonacci 数列：1，1，2，3，5，8，… 即： ，（ 1，2，3…）现要求该数列中第一个大于10000 的元素。

a(1)=1;a(2)=1;i=2;

while a(i)<=10000

a(i+1)=a(i-1)+a(i);

i=i+1;

end;

i,a(i),

5（1）if-end语句，例：

cost=10;number=12;

if number>8

sums=number095cost;

end,

sums

（2）if-else-end语句，例：

cost=10;number=5; % 改变number的初值，看结果有何不同

if number>8

sums=number095cost;

else sums=number05cost;

end,

sums

6 例：用for 循环语句来寻找Fibonacc 数列中第一个大于10000 的元素。

n=100;a=ones(1,n); % a是一个一行，n列的所有元素为1的矩阵

for i=3:n

a(i)=a(i-1)+a(i-2);

if a(i)>=10000

a(i),

break; % 表示跳出循环

end;

end, i

7 练习：课本264页，参考例4右边的流程图114，编程序求解例4，自己设置误差，并与书上的结果比较。

五 拟合与插值

曲 线 拟 合 和 插 值 函 数

polyfit(x, y, n) 对描述n阶多项式y=f(x)的数据进行最小二乘曲线拟合

interp1(x, y, xo) 1维线性插值

interp1(x, y, xo, ' spline ') 1维3次样条插值

interp1(x, y, xo, ' cubic ') 1维3次插值

interp2(x, y, Z, xi, yi) 2维线性插值

interp2(x, y, Z, xi, yi, ' cubic') 2维3次插值

1 插值

看课本266页§112第一段，了解什么是插值。

例：考虑下列问题，12小时内，一小时测量一次室外温度。数据如下：

时间：1，2，3， 4， 5， 6， 7， 8， 9，10，11，12

温度：5，8，9，15，25，29，31，30，22，25，27，24

现在根据以上数据估计32，47等时刻的温度

hours=1:12;

temps=[5 8 9 15 25 29 31 30 22 25 27 24];

t=interp1(hours, temps, [32，47]) % 一阶线性插值，如果只估计一个点的值，则无须加方括号

改为t=interp1(hours, temps, [32，47], 'spline'） 则为三次样条插值

如果输入如下程序，则画出插值曲线

hours=1:12;

temps=[5 8 9 15 25 29 31 30 22 25 27 24];

h=1:01:12;

t=interp1(hours, temps, h) ; % h后加上'spline'则为三次样条插值

plot(hours, temps, ' + ' , h, t)

用一阶线性插值和三次样条插值做课本268页例2，与书上之结果比较，然后挑课后题做一两道。

2 拟合

看课本270页§113，曲线拟合，比较拟合与插值有什么区别。

例：两组数据如下：

x=[0 1 2 3 4 5 6 7 8 9 1];

y=[-447 1978 328 616 708 734 766 956 948 930 112];

n=8;

p1=polyfit(x,y,n); % n表示用n阶多项式拟合，n=1为线性拟合，即通常所说最小二乘法

poly2sym(p1) % 前面的拟合命令只给出多项式的系数，用此命令则将结果转化为真正的多项式。或用 vpa(poly2sym(p1),10) 即取数值形式,取10位有效数字

x1=0:01:1; % 由此以后三句是画出拟合曲线的图像

y1=polyval(p1,x1); %此句是在x1这些点处求出多项式的值，送给y1

plot(x,y,'o',x1,y1)

改变n的数字，即用不同的多项式拟合，看看哪个结果好。

当n=10时，数据点之间出现大的波动。当企图进行高阶曲线拟合时，这种波动现象经常发生，并不利于我们认识两组数据之间的规律，因此并不是阶数越高越好，实际问题当中，适当选一个即可。

用上面的指令做课本271页例1及例2，将结果与书上之结果比较。

从我的程序里截出来的，有问题再追问吧

switch operator

case 1

for j = 1:num%矩阵 r 的列数

temp(:,j) = min(r(:,j),weight);

end

gar = max(temp);

case 2

for j = 1:num

temp(:,j) =r(:,j)wei;

end

gar = max(temp)

case 3

for j = 1:num

temp(:,j) = min(r(:,j),wei);

end

gar = sum(temp)

case 4

for j = 1:num

temp(:,j) = r(:,j)wei;

end

gar = sum(temp)

end

以前写的程序

#include<stdioh>

#include<stringh>

#define MAX 200

typedef char datetype;

typedef struct Stack

{

datetype elem[MAX];

int top;

}S;//定义存放运算符的栈

typedef struct stack

{

int num[MAX];

int top;

}ss;//定义存放数字的栈

int degree(char member)

{

if(member==''||member=='/')

return 1;

if(member=='+'||member=='-')

return 0;

return -1;

}//返回运算符各自的优先级

int number(char ch)

{

if(ch>='0'&&ch<='9')

return 1;

return 0;

}//判断是否是数字

int prior(char member1,char member2)

{

return degree(member1)-degree(member2);

}//返回两个运算符那个优先

void Init(S s)

{

s->top=0;

}//初始化栈，即对top赋0。

int Full(S s)

{

if(s->top==MAX)

return 1;

return 0;

}//判断是否栈满

int Empty(S s)

{

if(s->top==0)

return 1;

return 0;

}//判断是否空栈

/入栈

1判断是否栈满

2未满则将元素入栈，top++

/

int Push(S s,datetype x)

{

if(Full(s))

return 0;

s->elem[s->top++]=x;

return 1;

}

int push(ss s,int x)

{

if(s->top==MAX)

return 0;

s->num[s->top++]=x;

return 1;

}

/出栈

1判断是否空栈

2非空则top--。

/

int Pop(S s)

{

if(Empty(s))

return 0;

s->top--;

return 1;

}

int pop(ss s)

{

if(s->top==0)

return 0;

s->top--;

return 1;

}

/获取栈顶元素

1判断是否空栈

2非空则返回栈顶元素。

/

datetype GetTop(S s)

{

if(Empty(s))

{

return 0;

}

return s->elem[s->top-1];

}

int gettop(ss s)

{

if(s->top==0)

return 0;

return s->num[s->top-1];

}

/

1判断得到的运算符是什么类型

2根据运算符来计算并返回结果。

/

int Count(int num1,int num2,char ch)

{

if(ch=='')

return num1num2;

if(ch=='/')

return num2/num1;

if(ch=='+')

return num1+num2;

if(ch=='-')

return num2-num1;

return 0;

}

/

1将所有'['和']'转化为'('和')'。(因为其效果一样)

2遇到(则入栈，遇到运算符若栈为空或栈顶元素为(则入栈，栈顶元素为其他运算符则按3的方法。

3两个运算符优先级相同或后者小于前者，则d出前者，并d出num最后两个元素进行运算得到结果入栈，否则入栈

4遇到)则不停d出运算符直至遇到一个(,同时运算得到结果入栈

/

int translate(char str[])

{

int ans=0,len,n,m;

int i;

S sign;

ss num;

numtop=0;

Init(&sign);

len=strlen(str);

printf("转化为后缀为 :");

for(i=0;i<len;i++)

{

if(str[i]=='[') str[i]='(';

if(str[i]==']') str[i]=')';

}

for(i=0;i<len;i++)

{

if(number(str[i]))

{

printf("%c",str[i]);

m=str[i]-'0';

if(i>0&&number(str[i-1]))

{

n=gettop(&num);

pop(&num);

push(&num,n10+m);

}

else

{

push(&num,m);

}

if(i<len-1&&!number(str[i+1]))

printf(" ");

}

else

{

int n1,n2,nn;

char ch;

if(str[i]=='(')

{

Push(&sign,str[i]);

}

else if(str[i]==')')

{

while(GetTop(&sign)!='(')

{

ch=GetTop(&sign);

n1=gettop(&num);

pop(&num);

n2=gettop(&num);

pop(&num);

nn=Count(n1,n2,ch);

push(&num,nn);

printf("%c ",ch);

Pop(&sign);

}

Pop(&sign);

}

else

{

if(!Empty(&sign)||GetTop(&sign)=='(')

{

char ch1=str[i];

int h;

while(!Empty(&sign)||GetTop(&sign)=='(')

{

char ch=GetTop(&sign);

h=prior(ch1,ch);

if(h==1)

{

Push(&sign,ch1);

break;

}

if(h==0||h==-1)

{

int n1,n2;

n1=gettop(&num);

pop(&num);

n2=gettop(&num);

pop(&num);

nn=Count(n1,n2,ch);

push(&num,nn);

printf("%c ",ch);

Pop(&sign);

}

}

if(h!=1) Push(&sign,ch1);

}

else

{

Push(&sign,str[i]);

}

}

}

}

while(!Empty(&sign))

{

char ch=GetTop(&sign);

int n1,n2,nn;

n1=gettop(&num);

pop(&num);

n2=gettop(&num);

pop(&num);

nn=Count(n1,n2,ch);

push(&num,nn);

printf("%c ",ch);

Pop(&sign);

}

printf("\n");

ans=gettop(&num);

return ans;

}

/匹配检查

1遇到(,[入栈。

2遇到]，检查栈顶是否为[,遇到),检查栈顶是否为(

3结束后看是否栈为空或(与)、[与]个数是否相等。

/

int right(char str[])

{

int i,len,n1=0,n2=0,m1=0,m2=0;

S s;

Init(&s);

len=strlen(str);

for(i=0;i<len;i++)

{

if(str[i]=='(')

{

Push(&s,str[i]);

n1++;

}

if(str[i]=='[')

{

Push(&s,str[i]);

m1++;

}

if(str[i]==')')

{

if(GetTop(&s)!='(')

return 0;

Pop(&s);

n2++;

}

if(str[i]==']')

{

if(GetTop(&s)!='[')

return 0;

Pop(&s);

m2++;

}

}

if(n1!=n2||m1!=m2) return 0;

return 1;

}

int main()

{

char str[MAX]={"3[2+2(1+5)]-210+4"};

printf("举例如下：\n计算3[2+2(1+5)]-210+4\n");

printf("计算结果为%d\n",translate(str));

while(gets(str))

{

int ans;

if(!right(str))

{

printf("输入有误!\n");

continue;

}

ans=translate(str);

printf("计算结果为%d\n",ans);

}

return 0;

}

以上就是关于下面是关于C++ Sobel算子的程序代码 相关知识了解 语句看不懂 希望大神能帮我注释下 方便我理解 谢谢全部的内容，包括:下面是关于C++ Sobel算子的程序代码 相关知识了解 语句看不懂 希望大神能帮我注释下 方便我理解 谢谢、C语言编写程序，进行位运算。、帮忙写个Matlab小程序等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！