r语言自定义函数对matrix()函数的调用

怎么在函数里调用matrix,可以实现我想要的矩阵,并为两个矩阵的代数积做准备。

1个回答

Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!
其他相关推荐
Eigen在调用函数里面的问题
C++中声明一些矩阵时,使用eigen很方便,而且eigen里面的一些函数用起来也很方便。我想问一下对于调用的一些函数可以使用eigen声明吗?像可以这样写:double a[3],b[3][3];然后调用函数ma(a,b)。但是声明如果换成Vector,Matrix,用eigen的话,后面调用再这么写的话就会出错,求大神指点,能否给个实例程序!
OpenGL绘图颜色设置无效(MFC环境下)
用MFC做了一个小程序,大致功能是先显示一个读入的图形,然后在条件变化时显示另一个图形,最初的显示是正常的,然后消息处理函数也没问题,自定义了一个draw函数如下 void CReadView::Draw(){ glClearColor(1.0, 1.0,1.0, 0.0); glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT); // Clear Screen And Depth Buffer glEnable(GL_COLOR_MATERIAL); glLoadIdentity(); // Reset The Current Modelview Matrix glTranslatef(-1.5f, 0.0f, -6.0f);//物体左移1.5,向内移6,相当于移动镜头一样,让物体进入镜头中 glColor3f(1.0, 0.0, 0.0); glBegin(GL_TRIANGLES); // 绘制三角形(逆时针画出来的三角形才是正面朝着我们的) glVertex3f(-1.0, -1.0, 0.0); // 左下顶点(//glVertex 的第一个参数是X坐标,然后依次是Y坐标和Z坐标) glVertex3f(1.0, -1.0, 0.0); // 右下顶点 glVertex3f(0.0, 1.0, 0.0); // 上顶点 glEnd(); // 三角形绘制结束 glColor3f(0.4, 0.75, 1.0); glRectf(1.5f, 1.5f, 2.5f, 2.5f); glLineWidth(3); glBegin(GL_LINES); glColor3f(1.0, 0.0, 0.0); glVertex2f(3.0f, 0.0f); glVertex2f(3.5f, 0.0f); //在这里使用glVertex*系列函数 指定你所希望的顶点位置 glEnd(); glLineWidth(1); glColor3f(1.0f, 1.0f, 0.0f); glBegin(GL_LINES); glVertex3f(-3.0f, 2.0f, 3.0f); glVertex3f(3.0f, 1.0f, 0.0f); glEnd(); glColor3f(1.0f,0.0f,0.0f); glBegin(GL_LINES); glVertex2i(1, 1); glVertex2i(3, 3); glEnd(); glFlush(); } 把它放到的OnDraw函数里进行调用,程序可以正常执行,只是画出来的这些线段、图形都没有颜色(是默认的黑色),也就是glColor3f似乎失效了。 这些绘制的函数都是从网上找来做测试的,在正常单独的OpenGL环境里应该也都是没问题的。是MFC的问题吗?或者是OpenGL需要某些重新的预设置?求解答
迷你opengl的函数调用来计算坐标,怎么用C语言的程序的编写设计的代码的方式来实现的
Problem Description Jiaoshou selected a course about “openGL” this semester. He was quite interested in modelview, which is a part of “openGL”. Just using three functions, it could make the model to move, rotate and largen or lessen. But he was puzzled with the theory of the modelview. He didn’t know a vertex after several transformations where it will be. Now, He tells you the position of the vertex and the transformations. Please help Jiaoshou find the position of the vertex after several transformations. Input The input will start with a line giving the number of test cases, T. Each case will always begin with “glBegin(GL_POINTS);”.Then the case will be followed by 5 kinds of function. 1. glTranslatef(x,y,z); This function will translate the vertex(x’,y’,z’) to vertex(x+x’,y+y’,z+z’). 2. glRotatef(angle,x,y,z); This function will turn angle radians counterclockwise around the axis (0,0,0)->(x,y,z). 3. glScalef(x,y,z); This function wiil translate the vertex(x’,y’,z’) to vertex(x*x’,y*y’,z*z’). 4. glVertex3f(x,y,z); This function will draw an initial vertex at the position(x,y,z). It will only appear once in one case just before “glEnd();”. In openGL, the transformation matrices are right multiplied by vertex matrix. So you should do the transformations in the reverse order. 5. glEnd(); This function tells you the end of the case. In this problem angle,x,y,z are real numbers and range from -50.0 to 50.0. And the number of functions in each case will not exceed 100. Output For each case, please output the position of the vertex after several transformations x,y,z in one line rounded to 1 digits after the decimal point , separated with a single space. We guarantee that x,y,z are not very large. Sample Input 1 glBegin(GL_POINTS); glScalef(2.0,0.5,3.0); glTranslatef(0.0,1.0,0.0); glVertex3f(1.0,1.0,1.0); glEnd(); Sample Output 2.0 1.0 3.0
tf.matmul函数的用法请教
我查到tf.matmul函数的原型是: ```tf.matmul (matrix1,matrix2, transpose_matrix1=False, transpose_matrix2=False, a_is_spare=False, b_is_spare=False, name=None ) ``` 如果在调用该函数时,把 transpose_a=True 有什么用?
析构函数中使用delete报错
把析构函数里的delete去掉就不会报错了;但是搞不懂:我的构造函数里不是用了new吗,为什么会出现这种情况? 报错信息是:_BLOCK_TYPE_IS_VALID(pHead->nBlockUse),网上查了一下都说是重发调用delete的结果,但是我好像没有那么做。 //Matrix.h class Matrix { public: Matrix(); ~Matrix(); Matrix(int,double*); //矩阵行列数、数组指针 Matrix operator+(const Matrix&);//矩阵相加 friend ostream& operator<<(ostream&,const Matrix&);// private: int n; //矩阵的行列数 double* pt; }; //Matrix.cpp Matrix::Matrix():pt(0) { n = 0; } Matrix::Matrix(int r,double* p): n(r) { pt = new double[r*r]; for(int i=0; i<r*r;i++) pt[i] = p[i]; } Matrix::~Matrix() { delete []pt; pt = 0; } //矩阵相加函数实现 Matrix Matrix::operator+(const Matrix& m) { for(int i=0;i<n*n;i++) pt[i] = pt[i] + m.pt[i]; return *this; } //矩阵输出函数的实现 ostream& operator<<(ostream& output,const Matrix& m) { for(int i=0;i<m.n*m.n;i++) { if((i+1)%m.n == 0) output<<m.pt[i]<<endl; else output<<m.pt[i]<<" "; } return output; } 然后就是main函数中 double *p1,*p2,x,y; int n; cout<<"请输两个矩阵的行列数(n*n):"<<endl<<"n="; cin>>n; cout<<"请输入第一个矩阵m1的元素:"<<endl; p1=(double *)malloc(sizeof(double)*n*n); for(int i=0;i<n*n;i++) { cin>>x; p1[i]=x; } cout<<"请输入第二个矩阵m2的元素:"<<endl; p2=(double*)malloc(n*n*sizeof(double)); for(int j=0;j<n*n;j++) { cin>>y; p2[j]=y; } Matrix m1(n,p1),m2(n,p2),M; M=m1+m2; cout<<"m1+m2:"<<endl; cout<<M<<endl; free(p1); free(p2); return 0;
C++中返回数组供函数调用
如下是我编写的高斯顺序消元法程序,input函数是我自定义的,希望输入一个矩阵,然后给Gauss和zg matrix函数调用,但是输入的矩阵不知如何传递过去。 #include <stdio.h> //包含输入输出函数 #include <stdlib.h>//包含动态内存分配函数。 #define n 3 float input()//输入矩阵函数************出问题的地方 { float **a; //定义二维指针。 float *b;//定义一维指针。 int i,j; a = (float **) malloc(sizeof(float *) *n);//申请一组一维指针空间。 for(i = 0; i < n; i ++) a[i] = (float *)malloc(sizeof(float) *n); //对于每个一维指针,申请一行数据的空间。 b = (float *)malloc(sizeof(float) *n); printf("请输入系数矩阵:\n"); for(i = 0; i < n; i ++) { for(j = 0; j < n; j ++) { printf("a[%d][%d]: ", i+1, j+1); scanf("%f",&a[i][j]);//输入第i行第j列的数据。 } }return a; printf("请输入右端项数组:\n"); for(i=0; i<n; i++) { printf("b[%d]: ", i+1); scanf("%f", &b[i]); }return b; //释放内存 for(i = 0; i < n; i ++) free(a[i]); free(a);free(b); } int Gauss(float a[n][n],float b[n])//高斯消元过程 { int i,j,k,flag=1; float t; for(i=0;i<n-1;i++) { if(a[i][i]==0) { flag=0; break; } else { for(j=i+1;j<n;j++) { t=-a[j][i]/a[i][i]; b[j]=b[j]+t*b[i]; for(k=i;k<n;k++) a[j][k]=a[j][k]+t*a[i][k]; } } } return (flag); } void zg_matrix(float a[n][n],float b[n])//输出增广矩阵 { int i,j; for(i=0;i<n;i++) { for(j=0;j<n;j++) { printf("%10f",a[i][j]); } printf(" | "); printf("%10f\n",b[i]); } } void main() { float input(); float x[n]={0}; int i,j,flag; zg_matrix(a,b); flag=Gauss(a,b); if(flag==0) printf("Gauss Method Fail!!!\n");//求解 else { x[n-1]=b[n-1]/a[n-1][n-1]; for(i=n-2;i>=0;i--) { x[i]=b[i]; for(j=i+1;j<n;j++) x[i]=x[i]-a[i][j]*x[j]; x[i]=x[i]/a[i][i]; } for(i=0;i<n;i++) printf("x%d=%11.7f\n",i+1,x[i]); } }
android Bitmap的内存溢出,何处添加回收函数?
public Bitmap rotateBitmap(int degree, Bitmap bitmap) { Matrix matrix = new Matrix(); matrix.postRotate(degree); Bitmap bm = Bitmap.createBitmap(bitmap, 0, 0, bitmap.getWidth(), bitmap.getHeight(), matrix, true); return bm; } 在上面的函数中,如果不调用bm.recyle(),那添加几张图片就报内存溢出了,但是,我不知道回收函数应该放在哪?要是放在return前面那返回的值填充ImageView就成黑色的了
请问调用时候参数写什么?求助大神!我是个小白!
class Matrix:public MATRIX<double> { public: Matrix():MATRIX<double>(){} Matrix( int c, int r ):MATRIX<double>(c,r){} Matrix( const Matrix& m){ *this = m; } const Matrix& operator+=( const Matrix& m ); const Matrix& operator-=( const Matrix& m ); const Matrix& operator*=( const Matrix& m ); const Matrix& operator/=( const Matrix& m ); }; bool operator==( const Matrix& lhs, const Matrix& rhs ); // 重载操作符== bool operator!=( const Matrix& lhs, const Matrix& rhs ); // 重载操作符!= const Matrix operator+( const Matrix& lhs, const Matrix& rhs ); // 重载操作符+ const Matrix operator-( const Matrix& lhs, const Matrix& rhs ); // 重载操作符- const Matrix operator*( const Matrix& lhs, const Matrix& rhs ); // 重载操作符* const Matrix operator/( const Matrix& lhs, const Matrix& rhs ); // 重载操作符/ const double det( const Matrix& m ); // 计算行列式 const double det( const Matrix& m, int start, int end ); // 计算子矩阵行列式 const Matrix abs( const Matrix& m ); // 计算所有元素的绝对值 const double max( const Matrix& m ); // 所有元素的最大值 const double max( const Matrix& m, int& row, int& col); // 所有元素中的最大值及其下标 const double min( const Matrix& m ); // 所有元素的最小值 const double min( const Matrix& m, int& row, int& col); // 所有元素的最小值及其下标 const Matrix trans( const Matrix& m ); // 返回转置矩阵 const Matrix submatrix(const Matrix& m,int rb,int re,int cb,int ce); // 返回子矩阵 const Matrix inverse( const Matrix& m ); // 计算逆矩阵 const Matrix LU( const Matrix& m ); // 计算方阵的LU分解 const Matrix readMatrix( istream& in = std::cin ); // 从指定输入流读入矩阵 const Matrix readMatrix( string file ); // 从文本文件读入矩阵 const Matrix loadMatrix( string file ); // 从二进制文件读取矩阵 void printMatrix( const Matrix& m, ostream& out = std::cout ); // 从指定输出流打印矩阵 void printMatrix( const Matrix& m, string file); // 将矩阵输出到文本文件 void saveMatrix( const Matrix& m, string file); // 将矩阵保存为二进制文件 #endif 这是.h文件我想在另一个.cpp文件中调用const Matrix readMatrix( istream& in = std::cin ); 这个函数,不知道怎么写?麻烦大神帮我一下,我真是个小白。。
定义一个矩阵类(采用顺序存储)
定义一个矩阵类(采用顺序存储),类中除了定义二维数组(最大50*50)用于存放矩阵元素,同时存放矩阵的行列值。 为该类重载矩阵的加法运算和乘法运算,以及初始化函数和show函数,可根据需要定义其他函数。 在主函数中应用该类实现两个矩阵的数据值输入,检查矩阵的行列相等时,调用加法运算,并显示输出结果。 下面是我的代码 为什么运行错误 求大神讲解!!!! #include<iostream> using namespace std; class Matrix{ private: int mar[50][50]; int row,col; public: Matrix(); void Input(int r,int c); void show(); Matrix operator+(Matrix &); Matrix operator*(Matrix &); }; Matrix::Matrix() { int r,c; row=r;col=c; cin>>r>>c; for(int i=0;i<row;i++) for( int j=0;j<col;j++) mar[i][j]=0; } void Matrix::Input(int r,int c) { row=r;col=c; cin>>r>>c; for(int i=0;i<row;i++) for( int j=0;j<col;j++) cin>>mar[i][j]; } void Matrix::show() { cout<<"The sum of two matrices is:"<<endl; for(int i=0;i<row;i++) for( int j=0;j<col;j++) { cout<<mar[i][j]<<" "; cout<<endl; } } Matrix Matrix::operator+(Matrix &m) { int r,c; row=r;col=c; Matrix temp; for(int i=0;i<row;i++) for(int j=0;j<col;j++) temp.mar[i][j]=mar[i][j]+m.mar[i][j]; return temp; } Matrix Matrix:: operator*(Matrix &m) { int r,c; row=r;col=c; Matrix stemp; for(int i=0;i<row;i++) for(int j=0;j<col;j++) stemp.mar[i][j]=mar[i][j]*m.mar[i][j]; return stemp; } int main() { int r1,c1,r2,c2; Matrix a,b,c; a.Input(r1,c1); b.Input(r2,c2); if(r1==r2&&c1==c2) { c=a+b; c.show(); } else cout<<"Fail to add"<<endl; return 0; }
运行到广度优先遍历函数调用队列时就出问题
#include<iostream> using namespace std; /*********************各类定义******************************/ class Node//基本抽象数据类型 { public: char ch; //记录名称,如果将这里改成数组,结点名称可以是多个字符 int flag;//记录结点是否被访问 }; class Graph //图类,此类中,封装了图的一些成员和一些必须的成员函数 { private: int getSub(char); //获取某名称的下标 Node* arrNode; //记录名称和是否访问的数组 int numVertex,numEdge;//记录图的顶点数和边数 int **matrix; //用一个二维数组记录两点间是否相连,1相连,0断开 public: void setCh(char,int); //将数组的arrNode的每一个单元设置一个结点名称 Graph(int); ~Graph(); int getNumVertex();//获得图的顶点数 char first(char ch);//获得相邻结点 char next(char ch1,char ch2);//获得隔着ch2,但与ch2相邻的结点 void setEdge(char,char,int w=1);//设置两顶点的边和权重(权重默认为1) void eraserEdge(char,char); void outPut(); int getMark(char);//获取是否被访问的记录,已访问返回1,未访问返回0 void setMark(char);//把已访问的结点,设置标记 }; template<class T> class myQueue { private: int theFront; // 1 counterclockwise from theFront element int theBack; // position of theBack element int arrayLength; // queue capacity T *queue; // element array public: myQueue(int initialCapacity = 10); ~myQueue() {delete [] queue;} bool empty() const {return theFront == theBack;} int size() const {return (theBack - theFront + arrayLength) % arrayLength;} T& front() {// return front element /*if (theFront == theBack) throw queueEmpty();*/ return queue[(theFront + 1) % arrayLength]; } T& back() {// return theBack element /*if (theFront == theBack) throw queueEmpty();*/ return queue[theBack]; } void pop() {// remove theFront element /*if (theFront == theBack) throw queueEmpty();*/ theFront = (theFront + 1) % arrayLength; queue[theFront].~T(); // destructor for T } void push(const T& theElement); }; /*************各函数的函数体*****************************/ Graph::Graph(int n) { int i,j; arrNode=new Node[n]; numVertex=n; numEdge=0; for(i=0;i<numVertex;i++) arrNode[i].flag=0; matrix=new int*[numVertex]; for(i=0;i<numVertex;i++) matrix[i]=new int[numVertex]; for(i=0;i<numVertex;i++) for(j=0;j<numVertex;j++) matrix[i][j]=matrix[j][i]=0; } Graph::~Graph() { delete [] arrNode; for(int i=0;i<numVertex;i++) delete [] matrix[i]; delete [] matrix; } int Graph::getSub(char ch) //获取下标 { for(int i=0;i<numVertex;i++) if(ch==arrNode[i].ch) return i; } int Graph::getNumVertex() //获取顶点数 { return numVertex; } char Graph::first(char ch) //获取相邻结点 { for(int i=0;i<numVertex;i++) if(matrix[getSub(ch)][i]!=0) return arrNode[i].ch; return ch; } char Graph::next(char ch1,char ch2) //获取与ch2相邻的结点 { for(int i=getSub(ch2)+1;i<numVertex;i++) if(matrix[getSub(ch1)][i]!=0) return arrNode[i].ch; return ch1; } void Graph::setEdge(char ch1,char ch2,int w) //设置边 { if(matrix[getSub(ch1)][getSub(ch2)]==0) { numEdge++; matrix[getSub(ch1)][getSub(ch2)]=matrix[getSub(ch2)][getSub(ch1)]=1; } } void Graph::eraserEdge(char ch1,char ch2) { if(matrix[getSub(ch1)][getSub(ch2)]==1) { numEdge--; matrix[getSub(ch1)][getSub(ch2)]=matrix[getSub(ch2)][getSub(ch1)]=0; } } void Graph::outPut() { for(int i=0;i<numVertex;i++) { for(int j=0;j<numVertex;j++) { cout<<matrix[i][j]<<" "; } cout<<endl<<endl;; } } int Graph::getMark(char ch)//获取访问信息 { for(int i=0;i<numVertex;i++) if(ch==arrNode[i].ch) return arrNode[i].flag; return -1; } void Graph::setMark(char ch) //设置标记 { for(int i=0;i<numVertex;i++) if(ch==arrNode[i].ch) arrNode[i].flag=1; } void Graph::setCh(char ch,int n)//将结点名称设置在字符中 { arrNode[n].ch=ch; } template<class T> myQueue<T>::myQueue(int initialCapacity) {// Constructor. /*if (initialCapacity < 1) {ostringstream s; s << "Initial capacity = " << initialCapacity << " Must be > 0"; throw illegalParameterValue(s.str()); }*/ arrayLength = initialCapacity; queue = new T[arrayLength]; theFront = 0; theBack = 0; } template<class T> void myQueue<T>::push(const T& theElement) {// Add theElement to queue. // increase array length if necessary if ((theBack + 1) % arrayLength == theFront) {// double array length // allocate a new array T* newQueue = new T[2 * arrayLength]; // copy elements into new array int start = (theFront + 1) % arrayLength; if (start < 2) // no wrap around copy(queue + start, queue + start + arrayLength - 1, newQueue); else { // queue wraps around copy(queue + start, queue + arrayLength, newQueue); copy(queue, queue + theBack + 1, newQueue + arrayLength - start); } // switch to newQueue and set theFront and theBack theFront = 2 * arrayLength - 1; theBack = arrayLength - 2; // queue size arrayLength - 1 arrayLength *= 2; queue = newQueue; } // put theElement at the theBack of the queue theBack = (theBack + 1) % arrayLength; queue[theBack] = theElement; } /***************主函数*********************/ void DFS(Graph*,char); //深度优先遍历函数声明 void BFS(Graph*,char,myQueue*); //广度优先遍历函数声明 int main() { char option; int numVer,numE,i; char temp1,temp2; //temp1,temp2作临时变量,记录输入的值 // myQueue* q=new myQueue; cout << "输入定点数和弧数:"; cin >> numVer >> numE; Graph myGraph1(numVer); Graph myGraph2(numVer); cout << "请输入" << numVer << "个顶点:\n"; for(i=0;i<numVer;i++) //将结点名称设置在数组中 { cout << "输入顶点" << i << ":"; cin >> temp1; myGraph1.setCh(temp1,i); myGraph2.setCh(temp1,i); } cout << "请输入" << numE << "条弧:\n"; for(i=0;i<numE;i++) //设置边 { cout << "输入弧" << i << ":"; cin >> temp1 >> temp2; myGraph1.setEdge(temp1,temp2); myGraph2.setEdge(temp1,temp2); } myGraph1.outPut(); while(1) { cout<<"2:删除边 "<<endl; cout<<"3: 深度遍历"<<endl; cout<<"4: 广度遍历"<<endl; cout<<"0: 退出"<<endl; cout<<"\n请输入你的选择:"; cin>>option; switch(option) { case'0': exit(0); break; case'2': { cout << "输入弧的顶点:" << endl; cin >> temp1 >> temp2; myGraph1.eraserEdge(temp1,temp2); myGraph2.eraserEdge(temp1,temp2); myGraph1.outPut(); break; } case'3': { cout << "深度优先结果:"; DFS(&myGraph1,'a'); cout<<endl; break; } case'4': { cout << "\n广度优先结果:"; //BFS(&myGraph2,'a',q); cout<<endl; break; } cout<<endl; } }cout<<endl; } void DFS(Graph* G,char ch) //深度优先遍历函数体 { int i=0; cout << ch << " "; G->setMark(ch); for(char w=G->first(ch);i<G->getNumVertex();w=G->next(ch,w),i++) if(G->getMark(w)==0) DFS(G,w); } void BFS(Graph* G,char ch,myQueue* q) //广度优先遍历函数体 { char v,w; q->push(ch); G->setMark(ch); while(q->size()!=0) { q->pop(v); cout << v << " "; int i=0; for(w=G->first(v);i<G->getNumVertex();w=G->next(v,w),i++) if(G->getMark(w)==0) { G->setMark(w); q->push(w); } } } /*输入顶点数和弧数:8 9 输入8个顶点. 输入顶点0:a 输入顶点1:b 输入顶点2:c 输入顶点3:d 输入顶点4:e 输入顶点5:f 输入顶点6:g 输入顶点7:h 输入9条弧. 输入弧0:a b 1 输入弧1:b d 1 输入弧2:b e 1 输入弧3:d h 1 输入弧4:e h 1 输入弧5:a c 1 输入弧6:c f 1 输入弧7:c g 1 输入弧8:f g 1 */
有点奇怪的矩阵问题,C++的哦,快来位大佬
题目描述 编写以下函数: (1)在一个二维数组中形成以下形式的n阶矩阵: [1 1 1 1 1 2 1 1 1 1 3 2 1 1 1 4 3 2 1 1 5 4 3 2 1] (2)去掉靠边的元素,生成新的n-2阶矩阵; (3)求生成的n阶矩阵主对角线上的元素之和; (4)以方阵形式输出数组。 在main函数中调用以上函数进行测试。 输入 输入生成矩阵的阶数(n>=2) 输出 以方阵形式输出生成的n阶矩阵、去掉靠边的元素生成的新的n-2阶矩阵、以及生成的n阶矩阵主对角线上的元素之和,最后一行要回车 样例输入 5 样例输出 Generated matrix: 1 1 1 1 1 2 1 1 1 1 3 2 1 1 1 4 3 2 1 1 5 4 3 2 1 del the elements on the side: 1 1 1 2 1 1 3 2 1 The sum of the diagonal:5
帮忙看看这个程序出了什么问题
本菜写了一个小程序 能够创建一个矩阵,规定长宽,并进行矩阵之间的加减运算 代码如下 ``` #include<iostream> using namespace std; class matrix //创建矩阵类 { private: int lines; int rows; int **num; //建立二重数组 public: matrix(int a,int b) //构造函数 { lines=a; rows=b; cout<<"矩阵初始化"<<endl; num=new int*[lines]; int x,y; for(y=0;y<lines;y++) { num[y]=new int[rows]; //循环完成new申请内存 } } /*matrix(matrix& father) //拷贝构造函数 { cout<<"调用拷贝构造函数"<<endl; lines=father.lines; rows=father.rows; int x,y; for(y=0;y<lines;y++) { for(x=0;x<rows;x++) { num[x][y]=father.num[x][y]; //循环读入矩阵元素 } } }*/ void in(); //输入函数 void out(); //输出函数 int get_num(int x,int y) //取用指定元素的函数(内联) { return num[x][y]; } void operate(int mode,matrix m1,matrix m2); //运算函数,包含加和减 ~matrix() //析构函数 释放矩阵空间 { } void matrix_free(); //释放矩阵空间函数 }; void matrix::matrix_free() { int i; for(i=0;i<lines;i++) { delete []num[i]; } delete []num; } void matrix::in() { cout<<"请按顺序输入矩阵元素:"<<endl; int x,y; for(y=0;y<lines;y++) { for(x=0;x<rows;x++) { cin>>num[x][y]; //循环输入矩阵元素 } } } void matrix::out() //此函数可以按格式输出矩阵 { int x,y; for(y=0;y<lines;y++) { for(x=0;x<rows;x++) { cout<<num[x][y]; if(x!=rows-1) { cout<<" "; } } cout<<endl; } } void matrix::operate(int mode,matrix m1,matrix m2) //运算函数,mode=1代表加,mode=-1代表减,m1,m2为待处理矩阵 { int x,y; for(y=0;y<lines;y++) { for(x=0;x<rows;x++) { num[x][y]=m1.get_num(x,y)+mode*m2.get_num(x,y); //运算部分 } } } int main(void) { int a=0,b=0; cout<<"输入矩阵大小(行 列):"<<endl; while(a<=0||b<=0) { cin>>a>>b; if(a<=0||b<=0) { cout<<"输入错误,重新输入:"<<endl; } } matrix A1(a,b),A2(a,b),A3(a,b); //创建三个矩阵 A1.in(); //输入初始化 A2.in(); A3.operate(1,A1,A2); //A3=A1+A2 cout<<"矩阵和:"<<endl; A3.out(); //输出 A3.operate(-1,A1,A2); //A3=A1-A2 cout<<"矩阵差:"<<endl; A3.out(); //输出 A1.matrix_free(); //释放矩阵空间 A2.matrix_free(); A3.matrix_free(); } ``` 此程序目前问题有三个 1.拷贝构造函数出错,只能注释掉,否则输入完两个矩阵后回车程序未响应 2.释放矩阵空间的不能写在析构函数里,否则输出结果不正确(A3=A1-A2的最后一行) 3.输入矩阵的列不能大于行,否则输入完第一行数字并回车后程序未响应 请前辈指导,谢谢
求问大神此程序哪里有问题
报错显示 调用 "load_plaintext" 时,未对输出参数 "plaintext_matrix" (可能还包括其他参数)赋值 代码如下 ``` key_dir = 'D:/Users/Fault/aes/AES/'; keyfile_name = 'key.mat'; Kname = sprintf('%s\\%s', key_dir, keyfile_name); load(Kname); plaintext_dir = 'D:/Users/Fault/heihei/'; plaintext_name = 'plaintext'; plaintext_length = 128; amount = 5000; unit_length = 8; plaintext_matrix = load_plaintext( plaintext_dir, plaintext_name, plaintext_length, amount, unit_length ); plaintext_hex = P; plaintext_hex = {'00' '11' '22' '33' '44' '55' '66' '77' ... '88' '99' 'aa' 'bb' 'cc' 'dd' 'ee' 'ff'}; [s_box, inv_s_box, w, poly_mat, inv_poly_mat] = aes_init(key, 14); num_P = size(P,1); plaintext = cell(num_P,1); for i = 1:num_P plaintext{i,:} = (hex2dec (plaintext_hex(i,:)))'; end Nr = 14; inductor = round(rand(1,1+(Nr-1)*4+3)); x = inductor(1); fault_point = randi(3); selector = randi(1+(Nr-1)*4+3); bitmap= inductor; for i =1:fault_point bitmap(selector)=1; bitmap(); end for i = 1:num_P fault = 0; C_normal{i,:} = cipher (plaintext{i,:}, w, s_box, poly_mat, 14, fault); fault = 1; C_fault{i,:} = cipher (plaintext{i,:}, w, s_box, poly_mat, 14, fault); T_base_start{i,:} = bitxor(C_normal{i,:}, C_fault{i,:}); end w_base = w; w_base(57:60,:) = 0; for i = 1:num_P [P_normal, normal_base] = inv_cipher (C_normal{i,:}, w_base, inv_s_box, inv_poly_mat, 14); inter_normal_base{i,:} = normal_base; [P_fault, fault_base] = inv_cipher (C_fault{i,:}, w_base, inv_s_box, inv_poly_mat, 14); inter_fault_base{i,:} = fault_base; T_base_end{i,:} = bitxor(inter_normal_base{i,:}, inter_fault_base{i,:}); [P_normal, inter_normal] = inv_cipher (C_normal{i,:}, w, inv_s_box, inv_poly_mat, 14); T_k_normal{i,:} = bitxor(inter_normal, inter_normal_base{i,:}); [P_normal, inter_fault] = inv_cipher (C_fault{i,:}, w, inv_s_box, inv_poly_mat, 14); T_k_fault{i,:} = bitxor(inter_fault, inter_fault_base{i,:}); end for j = 1:num_P for i = 1:pow2(8) guess = i - 1; w_guess = w_base; w_guess(57,1) = guess; [P_normal, inter_normal_guess] = inv_cipher (C_normal{j,:}, w_guess, inv_s_box, inv_poly_mat, 14); T_k_normal_guess{i,j} = bitxor(inter_normal_guess, inter_normal_base{j,:}); [P_normal, inter_fault_guess] = inv_cipher (C_fault{j,:}, w_guess, inv_s_box, inv_poly_mat, 14); T_k_fault_guess{i,j} = bitxor(inter_fault_guess, inter_fault_base{j,:}); end end for j = 1:num_P for i = 1:pow2(8) temp = bitxor(bitxor(T_k_normal{j,:},T_k_normal_guess{i,j}),bitxor(T_k_fault{j,:},T_k_fault_guess{i,j})); T_k_diff(i,j) = temp(1,1); end end product = 0; Q = zeros(1,num_P); for j = 1:num_P tmp = T_k_diff(:,j); ind = find(tmp == 0); product = product + size(ind,1); Q(1,j) = product/num_P; end product = product/num_P; plot(Q); load('D:/Users/Fault/AES DATA/Q.mat'); plot(Q); ``` load_plaintext代码如下 ``` %UNTITLED4 此处显示有关此函数的摘要load_plaintext % index_dir: 索引文件所在目录 % name:索引文件名 % format: 索引文件格式 % -- key_length:密钥比特长度 % -- plaintext_length:明文比特长度 % -- offset_length:偏移量比特长度 % amount: 所需的明文数据量 %UNTITLED Summary of this function goes here % Detailed explanation goes here if ~(exist('index_dir', 'var')) error('The directory of index files must be defined'); end if ~(exist('name', 'var')) error('The name of index files must be defined'); end if ~(exist('length', 'var')) error('The length of index must be defined'); end if ~(exist('amount', 'var')) error('The amount of index must be defined'); end if exist(index_dir, 'dir') index_file = sprintf('%s%s.txt',index_dir, name); if exist(index_file, 'file') fid = fopen(index_file); %tmp = textscan(fid, '%s %s %s %s %s %s'); %plaintext_origin = tmp{:,2}; tmp = textscan(fid, '%s'); plaintext_origin = tmp{:,1}; clear tmp; num = size(plaintext_origin,1); for i = 1:num for j = 1:(length/unit_length) A = plaintext_origin{i,:}; plaintext(i,j) = {A(j*2-1:j*2)}; end end if amount == 'all' plaintext_matrix = plaintext; else plaintext_matrix = plaintext(1:amount,:); end end else error('index file is not exist'); end end ``` ``` ```
Java学习的正确打开方式
在博主认为,对于入门级学习java的最佳学习方法莫过于视频+博客+书籍+总结,前三者博主将淋漓尽致地挥毫于这篇博客文章中,至于总结在于个人,实际上越到后面你会发现学习的最好方式就是阅读参考官方文档其次就是国内的书籍,博客次之,这又是一个层次了,这里暂时不提后面再谈。博主将为各位入门java保驾护航,各位只管冲鸭!!!上天是公平的,只要不辜负时间,时间自然不会辜负你。 何谓学习?博主所理解的学习,它是一个过程,是一个不断累积、不断沉淀、不断总结、善于传达自己的个人见解以及乐于分享的过程。
程序员必须掌握的核心算法有哪些?
由于我之前一直强调数据结构以及算法学习的重要性,所以就有一些读者经常问我,数据结构与算法应该要学习到哪个程度呢?,说实话,这个问题我不知道要怎么回答你,主要取决于你想学习到哪些程度,不过针对这个问题,我稍微总结一下我学过的算法知识点,以及我觉得值得学习的算法。这些算法与数据结构的学习大多数是零散的,并没有一本把他们全部覆盖的书籍。下面是我觉得值得学习的一些算法以及数据结构,当然,我也会整理一些看过
大学四年自学走来,这些私藏的实用工具/学习网站我贡献出来了
大学四年,看课本是不可能一直看课本的了,对于学习,特别是自学,善于搜索网上的一些资源来辅助,还是非常有必要的,下面我就把这几年私藏的各种资源,网站贡献出来给你们。主要有:电子书搜索、实用工具、在线视频学习网站、非视频学习网站、软件下载、面试/求职必备网站。 注意:文中提到的所有资源,文末我都给你整理好了,你们只管拿去,如果觉得不错,转发、分享就是最大的支持了。 一、电子书搜索 对于大部分程序员...
linux系列之常用运维命令整理笔录
本博客记录工作中需要的linux运维命令,大学时候开始接触linux,会一些基本操作,可是都没有整理起来,加上是做开发,不做运维,有些命令忘记了,所以现在整理成博客,当然vi,文件操作等就不介绍了,慢慢积累一些其它拓展的命令,博客不定时更新 free -m 其中:m表示兆,也可以用g,注意都要小写 Men:表示物理内存统计 total:表示物理内存总数(total=used+free) use...
比特币原理详解
一、什么是比特币 比特币是一种电子货币,是一种基于密码学的货币,在2008年11月1日由中本聪发表比特币白皮书,文中提出了一种去中心化的电子记账系统,我们平时的电子现金是银行来记账,因为银行的背后是国家信用。去中心化电子记账系统是参与者共同记账。比特币可以防止主权危机、信用风险。其好处不多做赘述,这一层面介绍的文章很多,本文主要从更深层的技术原理角度进行介绍。 二、问题引入 假设现有4个人...
程序员接私活怎样防止做完了不给钱?
首先跟大家说明一点,我们做 IT 类的外包开发,是非标品开发,所以很有可能在开发过程中会有这样那样的需求修改,而这种需求修改很容易造成扯皮,进而影响到费用支付,甚至出现做完了项目收不到钱的情况。 那么,怎么保证自己的薪酬安全呢? 我们在开工前,一定要做好一些证据方面的准备(也就是“讨薪”的理论依据),这其中最重要的就是需求文档和验收标准。一定要让需求方提供这两个文档资料作为开发的基础。之后开发...
网页实现一个简单的音乐播放器(大佬别看。(⊙﹏⊙))
今天闲着无事,就想写点东西。然后听了下歌,就打算写个播放器。 于是乎用h5 audio的加上js简单的播放器完工了。 演示地点演示 html代码如下` music 这个年纪 七月的风 音乐 ` 然后就是css`*{ margin: 0; padding: 0; text-decoration: none; list-...
Python十大装B语法
Python 是一种代表简单思想的语言,其语法相对简单,很容易上手。不过,如果就此小视 Python 语法的精妙和深邃,那就大错特错了。本文精心筛选了最能展现 Python 语法之精妙的十个知识点,并附上详细的实例代码。如能在实战中融会贯通、灵活使用,必将使代码更为精炼、高效,同时也会极大提升代码B格,使之看上去更老练,读起来更优雅。
数据库优化 - SQL优化
以实际SQL入手,带你一步一步走上SQL优化之路!
通俗易懂地给女朋友讲:线程池的内部原理
餐盘在灯光的照耀下格外晶莹洁白,女朋友拿起红酒杯轻轻地抿了一小口,对我说:“经常听你说线程池,到底线程池到底是个什么原理?”
经典算法(5)杨辉三角
杨辉三角 是经典算法,这篇博客对它的算法思想进行了讲解,并有完整的代码实现。
使用 Docker 部署 Spring Boot 项目
Docker 技术发展为微服务落地提供了更加便利的环境,使用 Docker 部署 Spring Boot 其实非常简单,这篇文章我们就来简单学习下。首先构建一个简单的 S...
英特尔不为人知的 B 面
从 PC 时代至今,众人只知在 CPU、GPU、XPU、制程、工艺等战场中,英特尔在与同行硬件芯片制造商们的竞争中杀出重围,且在不断的成长进化中,成为全球知名的半导体公司。殊不知,在「刚硬」的背后,英特尔「柔性」的软件早已经做到了全方位的支持与支撑,并持续发挥独特的生态价值,推动产业合作共赢。 而对于这一不知人知的 B 面,很多人将其称之为英特尔隐形的翅膀,虽低调,但是影响力却不容小觑。 那么,在...
面试官:你连RESTful都不知道我怎么敢要你?
干货,2019 RESTful最贱实践
刷了几千道算法题,这些我私藏的刷题网站都在这里了!
遥想当年,机缘巧合入了 ACM 的坑,周边巨擘林立,从此过上了"天天被虐似死狗"的生活… 然而我是谁,我可是死狗中的战斗鸡,智力不够那刷题来凑,开始了夜以继日哼哧哼哧刷题的日子,从此"读题与提交齐飞, AC 与 WA 一色 ",我惊喜的发现被题虐既刺激又有快感,那一刻我泪流满面。这么好的事儿作为一个正直的人绝不能自己独享,经过激烈的颅内斗争,我决定把我私藏的十几个 T 的,阿不,十几个刷题网...
白话阿里巴巴Java开发手册高级篇
不久前,阿里巴巴发布了《阿里巴巴Java开发手册》,总结了阿里巴巴内部实际项目开发过程中开发人员应该遵守的研发流程规范,这些流程规范在一定程度上能够保证最终的项目交付质量,通过在时间中总结模式,并推广给广大开发人员,来避免研发人员在实践中容易犯的错误,确保最终在大规模协作的项目中达成既定目标。 无独有偶,笔者去年在公司里负责升级和制定研发流程、设计模板、设计标准、代码标准等规范,并在实际工作中进行...
SQL-小白最佳入门sql查询一
不要偷偷的查询我的个人资料,即使你再喜欢我,也不要这样,真的不好;
redis分布式锁,面试官请随便问,我都会
文章有点长并且绕,先来个图片缓冲下! 前言 现在的业务场景越来越复杂,使用的架构也就越来越复杂,分布式、高并发已经是业务要求的常态。像腾讯系的不少服务,还有CDN优化、异地多备份等处理。 说到分布式,就必然涉及到分布式锁的概念,如何保证不同机器不同线程的分布式锁同步呢? 实现要点 互斥性,同一时刻,智能有一个客户端持有锁。 防止死锁发生,如果持有锁的客户端崩溃没有主动释放锁,也要保证锁可以正常释...
项目中的if else太多了,该怎么重构?
介绍 最近跟着公司的大佬开发了一款IM系统,类似QQ和微信哈,就是聊天软件。我们有一部分业务逻辑是这样的 if (msgType = "文本") { // dosomething } else if(msgType = "图片") { // doshomething } else if(msgType = "视频") { // doshomething } else { // doshom...
Nginx 原理和架构
Nginx 是一个免费的,开源的,高性能的 HTTP 服务器和反向代理,以及 IMAP / POP3 代理服务器。Nginx 以其高性能,稳定性,丰富的功能,简单的配置和低资源消耗而闻名。 Nginx 的整体架构 Nginx 里有一个 master 进程和多个 worker 进程。master 进程并不处理网络请求,主要负责调度工作进程:加载配置、启动工作进程及非停升级。worker 进程负责处...
“狗屁不通文章生成器”登顶GitHub热榜,分分钟写出万字形式主义大作
一、垃圾文字生成器介绍 最近在浏览GitHub的时候,发现了这样一个骨骼清奇的雷人项目,而且热度还特别高。 项目中文名:狗屁不通文章生成器 项目英文名:BullshitGenerator 根据作者的介绍,他是偶尔需要一些中文文字用于GUI开发时测试文本渲染,因此开发了这个废话生成器。但由于生成的废话实在是太过富于哲理,所以最近已经被小伙伴们给玩坏了。 他的文风可能是这样的: 你发现,...
程序员:我终于知道post和get的区别
是一个老生常谈的话题,然而随着不断的学习,对于以前的认识有很多误区,所以还是需要不断地总结的,学而时习之,不亦说乎
《程序人生》系列-这个程序员只用了20行代码就拿了冠军
你知道的越多,你不知道的越多 点赞再看,养成习惯GitHub上已经开源https://github.com/JavaFamily,有一线大厂面试点脑图,欢迎Star和完善 前言 这一期不算《吊打面试官》系列的,所有没前言我直接开始。 絮叨 本来应该是没有这期的,看过我上期的小伙伴应该是知道的嘛,双十一比较忙嘛,要值班又要去帮忙拍摄年会的视频素材,还得搞个程序员一天的Vlog,还要写BU...
加快推动区块链技术和产业创新发展,2019可信区块链峰会在京召开
11月8日,由中国信息通信研究院、中国通信标准化协会、中国互联网协会、可信区块链推进计划联合主办,科技行者协办的2019可信区块链峰会将在北京悠唐皇冠假日酒店开幕。   区块链技术被认为是继蒸汽机、电力、互联网之后,下一代颠覆性的核心技术。如果说蒸汽机释放了人类的生产力,电力解决了人类基本的生活需求,互联网彻底改变了信息传递的方式,区块链作为构造信任的技术有重要的价值。   1...
Java世界最常用的工具类库
Apache Commons Apache Commons有很多子项目 Google Guava 参考博客
程序员把地府后台管理系统做出来了,还有3.0版本!12月7号最新消息:已在开发中有github地址
第一幕:缘起 听说阎王爷要做个生死簿后台管理系统,我们派去了一个程序员…… 996程序员做的梦: 第一场:团队招募 为了应对地府管理危机,阎王打算找“人”开发一套地府后台管理系统,于是就在地府总经办群中发了项目需求。 话说还是中国电信的信号好,地府都是满格,哈哈!!! 经常会有外行朋友问:看某网站做的不错,功能也简单,你帮忙做一下? 而这次,面对这样的需求,这个程序员...
网易云6亿用户音乐推荐算法
网易云音乐是音乐爱好者的集聚地,云音乐推荐系统致力于通过 AI 算法的落地,实现用户千人千面的个性化推荐,为用户带来不一样的听歌体验。 本次分享重点介绍 AI 算法在音乐推荐中的应用实践,以及在算法落地过程中遇到的挑战和解决方案。 将从如下两个部分展开: AI算法在音乐推荐中的应用 音乐场景下的 AI 思考 从 2013 年 4 月正式上线至今,网易云音乐平台持续提供着:乐屏社区、UGC...
【技巧总结】位运算装逼指南
位算法的效率有多快我就不说,不信你可以去用 10 亿个数据模拟一下,今天给大家讲一讲位运算的一些经典例子。不过,最重要的不是看懂了这些例子就好,而是要在以后多去运用位运算这些技巧,当然,采用位运算,也是可以装逼的,不信,你往下看。我会从最简单的讲起,一道比一道难度递增,不过居然是讲技巧,那么也不会太难,相信你分分钟看懂。 判断奇偶数 判断一个数是基于还是偶数,相信很多人都做过,一般的做法的代码如下...
为什么要学数据结构?
一、前言 在可视化化程序设计的今天,借助于集成开发环境可以很快地生成程序,程序设计不再是计算机专业人员的专利。很多人认为,只要掌握几种开发工具就可以成为编程高手,其实,这是一种误解。要想成为一个专业的开发人员,至少需要以下三个条件: 1) 能够熟练地选择和设计各种数据结构和算法 2) 至少要能够熟练地掌握一门程序设计语言 3) 熟知所涉及的相关应用领域的知识 其中,后两个条件比较容易实现,而第一个...
Android 9.0 init 启动流程
阅读五分钟,每日十点,和您一起终身学习,这里是程序员Android本篇文章主要介绍Android开发中的部分知识点,通过阅读本篇文章,您将收获以下内容:一、启动流程概述一、 启动流程概述Android启动流程跟Linux启动类似,大致分为如下五个阶段。1.开机上电,加载固化的ROM。2.加载BootLoader,拉起Android OS。3.加载Uboot,初始外设,引导Kernel启动等。...
8年经验面试官详解 Java 面试秘诀
作者 |胡书敏 责编 | 刘静 出品 | CSDN(ID:CSDNnews) 本人目前在一家知名外企担任架构师,而且最近八年来,在多家外企和互联网公司担任Java技术面试官,前后累计面试了有两三百位候选人。在本文里,就将结合本人的面试经验,针对Java初学者、Java初级开发和Java开发,给出若干准备简历和准备面试的建议。 Java程序员准备和投递简历的实...
面试官如何考察你的思维方式?
1.两种思维方式在求职面试中,经常会考察这种问题:北京有多少量特斯拉汽车?某胡同口的煎饼摊一年能卖出多少个煎饼?深圳有多少个产品经理?一辆公交车里能装下多少个乒乓球?一个正常成年人有多少根头发?这类估算问题,被称为费米问题,是以科学家费米命名的。为什么面试会问这种问题呢?这类问题能把两类人清楚地区分出来。一类是具有文科思维的人,擅长赞叹和模糊想象,它主要依靠的是人的第一反应和直觉,比如小孩...
前后端分离,我怎么就选择了 Spring Boot + Vue 技术栈?
前两天又有小伙伴私信松哥,问题还是职业规划,Java 技术栈路线这种,实际上对于这一类问题我经常不太敢回答,每个人的情况都不太一样,而小伙伴也很少详细介绍自己的情况,大都是一两句话就把问题抛出来了,啥情况都不了解,就要指出一个方向,这实在是太难了。 因此今天我想从我学习 Spring Boot + Vue 这套技术栈的角度,来和大家聊一聊没有人指导,我是如何一步一步建立起自己的技术体系的。 线上大...
17张图带你解析红黑树的原理!保证你能看懂!
二叉查找树 由于红黑树本质上就是一棵二叉查找树,所以在了解红黑树之前,咱们先来看下二叉查找树。 二叉查找树(Binary Search Tree),也称有序二叉树(ordered binary tree),排序二叉树(sorted binary tree),是指一棵空树或者具有下列性质的二叉树: 若任意结点的左子树不空,则左子树上所有结点的值均小于它的根结点的值; 若任意结点的...
so easy! 10行代码写个"狗屁不通"文章生成器
前几天,GitHub 有个开源项目特别火,只要输入标题就可以生成一篇长长的文章。 背后实现代码一定很复杂吧,里面一定有很多高深莫测的机器学习等复杂算法 不过,当我看了源代码之后 这程序不到50行 尽管我有多年的Python经验,但我竟然一时也没有看懂 当然啦,原作者也说了,这个代码也是在无聊中诞生的,平时撸码是不写中文变量名的, 中文...
知乎高赞:中国有什么拿得出手的开源软件产品?(整理自本人原创回答)
知乎高赞:中国有什么拿得出手的开源软件产品? 在知乎上,有个问题问“中国有什么拿得出手的开源软件产品(在 GitHub 等社区受欢迎度较好的)?” 事实上,还不少呢~ 本人于2019.7.6进行了较为全面的回答,对这些受欢迎的 Github 开源项目分类整理如下: 分布式计算、云平台相关工具类 1.SkyWalking,作者吴晟、刘浩杨 等等 仓库地址: apache/skywalking 更...
相关热词 c#选择结构应用基本算法 c# 收到udp包后回包 c#oracle 头文件 c# 序列化对象 自定义 c# tcp 心跳 c# ice连接服务端 c# md5 解密 c# 文字导航控件 c#注册dll文件 c#安装.net
立即提问