实验源代码 //Warshall.cpp #include<stdio.h> void warshall(int k,int n) { int i , j, t; int temp[20][20]; for(int a=0;a<k;a++) { printf("请输入矩阵第%d 行元素:",a); for(int b=0;b<n;b++) { scanf ("%d",&temp[a][b]); } } for(i=0;i<k;i++){ for( j=0;j<k;j++){ if(temp[ j][i]==1) { for(t=0;t<n;t++) { temp[ j][t]=temp[i][t]||temp[ j][t]; } } } } printf("可传递闭包关系矩阵是:\n"); for(i=0;i<k;i++) { for( j=0;j<n;j++) { printf("%d", temp[i][ j]); } printf("\n"); } } void main() { printf("利用 Warshall 算法求二元关系的可传递闭包\n"); void warshall(int,int); int k , n; printf("请输入矩阵的行数 i: "); scanf("%d",&k); 四川大学实验报告 printf("请输入矩阵的列数 j: "); scanf("%d",&n); warshall(k,n); }
上传时间: 2016-06-27
上传用户:梁雪文以
#include "iostream" using namespace std; class Matrix { private: double** A; //矩阵A double *b; //向量b public: int size; Matrix(int ); ~Matrix(); friend double* Dooli(Matrix& ); void Input(); void Disp(); }; Matrix::Matrix(int x) { size=x; //为向量b分配空间并初始化为0 b=new double [x]; for(int j=0;j<x;j++) b[j]=0; //为向量A分配空间并初始化为0 A=new double* [x]; for(int i=0;i<x;i++) A[i]=new double [x]; for(int m=0;m<x;m++) for(int n=0;n<x;n++) A[m][n]=0; } Matrix::~Matrix() { cout<<"正在析构中~~~~"<<endl; delete b; for(int i=0;i<size;i++) delete A[i]; delete A; } void Matrix::Disp() { for(int i=0;i<size;i++) { for(int j=0;j<size;j++) cout<<A[i][j]<<" "; cout<<endl; } } void Matrix::Input() { cout<<"请输入A:"<<endl; for(int i=0;i<size;i++) for(int j=0;j<size;j++){ cout<<"第"<<i+1<<"行"<<"第"<<j+1<<"列:"<<endl; cin>>A[i][j]; } cout<<"请输入b:"<<endl; for(int j=0;j<size;j++){ cout<<"第"<<j+1<<"个:"<<endl; cin>>b[j]; } } double* Dooli(Matrix& A) { double *Xn=new double [A.size]; Matrix L(A.size),U(A.size); //分别求得U,L的第一行与第一列 for(int i=0;i<A.size;i++) U.A[0][i]=A.A[0][i]; for(int j=1;j<A.size;j++) L.A[j][0]=A.A[j][0]/U.A[0][0]; //分别求得U,L的第r行,第r列 double temp1=0,temp2=0; for(int r=1;r<A.size;r++){ //U for(int i=r;i<A.size;i++){ for(int k=0;k<r-1;k++) temp1=temp1+L.A[r][k]*U.A[k][i]; U.A[r][i]=A.A[r][i]-temp1; } //L for(int i=r+1;i<A.size;i++){ for(int k=0;k<r-1;k++) temp2=temp2+L.A[i][k]*U.A[k][r]; L.A[i][r]=(A.A[i][r]-temp2)/U.A[r][r]; } } cout<<"计算U得:"<<endl; U.Disp(); cout<<"计算L的:"<<endl; L.Disp(); double *Y=new double [A.size]; Y[0]=A.b[0]; for(int i=1;i<A.size;i++ ){ double temp3=0; for(int k=0;k<i-1;k++) temp3=temp3+L.A[i][k]*Y[k]; Y[i]=A.b[i]-temp3; } Xn[A.size-1]=Y[A.size-1]/U.A[A.size-1][A.size-1]; for(int i=A.size-1;i>=0;i--){ double temp4=0; for(int k=i+1;k<A.size;k++) temp4=temp4+U.A[i][k]*Xn[k]; Xn[i]=(Y[i]-temp4)/U.A[i][i]; } return Xn; } int main() { Matrix B(4); B.Input(); double *X; X=Dooli(B); cout<<"~~~~解得:"<<endl; for(int i=0;i<B.size;i++) cout<<"X["<<i<<"]:"<<X[i]<<" "; cout<<endl<<"呵呵呵呵呵"; return 0; }
标签: 道理特分解法
上传时间: 2018-05-20
上传用户:Aa123456789
VIP专区-嵌入式/单片机编程源码精选合集系列(60)资源包含以下内容:1. 16 * 16 的127个ascii码的点阵字库.2. ARM MP3解码源代码 实现MP3播放应用.3. 本文介基于CPLD和USB的多路温度数据采集系统.4. 89c51+sj1000 的发送例程 can2.0 扩展帧.5. 一个很好的串口类.6. cpld状态及设计。 很好的文章。 要设计vhdl状态机的话.7. c51formenu.c该文件是一个用c51编写一个menu的例子..8. MSP430F449 子程序,MSP-FET430P440 Demo - USART1 UART 38400 Echo ISR, DCO SMCLK..9. MSP-FET430P440 Demo - WDT Toggle P5.1 Interval overflow ISR, DCO SMCLK.10. MSP-FET430P440 Demo - WDT Toggle P5.1 Interval overflow ISR, 32kHz ACLK.11. MSP-FET430P440 Demo - BasicTimer Toggle P5.1 using ISR, DCO SMCLK.12. MSP-FET430P440 Demo - BasicTimer Toggle P5.1 using ISR, 32kHz SMCLK.13. MSP-FET430P440 Demo - Timer_A PWM TA1-2 upmode, DCO SMCLK.14. 恒颐arm9原理图(protel),直接可用.15. protel中CPLD器件的库可以方便的放进protel中.16. 这是一个以太网接口RTL8019AS和电路图一份.希望对大家有点参考.17. 关于lpc2292外接ds1302时钟芯片的驱动程序.18. 本文分别以GPI0口直接连接、串并转换连接、CPLD分部连接三种方法阐述了无外部总线的Philips ARM微控制器LPC2l0X与点阵图形液晶显示器的接口设计,并给出了硬件电路框图和主要程序。.19. arm 板的原理图 是三星的.20. BSP-15:高性能媒体处理DSP芯片开放指南说明文档.21. 如果您刚接触嵌入式开发.22. TMS320LF2407A原理图.23. 嵌入式arm-linux 2.6内核 s3c2410 的i2c总线驱动程序.24. 输入两组数.25. 应用VHDL、CPLD、EDA开发软件设计数字系统,能够显著增强设计的灵活性,提高产品的性能,减轻设计的工作量,缩短设计周期。传统的“固定功能集成块+连线”的设计方法正逐步地缩小应用范围,而基于芯片的.26. 使用lpc2114arm7内核的两个工程实例,采用ads编译器进行编译.27. 嵌入式系统开发的中英文词汇表并有详细注解,具有参考价值.28. 本程序仅供广大电子爱好者制作MP3学习和参考使用.29. arm-linux下.30. 适用于ucosii操作系统的文件系统.31. LCD 控制源代码.32. 用java写的GUI Analog Clock,用上了Observer,可以更改时间.33. s3c2410平台 直流电机有用的demo!大家共享.34. s3c2410平台 步进电机的demo.35. s3c2410平台can总线的一个例子.36. 音乐程序主要是控制频率而达到效果.37. 240x128 graphic mode/Driver ic:t6963C writed by hawk liu 2001 Aug. 22 1/128 duty 1/12.3 bia.38. USB芯片CH375EVT 资料源程序.39. 51 模拟SPI 51 模拟SPI 51 模拟SPI 51 模拟SPI.40. 基于时钟芯片DSC12C887的实时时钟系统.
上传时间: 2013-04-15
上传用户:eeworm
VIP专区-嵌入式/单片机编程源码精选合集系列(73)资源包含以下内容:1. 定标电压表程序.2. 电路原理图 电路原理图 电路原理图 电路原理图.3. 红外线检测模块.4. 一本老师推荐的经典的VHDL覆盖基础的入门书籍.5. 一本有关FPGA设计的经验的汇总.6. 下载的一些重要的电子设计文档.7. 搜集过来的相关cpld 的书籍 比较经典的几篇文章.8. 2410原理图1 2410原理图1 2410原理图1 2410原理图1.9. CH452驱动程序.10. s3c2440开发板的wince源代码完全无需修改应用在开发板上.11. s3c2440开发板的usb驱动程序经过测试完全适用.12. NIOSII的7个c语言源码.13. cheertek952原理图和PCB图跟BOM清单.14. 简单的m序列产生办法(存储m序列.15. 一种任意方波信号的产生方法.16. cpld 与8051的总线接口VHDL源码 非常好用 试过的.17. altera的ip核.18. altera的ip核.19. altera的ip核.20. 控制器为NT7532的LCD-cgg128064M02显示屏的示例程序.21. 压缩包中含有DS10B20,JM12864液晶模块的资料,以及DS18B201有LCD显示温度的2个KEIL C程序.和电路原理图..22. TI公司的zigbee的文档说明.23. 单总线和i2c总线结合实现数字温度计的程序.24. IC卡身份识别开关的程序.25. 步进电机控制的程序.26. U盘电路原理图.27. LINUX44B0I2C驱动验证过 现在可以用 是中断方式的驱动.28. Atmel算法(pid)-步进电机驱动源码 HOWTO,不好用你找我,绝对ok!.29. 12bitAD-ThermocoupleSensorIC max6675(工业级标准) 不好用你找我.30. U盘和MP3上用的Fat文件系统源码.31. 精确定时+中断+进制转换+移位+LED显示+按关键加速+EEPROM 存储测试.32. CY7C63722的全部设计资料.33. Cygnal 8051F000 I2C Example Program.34. 一些常用串口的例子.35. PS2Key.c键盘扫描+液晶1602显示.36. 很好的51显示程序!做为一个人机对话,HAO.37. 完整的Linux 设备驱动第3版.38. 232同485转换的原理图.39. c51 单片红外通信 c51 单片红外通信.40. ALTER的CPLD/FPGA配置芯片的手册.
上传时间: 2013-07-22
上传用户:eeworm
非均匀有理b样条曲线nurbs的C++源码库。 3.0.11版
上传时间: 2015-04-24
上传用户:Miyuki
Haskell 中文教程1-3章(flw译) 纯粹的函数型编程语言。以著名逻辑学家 Haskell B. Curry 的名字命名。 最初的目的是想要设计出满足下面这些要求的语言: 1,它必须能够适合教学、研究、应用开发,包括一些大系统的构造。 2,它必须能够使用形式语言来准确描述。 3,它必须是自由免费的,任何人只要愿意都可以获取、使用和再次发布它。 4,它必须建立在大家一致认可的基础上。 5,它应该能够消除目前的函数型编程语言的差异。
上传时间: 2014-01-10
上传用户:685
Delphi2005 BS程序设计技巧集 (1-3) 现在delphi 2005 B/S的书籍资料太少,几乎没有可参考的资料,这在一定程度上限制了delphi 2005的使用,相反C#,ASP的书籍资料到处都是。通过几个月的学习也算是总结了一点经验,从现在起我将陆续将ASP和C#的例子、技巧翻译到delphi 2005下面,希望对大家的学习和工作有帮助,毕竟.Net是一种发展趋势。我将不定期写新的内容,同时,如果你在工作或学习中有什么问题,我也会将问题加入本篇文章,另外,一个人的能力毕竟有限,也希望大家共同来解决问题。 宋雨炫
上传时间: 2013-11-25
上传用户:lgnf
车牌定位---VC++源代码程序 1.24位真彩色->256色灰度图。 2.预处理:中值滤波。 3.二值化:用一个初始阈值T对图像A进行二值化得到二值化图像B。 初始阈值T的确定方法是:选择阈值T=Gmax-(Gmax-Gmin)/3,Gmax和Gmin分别是最高、最低灰度值。 该阈值对不同牌照有一定的适应性,能够保证背景基本被置为0,以突出牌照区域。 4.削弱背景干扰。对图像B做简单的相邻像素灰度值相减,得到新的图像G,即Gi,j=|Pi,j-Pi,j-1|i=0,1,…,439 j=0,1,…,639Gi,0=Pi,0,左边缘直接赋值,不会影响整体效果。 5.用自定义模板进行中值滤波 区域灰度基本被赋值为0。考虑到文字是由许多短竖线组成,而背景噪声有一大部分是孤立噪声,用模板(1,1,1,1,1)T对G进行中值滤波,能够得到除掉了大部分干扰的图像C。 6.牌照搜索:利用水平投影法检测车牌水平位置,利用垂直投影法检测车牌垂直位置。 7.区域裁剪,截取车牌图像。
上传时间: 2013-11-26
上传用户:懒龙1988
s平面中直接形式到级联形式的转换 %适合模拟滤波器的 %C为增益系数 %B为包含各bk的K乘3维实系数矩阵 %A为包含各ak的K乘3维实系数矩阵 %b为直接形式的分子多项式系数 %a为直接形式的分母多项式系数
上传时间: 2015-07-22
上传用户:sdq_123
%直接型到并联型的转换 % %[C,B,A]=dir2par(b,a) %C为当b的长度大于a时的多项式部分 %B为包含各bk的K乘2维实系数矩阵 %A为包含各ak的K乘3维实系数矩阵 %b为直接型分子多项式系数 %a为直接型分母多项式系数 %
上传时间: 2014-01-20
上传用户:lizhen9880