試著改編乙個fftw的傅利葉程式,但是老是出問題,希望有懂的的大神幫忙看哈 指出怎麼改!
#include "fftw3.h"
#include
#include
#include
#include
#include
//const int n=4240 ;
//const double pi = 3.1416;
//const double fn =4000;
void fftw (double fn,double max_a[2],int n)
fprintf (output, "%9.6f %e\n", i*fn/n, sqrt(out[i][0]*out[i][0] + out[i][1]*out[i][1]));
countout++;
}printf("output is %d\n", countout);
fprintf (output, "******************************=== \n");
if ( fclose (input))
if ( fclose (output))
fftw_destroy_plan(p);
fftw_free(in);
fftw_free(out);
}
void trans(double fn,double max_a, long pick,long n)
if (!(output = fopen (outputfile, "a")))
if (!(input2 = fopen(input2file, "a")))
for ( i=0;i
int count = 0;
for(i=0;i< n;i++)
fprintf (input2, "******************************=== \n");
printf("input is %d \n", count);
fftw (fn,max_a,n);
}void fft_test (double fn,long p_pick, long s_pick)
while ((fscanf (file1, "%lf %lf", &tem1, &tem1)) != eof)//讀取試驗資料
num_end=count;
s_long=num_end-s_pick+1;
trans(fn,s_pick,s_long,max_as);//s波傅利葉變換
s波傅利葉變換
file *output;
if (!(output = fopen ("max_amplitude_fre.dat", "w")))
fprintf (output, "frequency amplittude \n");
if(max_as[1]>max_ap[1])
else
if ( fclose (output))
}int main ()
問題如下:
1>------ 已啟動全部重新生成: 專案: fftw_p_s0, 配置: debug win32 ------
1>生成啟動時間為 2016-8-6 20:52:09。
1>_prepareforclean:
1> 正在刪除檔案「debug\fftw_p_s0.lastbuildstate」。
1>initializebuildstatus:
1> 正在對「debug\fftw_p_s0.unsuccessfulbuild」執行 touch 任務。
1>clcompile:
1> fftw.cpp
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(14): warning c4101: 「input2」: 未引用的區域性變數
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(77): warning c4996: 'fopen': this function or variable may be unsafe. consider using fopen_s instead. to disable deprecation, use _crt_secure_no_warnings. see online help for details.
1> d:\vs2010\vc\include\stdio.h(234) : 參見「fopen」的宣告
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(82): warning c4996: 'fopen': this function or variable may be unsafe. consider using fopen_s instead. to disable deprecation, use _crt_secure_no_warnings. see online help for details.
1> d:\vs2010\vc\include\stdio.h(234) : 參見「fopen」的宣告
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(87): warning c4996: 'fopen': this function or variable may be unsafe. consider using fopen_s instead. to disable deprecation, use _crt_secure_no_warnings. see online help for details.
1> d:\vs2010\vc\include\stdio.h(234) : 參見「fopen」的宣告
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(94): warning c4996: 'fscanf': this function or variable may be unsafe. consider using fscanf_s instead. to disable deprecation, use _crt_secure_no_warnings. see online help for details.
1> d:\vs2010\vc\include\stdio.h(253) : 參見「fscanf」的宣告
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(101): warning c4996: 'fscanf': this function or variable may be unsafe. consider using fscanf_s instead. to disable deprecation, use _crt_secure_no_warnings. see online help for details.
1> d:\vs2010\vc\include\stdio.h(253) : 參見「fscanf」的宣告
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(119): error c2664: 「trans」: 不能將引數 2 從「long」轉換為「double 」
1> 從整型轉換為指標型別要求 reinterpret_cast、c 樣式轉換或函式樣式轉換
1>c:\documents and settings\dell\桌面\fftw_p_s\fftw_p_s0\fftw.cpp(140): error c2664: 「trans」: 不能將引數 2 從「long」轉換為「double 」
1> 從整型轉換為指標型別要求 reinterpret_cast、c 樣式轉換或函式樣式轉換
1>
1>生成失敗。
1>
1>已用時間 00:00:00.29
********** 全部重新生成: 成功 0 個,失敗 1 個,跳過 0 個 **********
小菜鳥成長記
include define n 100 typedef struct book book int num 0,pnum int arrin book a,int pnum 輸入陣列元素,個數num由程式中自己確定。return i void arrout book a,int num 輸出陣列各元...
python小菜鳥入門
import os os.getcwd os.chdir os.getcwd name input micheal 接下來要執行的 部分需要首行縮排4位 官方標準,3位也能執行 當字串中包含 或者 的時候需要使用轉義字元,for instance i m ok 輸出就是 i m ok 其中 n表示換...
剛入行的小菜鳥
慢慢的,老大會讓做些客戶的小需求,參照別人的 進行一些縫縫補補,有時候做出來,還是蠻有成就感的額。但是這一兩個月來,我並沒有體會到,我到底能做什麼,前端 資料庫 後端,感覺啥都不懂,隨便難點的問題,就可以讓我卡上一天。想學東西,又不知道從何學,真是新人連下手的放向都沒有。前些日子,表哥讓我沒事多逛論...