www.gusucode.com > optimdemos 工具箱 matlab 源码程序 > optimdemos/airPollutionProblemData.m
function [xpos, ypos, nStacks, V, diam, Ts, Q, Te, theta, U] = airPollutionProblemData %AIRPOLLUTIONPROBLEMDATA Problem data for air pollution demo % % [XPOS, YPOS, NSTACKS, V, DIAM, TS, Q, TE, THETA, U] = AIRPOLLUTIONPROBLEMDATA % returns the problem specific data for the air pollution demo. A % description of each quantity returned is given below: % % XPOS, YPOS: Positions of chimney stacks % NSTACKS: Number of chimney stacks % V: Internal stack diameters % DIAM: Pollution source diameter % TS: Temperature of gas in source % Q: Uniform emission rate % TE: Environment temperature % THETA: Wind direction % U: Mean wind speed % % See also CONCSULFURDIOXIDE, AIRPOLLUTION % Copyright 2008 The MathWorks, Inc. % Stack positions xpos = [-3000;-2600;-1100;1000;1000;2700;3000;-2000;0;1500]; ypos = [-2500;-300;-1700;-2500;2200;1000;-1600;2500;0;-1600]; %xpos = xpos([2,3,5,9,10]); %ypos = ypos([2,3,5,9,10]); nStacks = length(xpos); % Internal stack diameters V = [19.245;19.245;17.690;17.690;23.404;23.404;27.128;27.128;22.293;22.293]; %V = V([2,3,5,9,10]); % Pollution source diameter diam = [8;8;7.6;7.6;6.3;6.3;4.3;4.3;5;5]; %diam = diam([2,3,5,9,10]); % Temperature of gas in source (Kelvin) Ts = 413*ones(nStacks, 1); % Uniform emission rate Q = [2882.6;2882.6;2391.3;2391.3;2173.9;2173.9;1173.9;1173.9;1304.3;1304.3]; %Q = Q([2,3,5,9,10]); % Environmental factors % Environment temperature Te = 283; % Wind direction in rad theta = 3.996; % Mean wind speed U = 5.64;