作业帮matlab编程,是不是r的类型要定义?
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/11/08 18:08:55
matlab编程,是不是r的类型要定义?
r1=1/6*10^6;
w1=21*pi;
w=-10:1:10;
r=1i*(w1^2-w^2-1i*(w.*r1))./(2*w);
Error using ==> mpower
Inputs must be a scalar and a square matrix.
r1=1/6*10^6;
w1=21*pi;
w=-10:1:10;
r=1i*(w1^2-w^2-1i*(w.*r1))./(2*w);
Error using ==> mpower
Inputs must be a scalar and a square matrix.
=1i*(w1^2-w.^2-1i*(w*r1))./(2*w);
% 第一个w后加点,第二个w后不用加点
再问： 那这个呢？！g=1.42*10^5; %光子-声子耦合常数 Ip=1 ; As=1; r1=1/6*10^6 ; w1=21*pi ; w=19*pi:0.01*pi:23*pi; r=1i*(w1^2-w.*w-1i*(r1.*w))./(2*w); z1=1; T=5*10^(-6); t0=6*10^(-6); z0=2 ; Vg=2*10^8; z=0;a=1;b=-1; t=0:100; t=t*10^(-6); as=b*(conj(b)-(conj(b)-a)*exp(-conj(r).*(t-t0-(2*z0-z)/Vg))).*(heaviside(t-t0-(2*z0-z)/Vg)-heaviside(t-t0-T-(2*z0-z)/Vg));
再答： r 和t 长度都不一样，怎么能相乘
再问： 谢谢！如果每个r都要和全部的t相乘，应该怎么改呀？！需要定义一个数组吗？谢谢！
再答： 用循环啊
再问： 试过用循环，但是无奈本人不会用数组和矩阵~每个r有一个as，每个as有一百个数~。~！麻烦大侠可否帮忙改一下~谢谢！
再答： 大哥，到此为止吧，别再追问了 ------------------------------------------ clear;clc Ip=1;As=1;r1=1/6*10^6;w1=21*pi;w=19*pi:0.01*pi:23*pi; r=1i*(w1^2-w.*w-1i*(r1*w))./(2*w); z1=1;T=5*10^(-6);t0=6*10^(-6);z0=2;Vg=2*10^8; z=0;a=1;b=-1; t=0:100; t=t*10^(-6); for ii=1:length(r) as(ii,:)=b*(conj(b)-(conj(b)-a)*exp(-conj(r(ii)).*(t-t0-(2*z0-z)/Vg))).*(heaviside(t-t0-(2*z0-z)/Vg)-heaviside(t-t0-T-(2*z0-z)/Vg)); end as
% 第一个w后加点,第二个w后不用加点
再问： 那这个呢？！g=1.42*10^5; %光子-声子耦合常数 Ip=1 ; As=1; r1=1/6*10^6 ; w1=21*pi ; w=19*pi:0.01*pi:23*pi; r=1i*(w1^2-w.*w-1i*(r1.*w))./(2*w); z1=1; T=5*10^(-6); t0=6*10^(-6); z0=2 ; Vg=2*10^8; z=0;a=1;b=-1; t=0:100; t=t*10^(-6); as=b*(conj(b)-(conj(b)-a)*exp(-conj(r).*(t-t0-(2*z0-z)/Vg))).*(heaviside(t-t0-(2*z0-z)/Vg)-heaviside(t-t0-T-(2*z0-z)/Vg));
再答： r 和t 长度都不一样，怎么能相乘
再问： 谢谢！如果每个r都要和全部的t相乘，应该怎么改呀？！需要定义一个数组吗？谢谢！
再答： 用循环啊
再问： 试过用循环，但是无奈本人不会用数组和矩阵~每个r有一个as，每个as有一百个数~。~！麻烦大侠可否帮忙改一下~谢谢！
再答： 大哥，到此为止吧，别再追问了 ------------------------------------------ clear;clc Ip=1;As=1;r1=1/6*10^6;w1=21*pi;w=19*pi:0.01*pi:23*pi; r=1i*(w1^2-w.*w-1i*(r1*w))./(2*w); z1=1;T=5*10^(-6);t0=6*10^(-6);z0=2;Vg=2*10^8; z=0;a=1;b=-1; t=0:100; t=t*10^(-6); for ii=1:length(r) as(ii,:)=b*(conj(b)-(conj(b)-a)*exp(-conj(r(ii)).*(t-t0-(2*z0-z)/Vg))).*(heaviside(t-t0-(2*z0-z)/Vg)-heaviside(t-t0-T-(2*z0-z)/Vg)); end as