如何用MATLAB的fzero求解给定方程零点？解决变量未定义报错

That "Undefined function or variable" error you're seeing almost always comes down to not properly passing your fixed parameters (C, Tₘ, Tᵢ, Lբ) to the function you're trying to solve with fzero. Let's break down exactly how to fix this, with two straightforward approaches:

Approach 1: Use an Anonymous Function (Quickest Method) #

This is the most common way to handle fixed parameters in MATLAB for solver functions like fzero. We'll wrap your equation into an anonymous function that binds your input parameters, so fzero only needs to handle the variable lambda.

Here's the full code:

% Step 1: Get user input for your parameters
C = input('Enter value of C: ');
Tm = input('Enter value of Tm: ');
Ti = input('Enter value of Ti: ');
Lf = input('Enter value of Lf: ');

% Step 2: Precompute the constant term to simplify the equation
constant_term = C*(Tm - Ti)/(Lf*sqrt(pi));

% Step 3: Define the anonymous function with lambda as the input variable
f = @(lambda) lambda * exp(lambda^2) * erfc(lambda) - constant_term;

% Step 4: Call fzero with an initial guess (choose a value close to the expected root)
% For your equation, lambda > 0 is typical, so start with x0=1 as a guess
lambda_solution = fzero(f, 1);

% Step 5: Display the result
disp(['Solution for lambda: ', num2str(lambda_solution)]);

Why this works: #

The anonymous function @(lambda) tells MATLAB that lambda is the variable to solve for, while all other parameters (C, Tm, etc.) are already defined and bound to the function. fzero can now evaluate the function for different lambda values without missing any variables.

Approach 2: Create a Separate Function File #

If you prefer to keep your equation logic in a dedicated file (good for more complex functions), create a .m file named lambda_equation.m with this code:

function y = lambda_equation(lambda, C, Tm, Ti, Lf)
    % Calculate the constant term from input parameters
    constant_term = C*(Tm - Ti)/(Lf*sqrt(pi));
    % Compute the value of your equation at the given lambda
    y = lambda * exp(lambda^2) * erfc(lambda) - constant_term;
end

Then call it from your main script/command window like this:

% Get user parameters
C = input('Enter value of C: ');
Tm = input('Enter value of Tm: ');
Ti = input('Enter value of Ti: ');
Lf = input('Enter value of Lf: ');

% Use an anonymous function wrapper to pass parameters to fzero
lambda_sol = fzero(@(x) lambda_equation(x, C, Tm, Ti, Lf), 1);

% Show the solution
disp(lambda_sol);

Critical Notes to Avoid Issues: #

• Initial Guess: fzero requires an initial guess (the second argument). For your equation, since λ·exp(λ²)·erfc(λ) ranges from 1/√π (as λ→0+) to 0 (as λ→∞), your constant term must be between 0 and 1/√π (~0.564) to have a real root. Choose an initial guess in this positive range (like 1) for best results.
• Check erfc: Make sure you're using MATLAB's built-in erfc function (complementary error function) — it should be available in all recent MATLAB versions.

内容的提问来源于stack exchange，提问作者Saurabh Das