ピクニック 091013
↓チェックせよ!!
●下記のプログラムで、s==0の所なのですがピポット操作の後に行いたいのですがどの....
下記のプログラムで、s==0の所なのですがピポット操作の後に行いたいのですがどのようなプログラムにすればよいのでしょうか?また、s<=0.00001にして、逆行列が存在しないと判断させたいと思っています。#include<stdio.h> #define M 5 int Input(double a[M][M]) { int i,j; int n; printf("行列の次数N入力(N=1~5)\\\\n",M); fflush(stdin); scanf("%d",&n); fflush(stdin); printf("%d行%d列\\\\n",n,n); fflush(stdin); for(i=0;i<n;i++){ for(j=0;j<n;j++){ printf("%d行%d列の値を入力\\\\n",i+1,j+1); scanf("%f",&a[i][j]); fflush(stdin); } } return n; } int Inverse(int n,double a[M][M],double b[M][M]) { int i,j,k; double s,t; for(i=0;i<n;i++){ for(j=0;j<n;j++){ if(i!=j) b[i][j]=0; else b[i][j]=1; } } for(k=0;k<n;k++){ s=a[k][k]; if(s==0){ printf("逆行列は存在しません。\\\\n"); return -1; } for(j=0;j<n;j++){ a[k][j]=a[k][j]/s; b[k][j]=b[k][j]/s; } for(i=0;i<n;i++){ if(i!=k){ t=a[i][k]; for(j=0;j<n;j++){ a[i][j]=a[i][j]-t*a[k][j]; b[i][j]=b[i][j]-t*b[k][j]; } } } } return 0; } void Output(int n,double b[M][M]) { int i,j; printf("入力した逆行列は、\\\\n"); for(i=0;i<n;i++){ for(j=0;j<n;j++){ printf("%d行%d列 = %f \\\\n",i+1,j+1,b[i][j]); fflush(stdin); } } } int main() { int err,n; double a[M][M]; double b[M][M]; n=Input(a); err=Inverse(n,a,b); if(!err)Output(n,b); return 0; }
続き
---
●このC++で書かれたプログラムをC言語に直してくれると嬉しいです。
このC++で書かれたプログラムをC言語に直してくれると嬉しいです。#include "iostream" #define _USE_MATH_DEFINES #include "cmath" using namespace std; #define PI M_PI #define n 16 #define m 4 int main (void){ int i=0,j=0,k=0,l=0,o=0,p=0,q=0,r=0,s=0; double xr[n]; double xi[n]; for(i=0;i<n;i++){ xr[i]=sin((double)2.0*1*PI*i/n)+cos((double)2.0*1*PI*i/n); xi[i]=0.0; } cout << "input Data:" << endl; for(i=0;i<n;i++){ cout << "xr[" << i << "]:" << xr[i] << ", xi[" << i << "]: " << xi[i] << endl; } double a=0.0; double b=PI*2/n; double xd=0.0; double S[n/2+1]; double C[n/2+1]; for(i=0;i<=n/2;i++){ S[i]=sin(a); C[i]=cos(a); a=a+b; } s=1; l=n; for(r=1;r<=m;r++){ l=l/2; k=0; for(q=1;q<=s;q++){ p=0; for(i=k;i<=l+k-1;i++){ j=i+l; a=xr[i]-xr[j]; b=xi[i]-xi[j]; xr[i]=xr[i]+xr[j]; xi[i]=xi[i]+xi[j]; if(p==0){ xr[j]=a; xi[j]=b; } else{ xr[j]=a*C[p]+b*S[p]; xi[j]=b*C[p]-a*S[p]; } p=p+s; } k=k+l+l; } s=s+s; } j=n/2; for(i=1;i<n-1;i++){ k=n; if(j<i){ xd=xr[i]; xr[i]=xr[j]; xr[j]=xd; xd=xi[i]; xi[i]=xi[j]; xi[j]=xd; } k=k/2; while(j>=k){ j=j-k; k=k/2; } j=j+k; } cout << "Output Data:" << endl; for(i=0;i<n;i++){ cout << "xr[" << i << "]:" << xr[i] << ", xi[" << i << "]: " << xi[i] << endl; } return 0; }
続き
---
●プログラミング C言語において、 行列の 行数 m 整数型 列数 n 整数型 要素 a...
プログラミング C言語において、 行列の 行数 m 整数型 列数 n 整数型 要素 a[M][N] double型の2次元配列 をまとめて扱うmatrixという構造体を定義し、行列の積の計算をするプログラムを作ってみる。 分かる方、教えてくださいm(_ _)m
続き
---
●随筆『Where is PICNIC的なモノ?』 : STRIKERS 「黄金の窓口」
この件に関しては、You tubeにも動画をあげたし、もうこのblogでも何度も書きましたし、今後も書く気がします。 で、そんなこんなで自動車という手段を失った自分は、もう完全に、BUS&TRAIN、そう、ブスと電車で移動。。。。違う違う。バスと電車を使っての 移動になってしまい、 ... 森的な場所も歩いたし、人の優しさにも触れたし、なんだかピクニックみたいで良かったじゃん!と自分に言い聞かせようとしましたが、坂下ったバス停にて主に通り過ぎるトラックの排気ガスがピクニック気分を吹き飛ばすのでした。 ...
続き
---
●随筆『Where is PICNIC的なモノ?』
... 件に関しては、You tubeにも動画をあげたし、もうこのblogでも何度も書きましたし、今後も書く気がします。 ... 場所も歩いたし、人の優しさにも触れたし、なんだかピクニックみたいで良かったじゃん!と自分に言い聞かせようとし ... 【友情出演】響(お笑い芸人)
続き
---
↓気になるサイト
これはいい!【パチスロデータシート】ver2.0&2.1(エクセル) 丸山直樹
これは使える!【競馬道GT8】競馬予想支援ソフト
初心者でも勝てるギャンブル商材徹底検証するブログ
クロニクル2010 2009年10月〜12月だけでも53万円もの利益をたたき出した、大好評!競馬で勝ち続ける男・大島和也の競馬予想配信サービスの結果検証報告ブログ
2010年春のG1を的中させながら競馬を楽しむための★情報★の結果 検証 報告プログ
2010年06月17日
億馬〜3連複2点的中法 的中率脅威50%超!
RSS取得

