![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?5a7ff87)
学校の課題で,二つのファイルから得たデータで相互相関関数を計算する,という課題がでました.
以下のようなプログラムを書いたのですが,コンパイルの時点でエラーはなく,実行したときに
malloc(): corrupted top size
Aborted (core dumped)
というエラーがでてきます.
この原因がもしわかれば教えていただきたいです.
また,他にもおかしな点があれば教えてください.
よろしくお願いいたします.
以下,プログラム
#include <stdio.h>
#include <stdlib.h>
#include <stddef.h>
#define BUF_SIZE 4000
int main()
{
/* open file */
FILE *fp_x = fopen("ecg400.txt", "r");
if (fp_x == NULL){
printf("Can't open data file.\n");
return 1;
}
FILE *fp_y = fopen("ecg4000.txt", "r");
if (fp_y == NULL){
printf("Can't open data file.\n");
return 1;
}
/* read data */
char buf[BUF_SIZE];
int N_x = 0, N_y = 0, k, y;
int a = 0, a2 = 2, b;
int line = 3000000;
float *etime_x, *resp_x, *etime_y, *resp_y, *resp2_x, *resp2_y;
float m = 0, *d, *e, *f, *peak;
float r, time_ave, bpm_ave;
etime_x = (float*)malloc(sizeof(float)*N_x);
resp_x = (float*)malloc(sizeof(float)*N_x);
etime_y = (float*)malloc(sizeof(float)*N_y);
resp_y = (float*)malloc(sizeof(float)*N_y);
fgets(buf, BUF_SIZE-1, fp_x);
sscanf(buf, "%d\n", &N_x);
fgets(buf, BUF_SIZE-1, fp_y);
sscanf(buf, "%d\n", &N_y);
for (int i = 0; i < N_x ; i++ ){
fgets(buf, BUF_SIZE-1, fp_x);
sscanf(buf, "%f,%f\n", &etime_x[i], &resp_x[i]);
}
for (int i = 0; i < N_y ; i++ ){
fgets(buf, BUF_SIZE-1, fp_y);
sscanf(buf, "%f,%f\n", &etime_y[i], &resp_y[i]);
}
b = N_y - N_x;
resp2_x = (float*)malloc(sizeof(float)*N_x);
resp2_y = (float*)malloc(sizeof(float)*N_y);
peak = (float*)malloc(sizeof(float)*b);
d = (float*)malloc(sizeof(float)*(b+1));
e = (float*)malloc(sizeof(float)*(b+1));
f = (float*)malloc(sizeof(float)*b);
printf("ノイズ除去(移動平均処理)しますか?(はい→1)\n");
scanf("%d", &y);
/* ノイズ除去するとき */
if(y == 1){
for(int i = a2-1; i <= N_x-1; i++ ){
for(int j = 0; j <= a2-1; j++ ){
resp2_x[i] += resp_x[i-j];
}
resp2_x[i] /= a2;
}
for(int i = a2-1; i <= N_y-1; i++ ){
for(int j = 0; j <= a2-1; j++ ){
resp2_y[i] += resp_y[i-j];
}
resp2_y[i] /= a2;
}
for(int i = a2-1; i <= N_x-1; i++ ){
resp_x[i] = resp2_x[i];
}
for(int i = a2-1; i <= N_y-1; i++ ){
resp_y[i] = resp2_y[i];
}
a = a2 - 1;
}
for(int i = a; i <= N_y-1; i++ ){
m += resp_y[i];
}
m /= (N_y - a);
for(int i = a; i <= N_x-1; i++){
resp_x[i] -= m;
}
for(int i = a; i <= N_y-1; i++){
resp_y[i] -= m;
}
for(int i = 0 ; i <= b ; i++){
for(int j = a; j <= N_x-1; j++){
d[i] += resp_x[j] * resp_y[j+i];
}
}
for(int i = 1; i <= b ; i++){
e[i] = d[i] - d[i-1];
}
for(int i = 1; i <= (b - 1) ; i++){
f[i] = e[i]*e[i+1];
}
for(int i = 1; i <= (b - 1) ; i++){
if(f[i] < 0){
if(d[i] >= line){
peak[k] = etime_y[i];
k++;
}
}
}
for(int i = 1; i <= k-1; i++){
r = peak[i] - peak[i-1];
printf("R波間隔%d:%f\n", i, r);
}
time_ave = (peak[k-1] - peak[0])/5;
bpm_ave = 60/time_ave;
printf("R波間隔の平均値:%f\n", time_ave);
printf("心拍数:%f\n", bpm_ave);
fclose(fp_x);
fclose(fp_y);
free(etime_x);
free(resp_x);
free(resp2_x);
free(etime_y);
free(resp_y);
free(resp2_y);
free(d);
free(e);
free(f);
free(peak);
}
No.2ベストアンサー
- 回答日時:
N_xとかに値いれる前に
メモリ確保に使ってるよね。
ご返信ありがとうございます!理解できました
etime_xなどを定義する前に,fgetsでN_xなどに値いれなければってことですよね!
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- C言語・C++・C# プログラミングの授業の課題です 1 2023/01/17 22:15
- C言語・C++・C# c言語の問題です 課題1 (二分探索木とセット) 大きさ size の配列 array を考える。す 2 2023/01/10 21:08
- C言語・C++・C# leetcode 155 minstack 1 2022/05/07 16:43
- C言語・C++・C# c言語の問題の説明、各所ごとに 5 2023/07/26 11:03
- C言語・C++・C# #include <stdio.h>int main(void) { int buf[100] = 6 2022/11/01 22:45
- C言語・C++・C# 10個の実数に対する降順ソート結果を出力するプログラムを作りたいのですが、以下のプログラムをどう直せ 1 2022/07/09 22:16
- C言語・C++・C# このプログラミング誰か教えてくれませんか 1 2022/06/02 15:27
- C言語・C++・C# C++ と、 1 2022/11/07 23:45
- C言語・C++・C# c言語の問題です 3 2023/01/10 16:15
- C言語・C++・C# C言語で再起関数とポインタを用いて文字列反転をする方法がわかりません。 4 2023/04/29 20:32
このQ&Aを見た人はこんなQ&Aも見ています
-
プロが教える店舗&オフィスのセキュリティ対策術
中・小規模の店舗やオフィスのセキュリティセキュリティ対策について、プロにどう対策すべきか 何を注意すべきかを教えていただきました!
-
C 開放してるのにエラー(double free or corruption (!prev))がでる
C言語・C++・C#
-
C言語 配列の長さの上限
C言語・C++・C#
-
セグメントエラー
C言語・C++・C#
-
-
4
メモリのセグメント違反の解決方法を教えてください。
C言語・C++・C#
-
5
C++でアボート(Abort)で処理が強制終了してしまう
C言語・C++・C#
-
6
C言語初心者の質問失礼します。
C言語・C++・C#
-
7
バッファとは何ですか
C言語・C++・C#
-
8
char*を初期化したいのですが
C言語・C++・C#
-
9
break文でループを一気に抜けるには
C言語・C++・C#
-
10
std::set<int> で、ある値が何番目の数なのか知る方法
C言語・C++・C#
-
11
e^(x^2)の積分に関して
数学
-
12
エラーの意味は? Lvalue required
C言語・C++・C#
-
13
変数名の付け方
Java
-
14
C言語 exitの使い方
C言語・C++・C#
-
15
fopne で失敗する原因
C言語・C++・C#
-
16
複数桁10進数の*桁目だけを抽出したい
C言語・C++・C#
-
17
malloc.cを見たい
UNIX・Linux
-
18
free()関数の多用 と Segmentation fault
C言語・C++・C#
-
19
define で 配列
C言語・C++・C#
-
20
gccのインクルードパスの確認
UNIX・Linux
関連するカテゴリからQ&Aを探す
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
複数桁10進数の*桁目だけを抽出...
-
「指定されたキャストは有効で...
-
【C++】関数ポインタの使い方
-
#define _CRT_SECURE_NO_WARNIN...
-
C言語 エラーの原因がわからな...
-
c++でテンプレートのコードでわ...
-
C言語での引数の省略方法
-
C言語の課題です
-
アスタリスクでダイヤ型を作る
-
困ってます…nCrを求めるC言語...
-
(int *)の意味
-
int型とfloat型の演算結果
-
acceptをalarmでタイムアウトさ...
-
答えがおかしい
-
任意の文字列のアルファベット...
-
C言語、行列の積を求めるプログ...
-
C++でRPGを作成する際のステー...
-
C言語に関する質問です フィボ...
-
日数チェックプログラム
-
ColorをRGBで指定する方法
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
「指定されたキャストは有効で...
-
C言語での引数の省略方法
-
#define _CRT_SECURE_NO_WARNIN...
-
複数桁10進数の*桁目だけを抽出...
-
へんな現象
-
【C++】関数ポインタの使い方
-
C言語 エラーの原因がわからな...
-
if と配列の組み合わせ
-
C言語での奇数の和
-
C言語 配列と関数の練習問題
-
ラップ関数とはどんなものですか?
-
(int *)の意味
-
C言語
-
実数の整数部,小数部の取得
-
足して100になるような乱数のア...
-
卒業研究でよく分からないとこ...
-
数字列を3桁ごとにカンマで区切...
-
c言語
-
std::set<int> で、ある値が何...
-
比較回数と交換回数表示について
おすすめ情報