MATLAB calcule la valeur de l'argument correspondant à la valeur extrême de la fonction complexe (y compris le nombre complexe, le produit conjugué complexe, pour la boucle)

Questions et réponses du rseoc 2021-12-30 03:14:49 阅读数:965

matlab calcule la valeur argument
Le phénomène rencontré par le problème et le contexte dans lequel il s'est produit

Pour obtenir l'argument correspondant au maximum d'une fonction complexexValeur de,EtxAvecwedgeLa relation entre,Voir ci - dessous

img

ArgumentsxChamp d'application:(0,0.2)
wedgeChamp d'application:(0.001,0.02)
La formule en cause est:

img

Codes liés au problème,Ne collez pas la capture d'écran
clc; clear all;close all; format long;%%%%%%%%%%%%%%%%Paramètres d'initialisation%%%%%%%%%%%%%%%syms x; %Définir les argumentsxw0=1; %Rayon de spot,Unité:mmz0=200; %La surface incidente est située à une distance de la taille,Unité:mmA=1; % Amplitude du Centre de la taille du Faisceau gaussien r=0.9; % Réflectivité du miroir de cavité k=A*(1-r); %Première foisx1 Amplitude de la lumière à la position du point lamda=632.8*10^(-6); % Longueur d'onde de crête ,Unité:mmL=75000*lamda; % Écart d'air du miroir ,Unité:mmf=pi*w0^2/lamda; % Paramètres confocaux du Faisceau gaussien beta=2*pi/lamda; % Constante de propagation du Faisceau gaussien %%%%%%%%%%%%%%%% Cycle de coin %%%%%%%%%%%%%%%wedge=(0.01:0.001:0.02); % Angle de coin de deux miroirs ,Unité:°Num=numel(wedge); % Donne le nombre de vecteurs %%%%%%%%%%%%%%% Calcul des paramètres pertinents %%%%%%%%%%%%%%%%%%for j=1:Num %CycleNumUne fois, Traverser tous les coins N=floor(pi./(4.*wedge(j))-0.5); %FloorArrondi vers le bas, C'est - à - dire pas plus de for m=0:N %S'il te plaît.m+1 Champ lumineux secondaire delt_xm=L./tand(wedge(j)).*(1./cosd((2*m+1).*wedge(j))-1./cosd(wedge(j))); % ..La distance parcourue par le Centre de la tache lumineuse sur le miroir d'éjection delt_zm=L./tand(wedge(j)).*tand((2*m+1).*wedge(j)); % La trajectoire du Centre de la tache xm=abs((x-delt_xm)*cosd((2*m+1).*wedge(j))); %pPoint à pointm+1 Distance verticale de l'axe de transmission du faisceau secondaire Secondaire sortant zm=z0+delt_zm+(x-delt_xm)*sind((2*m+1).*wedge(j)); % Distance parcourue par le Centre de la tache lumineuse R_zm=zm+f^2/zm; % Rayon de courbure de la surface isophasée w_zm=w0*sqrt(1+(zm/f)^2); %Rayon de spot Em(m+1)={@(x) k_r*r^m*w0/w_zm*exp(-xm^2/(w_zm^2))*exp(-(beta*zm-atan(zm/f))*i)*exp(-beta*xm^2/(2*R_zm)*i)}; endEt={@(x) sum(Em)}; %C'est exact.Em Somme de tous les éléments , Le résultat est l'inclusion d'arguments xAu pluriel deEt_c={@(x) conj(Et)}; %EtConjugaison complexe de, Le résultat est l'inclusion d'arguments xAu pluriel deIntensity_total={@(x) Et*Et_c}; %EmAvecEtMultiplication, Le résultat est xNombre réel [xmax]=fminbnd(inline(-Intensity_total(x)),0,0.2); %S'il te plaît.Intensity_totalAu maximum,xValeur dex_p(j)=xmin; %Pour chaquewedge,Il y en a unxValeur de, Stocké dans un vecteur x_pMoyenne end
Résultats de l'opération et contenu du rapport d'erreur

img

img

Mes idées de solution et les méthodes que j'ai essayées

Ne cherche pas d'inconnu x Fonctionnement de la fonction

Le résultat que je voulais

Allez.xLes vecteurs qui composent,xAvecwedgeDiagramme pour




Adopter les réponses:

Résolu par lui - même



Autres réponses2:

Bonjour,Je suis un petit assistant.,Je suis vraiment désolé,Toutes les questions que vous avez posées cette fois,Le Groupe d'experts techniques a expiré sans réponse


Nombre de fois où vous devez répondre aux questions retenues pour cette question,Il y aura des questions et des réponsesVIPCarte d'expérience(1Une chance de répondre à une question.、Achat de livres physiques dans les centres commerciaux95Réduction)Le formulaire est envoyé à votre compte.


Parce qu'il faut répondre aux questions.VIPCarte d'expérience valide seulement1Oh, mon Dieu.,Quand vous en avez besoin【Message privé】Contactez - moi.,Je vais le refaire pour toi..

版权声明:本文为[Questions et réponses du rseoc]所创,转载请带上原文链接,感谢。 https://primo.wiki/2021/12/202112231148286411.html