Introdution to 3D Game Programming With DirectX11 第2章 习题解答

8 .不是

9.不是

18.

// Win32Project1.cpp : Defines the entry point for the console application.//#include "stdafx.h"#include <iostream>using namespace std;int _tmain(int argc, _TCHAR* argv[]){int m;int n;cout << "Please type the row number:" << endl;cin >> m;cout << "Please type the col number:" << endl;cin >> n;cout << "Please type the Matrix[" << m << "]" << "[" << n << "] elemnts" << endl;double **Matrix = new double *[m];for(int i = 0; i < m; i++){Matrix[i] = new double[n];}for(int i = 0; i < m; i++){for(int j = 0; j < n; j++){cin >> Matrix[i][j];}}cout << "The Matrix is:" << endl;for(int i = 0; i < m; i++){for(int j = 0; j < n; j++){cout << Matrix[i][j] << " ";}cout << endl;}cout << "The Transform Matrix is:" << endl;for(int j = 0; j < n; j++){for(int i = 0; i < m; i++){cout << Matrix[i][j] << " ";}cout << endl;}for(int i = 0; i < m; i++){delete []Matrix[i];}delete []Matrix;return 0;}

19.

// Win32Project2.cpp : Defines the entry point for the console application.//#include "stdafx.h"#include <iostream>#include <math.h>using namespace std;double det4(double M[4][4]);double det3(double M[3][3]);double (*mirrorA(double M[4][4], int r, int c,double (&Matrix)[3][3]))[3];double (*transpose(double M[4][4], double(&transposeMatrix)[4][4]))[4];double (*cofactorMatrix(double M[4][4], double (&cMatrix)[4][4]))[4];int _tmain(int argc, _TCHAR* argv[]){double Matrix[4][4];cout << "Please type the Matrix element:" << endl;for(int i = 0; i < 4; i++)for(int j = 0; j < 4; j++){cin >> Matrix[i][j];}cout << "The Matrix is:" << endl;for(int i = 0; i < 4; i++){for(int j = 0; j < 4; j++){cout << Matrix[i][j] <<" ";}cout << endl;}cout << "The Matrix det is:" << endl;double det;det = det4(Matrix);cout << det << endl;double coMatrix[4][4];cofactorMatrix(Matrix, coMatrix);double adjointMtric[4][4];transpose(coMatrix, adjointMtric);cout << "The Inverse Matrix det is:" << endl;double inverseMatrix[4][4];for(int i = 0; i < 4; i++)for(int j = 0; j < 4; j++){inverseMatrix[i][j] = adjointMtric[i][j] / det;}for(int i = 0; i < 4; i++){for(int j = 0; j < 4; j++){cout << inverseMatrix[i][j] << " ";}cout << endl;}return 0;}double det4(double M[4][4]){double det =  M[0][0]*(M[1][1]*(M[2][2]*M[3][3] -M[2][3]*M[3][2]) -M[1][2]*(M[2][1]*M[3][3] -M[2][3]*M[3][1]) +M[1][3]*(M[2][1]*M[3][2] -M[2][2]*M[3][1]))         -M[0][1]*(M[1][0]*(M[2][2]*M[3][3] -M[2][3]*M[3][2]) -M[1][2]*(M[2][0]*M[3][3] -M[2][3]*M[3][0]) +M[1][3]*(M[2][0]*M[3][2] -M[2][2]*M[3][0])) +M[0][2]*(M[1][0]*(M[2][1]*M[3][3] -M[2][3]*M[3][1]) -M[1][1]*(M[2][0]*M[3][3] -M[2][3]*M[3][0]) +M[1][3]*(M[2][0]*M[3][1] -M[2][1]*M[3][0])) -M[0][3]*(M[1][0]*(M[2][1]*M[3][2] -M[2][2]*M[3][1]) -M[1][1]*(M[2][0]*M[3][2] -M[2][2]*M[3][0]) +M[1][2]*(M[2][0]*M[3][1] -M[2][1]*M[3][0]));return det;}double det3(double M[3][3]){double det = M[0][0]*(M[1][1]*M[2][2] -M[1][2]*M[2][1]) -M[0][1]*(M[1][0]*M[2][2] -M[1][2]*M[2][0]) +M[0][2]*(M[1][0]*M[2][1] -M[1][1]*M[2][0]);return det;}double (*mirrorA(double M[4][4], int r, int c,double (&Matrix)[3][3]))[3]{int m,n;//double Matrix[3][3];for(int i = 0; i < 4; i++){if(i == r){continue;}else if(i < r){m = i;}else if(i > r){m = i-1;}for(int j = 0; j < 4; j++){if(j == c){continue;}else if(j < c){n = j;}else if(j > c){n = j-1;}Matrix[m][n] = M[i][j];}}return Matrix;}double (*transpose(double M[4][4], double(&transposeMatrix)[4][4]))[4]{//double transposeMatrix[4][4];for(int i = 0; i < 4; i++)for(int j = 0; j < 4; j++){transposeMatrix[i][j] = M[j][i];}return transposeMatrix;}double (*cofactorMatrix(double M[4][4], double (&cMatrix)[4][4]))[4]{//double cMatrix[4][4];double mMatrix[3][3];for(int i = 0; i < 4; i++)for(int j = 0; j < 4; j++){cMatrix[i][j] = pow((-1),(i+j))*det3(mirrorA(M, i, j, mMatrix));}return cMatrix;}