我是巴斯肯特大学(土耳其，安卡拉(计算机工程专业的学生。
我可以在visual studio中通过导入或创建库之类的方式使用matlab k波工具箱代码吗？我需要知道我的Graduation项目。
例如：
% 2D Time Reversal Reconstruction For A Line Sensor Example
%
% This example demonstrates the use of k-Wave for the time-reversal
% reconstruction of a two-dimensional photoacoustic wave-field recorded
% over a linear array of sensor elements. The sensor data is simulated and
% then time-reversed using kspaceFirstOrder2D. It builds on the 2D FFT
% Reconstruction For A Line Sensor Example.
%
% author: Bradley Treeby
% date: 6th July 2009
% last update: 25th July 2019
%
% This function is part of the k-Wave Toolbox (http://www.k-wave.org)
% Copyright (C) 2009-2019 Bradley Treeby
% This file is part of k-Wave. k-Wave is free software: you can
% redistribute it and/or modify it under the terms of the GNU Lesser
% General Public License as published by the Free Software Foundation,
% either version 3 of the License, or (at your option) any later version.
%
% k-Wave is distributed in the hope that it will be useful, but WITHOUT ANY
% WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
% FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for
% more details.
%
% You should have received a copy of the GNU Lesser General Public License
% along with k-Wave. If not, see <http://www.gnu.org/licenses/>.
clearvars;
% =========================================================================
% SIMULATION
% =========================================================================
% create the computational grid
PML_size = 20; % size of the PML in grid points
Nx = 128 - 2 * PML_size; % number of grid points in the x direction
Ny = 256 - 2 * PML_size; % number of grid points in the y direction
dx = 0.1e-3; % grid point spacing in the x direction [m]
dy = 0.1e-3; % grid point spacing in the y direction [m]
kgrid = kWaveGrid(Nx, dx, Ny, dy);
% define the properties of the propagation medium
medium.sound_speed = 1500; % [m/s]
% create initial pressure distribution using makeDisc
disc_magnitude = 5; % [Pa]
disc_x_pos = 60; % [grid points]
disc_y_pos = 140; % [grid points]
disc_radius = 5; % [grid points]
disc_2 = disc_magnitude * makeDisc(Nx, Ny, disc_x_pos, disc_y_pos, disc_radius);
disc_x_pos = 30; % [grid points]
disc_y_pos = 110; % [grid points]
disc_radius = 8; % [grid points]
disc_1 = disc_magnitude * makeDisc(Nx, Ny, disc_x_pos, disc_y_pos, disc_radius);
% smooth the initial pressure distribution and restore the magnitude
p0 = smooth(disc_1 + disc_2, true);
% assign to the source structure
source.p0 = p0;
% define a binary line sensor
sensor.mask = zeros(Nx, Ny);
sensor.mask(1, :) = 1;
% create the time array
kgrid.makeTime(medium.sound_speed);
% set the input arguements: force the PML to be outside the computational
% grid; switch off p0 smoothing within kspaceFirstOrder2D
input_args = {'PMLInside', false, 'PMLSize', PML_size, 'Smooth', false, 'PlotPML', false};
% run the simulation
sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor, input_args{:});
% reset the initial pressure
source.p0 = 0;
% assign the time reversal data
sensor.time_reversal_boundary_data = sensor_data;
% run the time reversal reconstruction
p0_recon = kspaceFirstOrder2D(kgrid, medium, source, sensor, input_args{:});
% add first order compensation for only recording over a half plane
p0_recon = 2 * p0_recon;
% repeat the FFT reconstruction for comparison
p_xy = kspaceLineRecon(sensor_data.', dy, kgrid.dt, medium.sound_speed, ...
'PosCond', true, 'Interp', '*linear');
% define a second k-space grid using the dimensions of p_xy
[Nx_recon, Ny_recon] = size(p_xy);
kgrid_recon = kWaveGrid(Nx_recon, kgrid.dt * medium.sound_speed, Ny_recon, dy);
% resample p_xy to be the same size as source.p0
p_xy_rs = interp2(kgrid_recon.y, kgrid_recon.x - min(kgrid_recon.x(:)), p_xy, kgrid.y, kgrid.x - min(kgrid.x(:)));
% =========================================================================
% VISUALISATION
% =========================================================================
% plot the initial pressure and sensor distribution
figure;
imagesc(kgrid.y_vec * 1e3, kgrid.x_vec * 1e3, p0 + sensor.mask * disc_magnitude, [-disc_magnitude, disc_magnitude]);
colormap(getColorMap);
ylabel('x-position [mm]');
xlabel('y-position [mm]');
axis image;
colorbar;
scaleFig(1, 0.65);
% plot the reconstructed initial pressure
figure;
imagesc(kgrid.y_vec * 1e3, kgrid.x_vec * 1e3, p0_recon, [-disc_magnitude, disc_magnitude]);
colormap(getColorMap);
ylabel('x-position [mm]');
xlabel('y-position [mm]');
axis image;
colorbar;
scaleFig(1, 0.65);
% apply a positivity condition
p0_recon(p0_recon < 0) = 0;
% plot the reconstructed initial pressure with positivity condition
figure;
imagesc(kgrid.y_vec * 1e3, kgrid.x_vec * 1e3, p0_recon, [-disc_magnitude, disc_magnitude]);
colormap(getColorMap);
ylabel('x-position [mm]');
xlabel('y-position [mm]');
axis image;
colorbar;
scaleFig(1, 0.65);
% plot a profile for comparison
figure;
plot(kgrid.y_vec * 1e3, p0(disc_x_pos, :), 'k-', ...
kgrid.y_vec * 1e3, p_xy_rs(disc_x_pos, :), 'r--', ...
kgrid.y_vec * 1e3, p0_recon(disc_x_pos, :), 'b:');
xlabel('y-position [mm]');
ylabel('Pressure');
legend('Initial Pressure', 'FFT Reconstruction', 'Time Reversal');
axis tight;
set(gca, 'YLim', [0, 5.1]);
该代码基本上创建了一个线传感器的2D时间反转重建图像
在这个例子中，我想用变量为这个代码创建一个函数
returnValue twoDreversalLineSensor(int scaler, int soundSpeed, int xPos, int yPos)-因为这是一个输入是随机的例子-
我能做这个吗？
如果我能，我该怎么办？
在视觉工作室之后，我可以接受用户的输入吗-因为这个代码来自其他地方-
谢谢你所做的一切。
顺便说一句，对不起我的英语，我写得不好。
这不是一条无故障的路径，但您可以使用matlab引擎，请参阅的示例
https://www.mathworks.com/help/matlab/matlab_external/calling-matlab-software-from-a-c-application.html
基本上，您可以调用engEvalString()在后端的一个不可见的matlab会话中运行matlab命令。
如果您只需要一个结果，您可以使用系统调用(ShellExecute或ShellExecuteEx(并调用
/path/to/matlab -nojvm -nodesktop < /path/to/yourscript.m > cmdoutput.txt
调用matlab会话。