| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 271 人关注过本帖
标题:高手帮忙看看这个程序错在哪里?万分感谢!
取消只看楼主 加入收藏
光学衍射
Rank: 1
等 级:新手上路
帖 子:5
专家分:0
注 册:2010-11-15
收藏
 问题点数:0 回复次数:0 
高手帮忙看看这个程序错在哪里?万分感谢!
可以运行,但是输出的值是一样的..请大神帮忙改下!万分感谢
#include<stdio.h>
#include<math.h>
#define PI 3.1415926

double d, lad, dltL,A0,dlty,dltx,R,sita;
int N;
struct complex
{
    double read,imag;
};
complex fadd(complex fs1,complex fs2)   
{
    complex fs3;
    fs3.read=fs1.read+fs2.read;
    fs3.imag=fs1.imag+fs2.imag;
    return fs3;
}

complex Aup(double pointx,double pointy)
{double pointxp,pointyp;
    complex rst;
    rst.read=0;
    rst.imag=0;
    lad=400;  
    lad*=1e-6;
    double k=2*PI/lad;           
    for(int ii=-49;ii<49;ii++)
    for(int jj=-49;jj<49;jj++)
    {
         pointx=ii*R/50;
         pointy=jj*R/50;         
        double r=sqrt((pointyp-pointy)*(pointyp-pointy)+(pointxp-pointx)*(pointxp-pointx)+d*d);
        if(sqrt((pointx*pointx)+(pointy*pointy))>=R)continue;
        complex lcpx;
        lcpx.read=cos(k*r)/r;
        lcpx.imag=sin(k*r)/r;
        rst=fadd(rst,lcpx);
        
    }
    //rst=rst*A0;
    rst.read=rst.read*A0;
    rst.imag=rst.imag*A0;
    return rst;
}

int main()
{
    double pointyp,pointxp;
       complex AP,IP;
     printf("input a0:\n");
    scanf("%lf",&A0);
    printf("input dltx:\n");
    scanf("%lf",&dltx);
    printf("input d:\n");
    scanf("%lf",&d);
    printf("input R:\n");
    scanf("%lf",&R);
    printf("input dlty:\n");
    scanf("%lf",&dlty);
    printf("input sita(弧度制):\n");
       scanf("%f",&sita);
    for(int mm=1;mm<=200;mm++)
    for(int ss=1;ss<=200;ss++)
    {
        pointxp=10+dltx*mm;
        pointyp=10+dlty*ss;
        AP=Aup(pointxp,pointyp);
        IP.read=AP.read*AP.read;
        IP.imag=AP.imag*AP.imag;
        
        printf("IP=%.5f+i%.5f\n",IP.read,IP.imag);
    }
    return 0;
}
搜索更多相关主题的帖子: 感谢 
2010-11-24 10:11
快速回复:高手帮忙看看这个程序错在哪里?万分感谢!
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.016568 second(s), 10 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved