三维坐标下的平面点阵（射线点源）

最新推荐文章于 2024-09-29 00:57:39 发布

原创最新推荐文章于 2024-09-29 00:57:39 发布 · 796 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#vector #include #delete #file #ios

计算电磁学专栏收录该内容

41 篇文章

订阅专栏

本文介绍了一个使用C++实现的三维模型光照模拟程序。该程序通过定义光线的方向、模型的包围球等参数，生成一系列用于模拟光照效果的数据点，并将这些数据写入文件，可用于进一步的光学实验或计算机图形学研究。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >


#include <math.h>
#include <stdlib.h>
#include <iostream>
#include <fstream>
using namespace std;

#include "Vector.h"

#define  PI 3.14159265
#define  CC 3E8

void main()
{
    Vector3d dir(-1,-1,-1);//光线照射方向
    
    dir = dir.normalize();//归一化入射方向

    Vector3d cen(0,0,0);//模型的包围球球心
    double R = 2;//半径
    
    Vector3d ori = cen - 2*R*dir;//入射光线原点在物体之外

    Vector3d x(dir.y, -dir.x, 0);

    if (x==Vector3d(0,0,0))
    {
        x = Vector3d(1,0,0);
    }

    x = x.normalize();

    Vector3d y = Cross(dir,x);
    y = y.normalize();
    
    double freq = 30e9;
    double lumda = CC / freq;

    double dx = 0.2;//lumda / 10 ;
    double dy = dx;

    int NX= (int) (2*R/dx);
    int NY= (int) (2*R/dy);
    
    Vector3d** points = new Vector3d*[NX];
    
    for(int i=0;i<NX;i++)
    {
        points[i] = new Vector3d[NY];
    }
        
    fstream file("outdata2.txt",ios::out);

    file<<"x"<<"\t"<<"y"<<"\t"<<"z"<<endl;

    for(int nx = 0;nx < NX;nx++)
    {
        for(int ny=0;ny<NY;ny++)
        {
            points[nx][ny] = ori+dx*(nx-NX/2)*x+dy*(ny-NY/2)*y;
        
            file<<points[nx][ny].x<<"\t"<<points[nx][ny].y<<"\t"<<points[nx][ny].z<<endl;
        }
    }
    
    for(int i=0;i<NX;i++)
    {
        delete []points[i];
    }
    
    file.close();

    cout<<"done!"<<endl;

    getchar();
}