版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

eamtbfs

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 528.6
帖子: 209
在线: 35.6小时
虫号: 533831

[交流] 【求助】帮忙看看哪错了error: invalid lvalue in assignment

我最近用到FFTW，就网上找了个简单例子来学习了下。但这个简单例子编译出错，不知道怎么回事。我网上查了查，也没搞清是FFTW语句的问题，还是gcc编译器的问题。请高手帮忙看看啊！

错误信息如下：
$gcc34 FFTW_example.c -lfftw3 -lm
FFTW_example.c: In function `main':
FFTW_example.c:35: error: invalid lvalue in assignment
FFTW_example.c:36: error: invalid lvalue in assignment

用gcc4.3也是一样出错。

这个FFTW的简单例子程序的网址在这里http://www.elisanet.fi/mnentwig/webroot/FFTW_example/index.html

程序比较小就贴在下面。是creal 和cimag那两行出错了
/* ****************************************************
* FFTW example
* Markus Nentwig, 2007
* This program is in the public domain
* ****************************************************/
#include
#include
#include
#include
#include

int main(void){
  /* ****************************************************
* Allocate memory
* ****************************************************/
  int n=128;
  fftw_complex* b1=fftw_malloc(sizeof(fftw_complex)*n);
  fftw_complex* b2=fftw_malloc(sizeof(fftw_complex)*n);

  /* ****************************************************
* Create forward FFT plan from b1 into b2
* ****************************************************/
  fftw_plan p1=fftw_plan_dft_1d(n, b1, b2, FFTW_FORWARD, FFTW_ESTIMATE);

  /* ****************************************************
* Create reverse FFT plan from b2 into b1
* ****************************************************/
  fftw_plan p2=fftw_plan_dft_1d(n, b2, b1, FFTW_BACKWARD, FFTW_ESTIMATE);

  /* ****************************************************
* Generate a single sine wave cycle
* ****************************************************/
  int i;
  for (i=0; i < n; ++i){
creal(b1)=sin((double)i/(double)n*2.0*M_PI);
cimag(b1)=0;
  }

  /* ****************************************************
* Execute FFT from b1 into b2
* ****************************************************/
  fftw_execute(p1);

  /* ****************************************************
* Dump the "spectrum".
* Result:
* First bin ("DC  term" is zero
* * i*n/2 in 2nd bin
* * -i*n/2 in last bin
* * all other bins are zero
* ****************************************************/
  for (i=0; i < n; ++i){
printf("%1.7f\t%1.7f\n", creal(b2), cimag(b2));
  };

  /* ****************************************************
* IFFT back into b1
* ****************************************************/
  fftw_execute(p2);
  printf("\n";

  /* ****************************************************
* Result: Original waveform scaled by n
* ****************************************************/
  for (i=0; i < n; ++i){
printf("%1.7f\t%1.7f\n", creal(b1), cimag(b1));
  };

  /* ****************************************************
* Clean up
* ****************************************************/
  fftw_destroy_plan(p1);
  fftw_destroy_plan(p2);

  fftw_free(b1);
  fftw_free(b2);
  return 0;
}

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】中科大分析化学的范围已经有2人回复
【求助】液相色谱，老是出现泵里面有空气的情况？已经有3人回复
【求助】液相色谱操作相关问题，请教！！已经有17人回复
【求助】按论坛方法安装Gaussian09出错已经有10人回复
【求助】各位前辈看过来，请教关于Dmol3的问题！急~~~~~ 已经有0人回复
【求助/交流】cd31抗体有哪位虫友用过，麻烦给推荐个质量可靠的兔抗已经有1人回复
【求助】请问什么是淋巴癌的前期征兆已经有2人回复
【整理】《提问的智慧》——献给那些不会提问的人已经有78人回复
【求助】帮忙看看这是什么报错已经有3人回复
【求助】一个fortran程序转换为VC++ 程序已经有2人回复
【求助】攻博该去哪？已经有28人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-02-24 14:44:59

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yalefield

金虫 (文坛精英)

★ ★ ★
eamtbfs(金币+1): 谢谢参与 2011-02-25 06:01:14
余泽成(金币+3): 谢谢参与应助！ 2011-02-27 15:06:29

抄都没抄对咯

CODE:

    creal(b1[i])=sin((double)i/(double)n*2.0*M_PI);

    cimag(b1[i])=0;

CODE:

printf("%1.7f\t%1.7f\n", creal(b1[i]), cimag(b1[i]));

赞一下(1人)

回复此楼

2楼2011-02-24 19:31:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

eamtbfs

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 528.6
帖子: 209
在线: 35.6小时
虫号: 533831

zzzzz

[ Last edited by eamtbfs on 2011-2-25 at 06:06 ]

回复此楼

3楼2011-02-25 06:04:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

eamtbfs

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 528.6
帖子: 209
在线: 35.6小时
虫号: 533831

引用回帖:

Originally posted by yalefield at 2011-02-24 19:31:36:
抄都没抄对咯

CODE:

    creal(b1[i])=sin((double)i/(double)n*2.0*M_PI);

    cimag(b1[i])=0;

CODE:

printf("%1.7f\t%1.7f\n", creal(b1[i]), cimag(b1[i]));

我从网页上帖过来的，不知怎么贴过来就变了。不过我编译的时候是用的网上下载的程序，我仔细看了，我编译的程序语法没有这样的错，但就是gcc报告error: invalid lvalue in assignment

请帮忙编译下，看看问题出在哪。多谢多谢！

赞一下

回复此楼

4楼2011-02-25 06:05:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lucktroy

禁虫 (小有名气)

★ ★ ★
余泽成(金币+3): 谢谢参与应助！ 2011-02-28 20:09:22

本帖内容被屏蔽

5楼2011-02-28 19:25:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

holmescn

金虫 (正式写手)

程序强帖: 37
应助: 1 (幼儿园)
金币: 1918.8
帖子: 699
在线: 102.6小时
虫号: 913482

★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-03-07 21:38:09

没用过FFTW，不过这两行确实比较奇怪

CODE:

  /* ****************************************************

   * Generate a single sine wave cycle

   * ****************************************************/

  int i;

  for (i=0; i < n; ++i){

    creal(b1[i])=sin((double)i/(double)n*2.0*M_PI);

    cimag(b1[i])=0;

  }

creal和cimag是两个函数还是两个宏啊。如果是两个函数，那这里肯定是有问题的。如果是两个宏，不知道会展开成什么。

赞一下(1人)

回复此楼

6楼2011-03-04 09:46:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖