matlab二元函数驻点,求matlab高手帮忙写个关于对二元函数积分的程序 - 程序语言 - 小木虫 - 学术 科研 互动社区...

clear, clc;

%% build matrice

m1=[0.8+0.1i,0.5i;0.3i,0.3+0.4i];

m2=[0.1+0.05i,0.2i;

0,0.4+0.9i];

S=[1 0.6+0.4i];

%% differentiate symbolic M,

syms kx ky

M = m1 + m2*exp((kx -ky)*1i) + m2'*exp(-(kx-ky)*1i);

dmdkx = diff(M,kx); % 符号函数偏微分

%% get function F

F = S*dmdkx*S';

global ki; ki=0;

I =dblquad('f', -1,1,-2,2)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% function 单独放在一个文件(和Untitle文件同目录)

function F=f(kx,ky)

global ki;

ki=ki+1;

F= exp(- kx.*1i + ky.*1i)*(- 319/500 - (97*1i)/250) + exp(kx.*1i - ky.*1i)*(- 319/500 + (97*1i)/250);