矩阵计算

最新推荐文章于 2021-08-19 16:33:23 发布
最新推荐文章于 2021-08-19 16:33:23 发布
阅读量99 收藏
点赞数
文章标签： python
原文链接：https://my.oschina.net/yanjianhai/blog/93315
版权

2019独角兽企业重金招聘Python工程师标准>>> hot3.png

#include "stdafx.h"
#include <time.h>

#include <windows.h>
#include<iostream>
using namespace std;

void Familiar_2(int m,int ROW)
{
		int chardata[26][26]={
					{4,-2,0,-2,-1,-2,0,-2,-1,0,-1,-1,-1,-2,0,-1,-1,-1,1,0,0,0,-3,0,-2,-1},
					{-2,4,-3,4,1,-3,-1,0,-3,0,0,-4,-3,3,0,-2,0,-1,0,-1,0,-3,-4,-1,-3,1},
					{0,-3,9,-3,-4,-2,-3,-3,-1,0,-3,-1,-1,-3,0,-3,-3,-3,-1,-1,0,-1,-2,-2,-2,-3},
					{-2,4,-3,6,2,-3,-1,-1,-3,0,-1,-4,-3,1,0,-1,0,-2,0,-1,0,-3,-4,-1,-3,1},
					{-1,1,-4, 2,5,-3,-2,0,-3,0,1,-3,-2,0,0,-1,2,0,0,-1,0,-2,-3,-1,-2, 4},
					{-2,-3,-2,-3,-3,6,-3,-1,0,0,-3,0,0,-3,0,-4,-3,-3,-2,-2,0,-1,1,-1,3,-3},
					{0,-1,-3,-1,-2,-3,6,-2,-4,0,-2,-4,-3,0,0,-2,-2,-2,0,-2,0,-3,-2,-1,-3,-2},
					{-2,0,-3,-1,0,-1,-2, 8,-3,0,-1,-3,-2,1,0,-2,0,0,-1,-2,0,-3,-2,-1,2,0},
					{-1,-3,-1,-3,-3,0,-4,-3,4,0,-3,2,1,-3,0,-3,-3,-3,-2,-1,0,3,-3,-1,-1,-3},
					{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
					{-1,0,-3,-1,1,-3,-2,-1,-3,0,5,-2,-1,0,0,-1,1,2,0,-1,0,-2,-3,-1,-2,1},
					{-1,-4,-1,-4,-3,0,-4,-3,2,0,-2,4,2,-3,0,-3,-2,-2,-2,-1,0,1,-2,-1,-1,-3},
					{-1,-3,-1,-3,-2,0,-3,-2,1,0,-1,2,5,-2,0,-2,0,-1,-1,-1,0,1,-1,-1,-1,-1},
					{-2,3,-3,1,0,-3,0,1,-3,0,0,-3,-2,6,0,-2,0,0,1,0,0,-3,-4,-1,-2,0},
					{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
					{-1,-2,-3,-1,-1,-4,-2,-2,-3, 0,-1,-3,-2,-2, 0, 7,-1,-2,-1,-1, 0,-2,-4,-2,-3,-1},
					{ -1, 0,-3, 0, 2,-3,-2, 0,-3, 0, 1,-2, 0, 0, 0,-1, 5, 1, 0,-1, 0,-2,-2,-1,-1, 3},
					{ -1,-1,-3,-2, 0,-3,-2, 0,-3, 0, 2,-2,-1, 0, 0,-2, 1, 5,-1,-1, 0,-3,-3,-1,-2, 0},
					{ 1, 0,-1, 0, 0,-2, 0,-1,-2, 0, 0,-2,-1, 1, 0,-1, 0,-1, 4, 1, 0,-2,-3, 0,-2, 0},
					{ 0,-1,-1,-1,-1,-2,-2,-2,-1, 0,-1,-1,-1, 0, 0,-1,-1,-1, 1, 5, 0, 0,-2, 0,-2,-1},
					{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
					{ 0,-3,-1,-3,-2,-1,-3,-3, 3, 0,-2, 1, 1,-3, 0,-2,-2,-3,-2, 0, 0, 4,-3,-1,-1,-2},
					{ -3,-4,-2,-4,-3, 1,-2,-2,-3, 0,-3,-2,-1,-4, 0,-4,-2,-3,-3,-2, 0,-3,11,-2, 2,-3},
					{ 0,-1,-2,-1,-1,-1,-1,-1,-1, 0,-1,-1,-1,-1, 0,-2,-1,-1, 0, 0, 0,-1,-2,-1,-1,-1},
					{ -2,-3,-2,-3,-2, 3,-3, 2,-1, 0,-2,-1,-1,-2, 0,-3,-1,-2,-2,-2, 0,-1, 2,-1, 7,-2},
					{ -1, 1,-3, 1, 4,-3,-2, 0,-3, 0, 1,-3,-1, 0, 0,-1, 3, 0, 0,-1, 0,-2,-3,-1,-2, 4}};

	int d=-8;
	int i,j,k,row,col,temp_1,temp_2;
	int n=m+1;
	int **M=new int *[n];
	for(i=0;i<n;i++)
	{
		M[i]=new int[n];
	}
	for(i=0;i<n;i++)
	{
		M[0][i]=i*d;
		M[i][0]=i*d;
	}

clock_t t1 = clock();

for(k=2;k<=2*m;k++)
{
		
	if(k<m+1)
	{	
		temp_1=1;
		temp_2=k;
	}
	
	else
	{
		temp_1=k-m;
		temp_2=m+1;
	}
		//加入并行
#pragma omp parallel for schedule(static,16)
		for(i=temp_1;i<temp_2;i++)
		{
			j=k-i;
			M[i][j]=M[i-1][j-1]+chardata[(i-1)%26][(j-1)%26];
			
			for(col=0;col<j;col++)//计算行最大
				if(M[i][col]+(j-col)*(j-col)>M[i][j])
					M[i][j]=M[i][col]+(j-col)*(j-col);

			
			
			for(row=0;row<i;row++)//计算列最大
				if(M[row][j]+(i-row)*(i-row)>M[i][j])
					M[i][j]=M[row][j]+(i-row)*(i-row);

			
			
			
		}


}
clock_t t2 = clock();
int t3=t2-t1;
cout<<t3<<endl;

	for(i=0;i<n;i++)
			cout<<M[ROW][i]<<" ";
	cout<<endl;



	
		for(i=0;i<n;i++)
		{
			delete [] M[i];
		}
		delete [] M;

}




void Familiar_1(int m,int ROW)
{
	
	
	
	int chardata[26][26]={
					{4,-2,0,-2,-1,-2,0,-2,-1,0,-1,-1,-1,-2,0,-1,-1,-1,1,0,0,0,-3,0,-2,-1},
					{-2,4,-3,4,1,-3,-1,0,-3,0,0,-4,-3,3,0,-2,0,-1,0,-1,0,-3,-4,-1,-3,1},
					{0,-3,9,-3,-4,-2,-3,-3,-1,0,-3,-1,-1,-3,0,-3,-3,-3,-1,-1,0,-1,-2,-2,-2,-3},
					{-2,4,-3,6,2,-3,-1,-1,-3,0,-1,-4,-3,1,0,-1,0,-2,0,-1,0,-3,-4,-1,-3,1},
					{-1,1,-4, 2,5,-3,-2,0,-3,0,1,-3,-2,0,0,-1,2,0,0,-1,0,-2,-3,-1,-2, 4},
					{-2,-3,-2,-3,-3,6,-3,-1,0,0,-3,0,0,-3,0,-4,-3,-3,-2,-2,0,-1,1,-1,3,-3},
					{0,-1,-3,-1,-2,-3,6,-2,-4,0,-2,-4,-3,0,0,-2,-2,-2,0,-2,0,-3,-2,-1,-3,-2},
					{-2,0,-3,-1,0,-1,-2, 8,-3,0,-1,-3,-2,1,0,-2,0,0,-1,-2,0,-3,-2,-1,2,0},
					{-1,-3,-1,-3,-3,0,-4,-3,4,0,-3,2,1,-3,0,-3,-3,-3,-2,-1,0,3,-3,-1,-1,-3},
					{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
					{-1,0,-3,-1,1,-3,-2,-1,-3,0,5,-2,-1,0,0,-1,1,2,0,-1,0,-2,-3,-1,-2,1},
					{-1,-4,-1,-4,-3,0,-4,-3,2,0,-2,4,2,-3,0,-3,-2,-2,-2,-1,0,1,-2,-1,-1,-3},
					{-1,-3,-1,-3,-2,0,-3,-2,1,0,-1,2,5,-2,0,-2,0,-1,-1,-1,0,1,-1,-1,-1,-1},
					{-2,3,-3,1,0,-3,0,1,-3,0,0,-3,-2,6,0,-2,0,0,1,0,0,-3,-4,-1,-2,0},
					{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0},
					{-1,-2,-3,-1,-1,-4,-2,-2,-3, 0,-1,-3,-2,-2, 0, 7,-1,-2,-1,-1, 0,-2,-4,-2,-3,-1},
					{ -1, 0,-3, 0, 2,-3,-2, 0,-3, 0, 1,-2, 0, 0, 0,-1, 5, 1, 0,-1, 0,-2,-2,-1,-1, 3},
					{ -1,-1,-3,-2, 0,-3,-2, 0,-3, 0, 2,-2,-1, 0, 0,-2, 1, 5,-1,-1, 0,-3,-3,-1,-2, 0},
					{ 1, 0,-1, 0, 0,-2, 0,-1,-2, 0, 0,-2,-1, 1, 0,-1, 0,-1, 4, 1, 0,-2,-3, 0,-2, 0},
					{ 0,-1,-1,-1,-1,-2,-2,-2,-1, 0,-1,-1,-1, 0, 0,-1,-1,-1, 1, 5, 0, 0,-2, 0,-2,-1},
					{ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0},
					{ 0,-3,-1,-3,-2,-1,-3,-3, 3, 0,-2, 1, 1,-3, 0,-2,-2,-3,-2, 0, 0, 4,-3,-1,-1,-2},
					{ -3,-4,-2,-4,-3, 1,-2,-2,-3, 0,-3,-2,-1,-4, 0,-4,-2,-3,-3,-2, 0,-3,11,-2, 2,-3},
					{ 0,-1,-2,-1,-1,-1,-1,-1,-1, 0,-1,-1,-1,-1, 0,-2,-1,-1, 0, 0, 0,-1,-2,-1,-1,-1},
					{ -2,-3,-2,-3,-2, 3,-3, 2,-1, 0,-2,-1,-1,-2, 0,-3,-1,-2,-2,-2, 0,-1, 2,-1, 7,-2},
					{ -1, 1,-3, 1, 4,-3,-2, 0,-3, 0, 1,-3,-1, 0, 0,-1, 3, 0, 0,-1, 0,-2,-3,-1,-2, 4}};

	int d=-8;
	int i,j,k,row,col,temp_1,temp_2;
	int n=m+1;
	int **M=new int *[n];
	for(i=0;i<n;i++)
	{
		M[i]=new int[n];
	}
	for(i=0;i<n;i++)
	{
		M[0][i]=i*d;
		M[i][0]=i*d;
		
	}

	



clock_t t1 = clock();

for(k=2;k<=2*m;k++)
{
		
	if(k<m+1)
	{	
		temp_1=1;
		temp_2=k;
	}
	
	else
	{
		temp_1=k-m;
		temp_2=m+1;
	}
		//加入并行
#pragma omp parallel for 
		for(i=temp_1;i<temp_2;i++)
		{
			j=k-i;
			M[i][j]=M[i-1][j-1]+chardata[(i-1)%26][(j-1)%26];
			
			for(col=0;col<j;col++)//计算行最大
				if(M[i][col]+d>M[i][j])
					M[i][j]=M[i][col]+d;

			
			
			for(row=0;row<i;row++)//计算列最大
				if(M[row][j]+d>M[i][j])
				M[i][j]=M[row][j]+d;

		
		}


}
clock_t t2 = clock();
int t3=t2-t1;
cout<<t3<<endl;

			
	for(i=0;i<n;i++)
			cout<<M[ROW][i]<<" ";
	cout<<endl;



		for(i=0;i<n;i++)
		{
			delete [] M[i];
		}
		delete [] M;
	
		
}




int main(){
	int demension,row;
	cout<<"这个程序输出两个序列的相似度，并且算出相似度矩阵，请输入维数："<<endl;
	cin>>demension;
	cout<<"请输入要打印出来的行数："<<endl;
	cin>>row;
	Familiar_1(demension,row);
	//Familiar_2(demension,row);
	return 0;
}

转载于:https://my.oschina.net/yanjianhai/blog/93315

weixin_33697898
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Golub《Matrix Computations》Chapter 5 学习
weixin_43745120的博客
03-04 474
  因为在看SAM和ISAM，其中有些推导不太懂，所以想学一下Golub这本书中的第5章“Orthogonalization and Least Squares”。 5.1 Householder and Givens Matrices   首先回顾正交(orthogonal)的定义,如果QTQ=QQT=In，Q∈Rn∗nQ^TQ=QQ^T=I_n，Q\in R^{n*n}QTQ=QQT=In​，...
矩阵计算第3版 [（美）戈卢布，（美）范洛恩 著] 2011年版.rar
06-01
本书是数值计算领域的名著，系统地介绍了矩阵计算的基本理论和方法。内容包括：矩阵乘法、矩阵分析、线性方程组、正交化和最小二乘法、特征值问题、Lanczos方法、矩阵函数及专题讨论等。书中的许多算法都由现成的软件包来实现，每节后还附有习题，并有注释和大量参考文献
矩阵计算-中文版-Gene H. Golub-袁亚湘-人民邮电出版社.pdf
07-18
本书是数值计算领域的名著，系统介绍了矩阵计算的基本理论和方法。内容包括：矩阵乘法、矩阵分析、线性方程组、正交化和最小二乘法、特征值问题、Lanczos 方法、矩阵函数及专题讨论等。书中的许多算法都有现成的软件包实现，每节后附有习题，并有注释和大量参考文献。新版增加约四分之一内容，反映了近年来矩阵计算领域的飞速发展。
增加动手环节-将Golub的矩阵计算代码封装在matlab 函数文件和CUDA lib-第一章
eloudy的专栏
08-19 214
一本值得百炼的计算数学书籍 突发奇想，这次把将Golub的矩阵计算代码封装在matlab 函数文件中应用。 书中代码是为普通代码场景预备的，并不适合直接封装成函数，但是可以在函数中使用它，虽然会略显不自然。 algorithm1.1.1: function c = mc_innerproduct(x,y) %Matrix Computations 4th edition Chapter1 algorithm 1.1.1 c =0; for i=1:numel(x)
矩阵计算_labview矩阵计算_labview矩阵运算_
10-04
在LabVIEW中，矩阵计算是核心功能之一，因为矩阵和数组是处理复杂数据和执行高级数学运算的基础。这篇描述提到的是一个LabVIEW项目，其中包含多个子VI（虚拟仪器），专门用于矩阵的计算和操作。 1. **矩阵创建与...
矩阵计算matrix实用.pdf
12-05
矩阵计算Matrix实用.pdf 矩阵计算Matrix实用.pdf是Maplesoft公司开发的一款功能强大的矩阵计算工具，提供了多种矩阵操作和计算功能，包括矩阵的基本操作、线性代数计算、矩阵分解、矩阵乘法、矩阵求逆、矩阵行列式...
test_矩阵计算小程序_
10-03
在IT领域，矩阵计算是计算机科学中的一个重要组成部分，特别是在数值分析、线性代数、图形学、机器学习等众多领域有着广泛的应用。本程序“test_矩阵计算小程序_”是用Java语言编写的一个小型应用程序，专门用于进行...
cav_矩阵计算_
09-29
在IT领域，矩阵计算是一种基础且重要的计算方法，特别是在数据科学、机器学习、图像处理以及高性能计算中。本文将深入探讨“cav_矩阵计算_”这一主题，讲解如何将网络上数据网站使用的数据文件格式转换为传统的矩阵...
矩阵计算EJML
08-19
用于Java开发中的矩阵计算
矩阵计算讲义
11-08
矩阵计算是数学、计算机科学和工程领域中的一个重要分支，它涉及到多维度数组的运算，尤其在处理线性方程组、数据分析、图像处理、机器学习等领域有着广泛应用。本讲义将深入探讨矩阵计算的基础概念、理论及其实际...
Desktop_由ABCD矩阵计算RLCG分布参数_
09-30
这里，我们将深入探讨这两个概念及其关系，并了解如何通过ABCD矩阵计算RLCG分布参数。 首先，ABCD矩阵是一种描述线性二端口网络特性的数学工具。这个2x2矩阵包含了网络对电压和电流比例的关系，即输入和输出之间的...
矩阵计算jsp
03-30
jsp 实现的矩阵计算器，转置，逆矩阵，乘法等等。jsp 实现的矩阵计算器，转置，逆矩阵，乘法等等。
vb矩阵计算
12-31
在VB（Visual Basic）编程环境中，矩阵计算是一个重要的数学应用领域，主要涉及到线性代数的概念和算法。本文将深入探讨VB中矩阵计算的相关知识点，包括矩阵的基本操作、行列式计算、逆矩阵求解以及矩阵的转置。此外...
潘建瑜老师的矩阵计算讲义(带有书签)
11-16
《潘建瑜老师的矩阵计算讲义》是一份珍贵的学习资源，专为深入理解矩阵计算而设计。矩阵计算是数学、工程、计算机科学等领域中的基础工具，尤其在数据分析、机器学习和信号处理等现代科技中扮演着至关重要的角色。这...
C++矩阵计算类
01-19
在C++编程中，矩阵计算是数值分析、线性代数和计算机图形学等领域常见的操作。为了方便地处理矩阵，我们可以创建一个专门的C++类来封装这些功能。"C++矩阵计算类"的设计和实现旨在提供高效且灵活的矩阵运算，包括...
MatLab矩阵计算详解
"ChMatLab中的矩阵计算完美版资料.ppt" 本文主要介绍ChMatLab中的矩阵计算基础知识，包括变量及数据操作、矩阵定义、常用数学函数以及数据输出格式。 一、变量及数据操作 在ChMatLab中，变量是用来存储数据的，...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值