matlab2c使用c++实现matlab函数系列教程-sawtooth函数

全栈工程师开发手册 （作者：栾鹏）

matlab2c动态链接库下载
matlab库函数大全
matlab2c基础教程
matlab2c开发全解教程

matlab2c调用方法：

1、下载动态链接库
2、将Matlab2c.dll拷贝到exe同目录下
3、将Matlab2c.h、Matlab2c.lib放到项目头文件目录下
4、在cpp文件中引入下面的代码

#include "Matlab2c.h"
#pragma comment(lib,"Matlab2c.lib")  
using namespace Matlab2c;

---

matlab中sawtooth函数简介

1、sawtooth函数：

产生锯齿波或三角波的函数

2、用法说明

x = sawtooth(t) 第一种调用方式，将产生周期为2π的锯齿波。以0_{2π这个周期内为例，当t=0时，x=-1，当t=2π时，x=1。由此可见，在0}2π这个周期内，x(t)是关于t的以1/π为斜率的线段。
x = sawtooth(t,width) 第二种调用方式中，width是0到1之间的标量。在0到2π×width区间内，x的值从-1线性变化到1；在2π×width~2π区间内，x的值又从1线性变化到-1。sawtooth(t,1)和sawtooth(t)是等价的。

sawtooth的c++源码实现

产生锯齿波或者三角波，幅度为正负1
输入：锯齿波信号时域取值点a，上升段比例width
输出：锯齿波信号在取值点上的值
提示：幅度为正负1，周期为2pi，width为上升段比例,默认为1

Matrix Matlab2c::sawtooth(Matrix& a,double width)
{
	int i,j;
	Matrix p(a.row,a.column);
	int n;
	double data;
	for (i=0;i<a.row;i++)
		for (j=0;j<a.column;j++)
		{
			n=std::floor(a(i,j)/(2*PI));  //向下取整
			data = a(i,j)-n*2*PI;
			if (data<=(2*PI*width) && data>0)
				p(i,j) = 1.0/(PI*width)*data-1;
			else if (data>(2*PI*width) && data<=(2*PI))
				p(i,j) = -1.0/(PI*(1-width))*data+(1+width)/(1-width);
		}
		return p;
}

sawtooth函数的使用测试

#include "Matlab2c.h"
#pragma comment(lib,"Matlab2c.lib")  
using namespace Matlab2c;


int main()
{
	double a[]={1,2,3,4,5,6};
	Matrix aa=Matrix(1,6,a);
	Matrix cc=Matlab2c::sawtooth(aa,0.5);
	cout<<cc.toString()<<endl;

	system("pause");
	return 0;
}